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Abstract  (corn.) 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode,  and  the  perfor¬ 
mance  has  been  very  satisfactory.  Since  October  1991,  a  new  version  of  the  IMS  that 
accepts  data  from  an  arbitrary  number  of  arrays  and  single  3 -component  stations  has  been 
operated. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  97.8%  as  compared  to  96.7%  for  the  previous  reporting  period.  A 
total  of  1733  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  The  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  tested  successfully  for  acquiring  such  data 
by  AFTAC.  Processing  of  requests  for  full  NORSAR  and  regional  array  data  on  magnetic 
tapes  has  progressed  according  to  established  schedules.  There  have  been  no  modifica¬ 
tions  made  to  the  NORSAR  data  acquisition  system. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS,  FINESA  and  GERESS  data  have  been  conducted 
throughout  the  period.  Data  from  two  experimental  small-aperture  arrays  at  sites  in  Spits¬ 
bergen  and  Apatity,  Kola  Peninsula,  have  been  recorded  and  processed  in  an  experimental 
mode.  Monthly  processing  statistics  for  the  arrays  as  well  as  results  of  the  IMS  analysis 
for  the  reporting  period  are  given. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  In  addition,  the  mainte¬ 
nance  center  has  been  involved  with  occasional  maintenance  of  equipment  for  FINESA. 
Other  activities  have  involved  testing  of  the  NORSAR  communications  systems  and  work 
in  connection  with  the  experimental  small-aperture  arrays  in  Spitsbergen  and  the  Kola 
Peninsula. 

Starting  1  October  1991,  an  effort  began  to  carry  out  a  complete  technical  refurbishment 
of  the  NORSAR  array.  This  project  is  funded  jointly  by  AFTAC,  ARPA  and  NTNF.  Dur¬ 
ing  the  reporting  period,  efforts  have  focused  upon  evaluation  and  laboratory  testing  of 
technical  options  for  field  instrumentation,  in  particular  state-of-the-art  A/D  converters, 
data  acquisition  and  synchronization  devices.  During  the  reporting  period,  we  have  also 
been  testing  several  such  systems  under  realistic  operating  conditions  in  the  field.  When 
these  studies  have  been  completed,  a  recommendation  for  a  system  to  be  installed  will  be 
presented  to  the  funding  agencies. 

Summaries  of  six  scientific  contributions  are  presented  in  Chapter  7  of  this  report 

Section  7. 1-7.3  describe  the  basic  principles  and  give  examples  of  application  of  a  promis¬ 
ing  new  approach  to  automatic  analysis  of  events  from  sites  with  recurring  seismic  activ¬ 
ity.  This  method,  which  we  have  called  “intelligent  post-processing  of  seismic  events”  is 
particularly  well  suited  to  implement  as  a  post-event  processing  technique  in  the  Intelli¬ 
gent  Monitoring  System. 
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Section  7.1  describes  the  basic  principles  of  this  approach.  From  experience  with  analyst 
review  of  events  automatically  defined  by  the  Intelligent  Monitoring  System  (IMS),  we 
have  realized  that  the  quality  of  the  automatic  event  locations  can  be  significantly 
improved  if  the  event  intervals  are  reprocessed  with  signal  processing  parameters  tuned  to 
phases  from  events  in  the  given  region.  The  tuned  processing  parameters  are  obtained 
from  off-line  analysis  of  events  located  in  the  region  of  interest  The  primary  goal  of  such 
intelligent  post-processing  is  to  provide  event  definitions  of  a  quality  that  minimizes  the 
need  for  subsequent  analysis  and  review. 

The  first  step  in  this  post-processing  is  to  subdivide  the  area  to  be  monitored  in  order  to 
identify  sites  of  interest  Clearly,  calibration  will  be  easiest  and  potential  savings  in  man¬ 
power  are  largest  for  areas  with  high,  recurring  seismic  activity.  We  have  identified  8  min¬ 
ing  sites  in  Fennoscandia/W.  Russia/Estonia  and  noted  that  65.6  per  cent  of  the  events  of 
Ml  >  2.0  in  this  region  can  be  associated  with  one  of  these  sites.  This  result  is  based  on  1 
1/2  years  of  data. 

Section  12  describes  the  second  step  in  the  post-processing  approach.  This  second  step  is 
to  refine  the  phase  arrival  time  and  azimuth  estimates  using  frequency  filters  and  process¬ 
ing  parameters  that  are  tuned  to  the  initial  event  location  provided  by  the  IMS.  We  have 
studied,  as  a  first  example,  a  set  of  58  events  from  the  Khibiny  Massif  in  the  Kola  penin¬ 
sula.  Very  accurate  locations  of  these  events  has  been  provided  by  the  Kola  Regional  Seis- 
mological  Centre.  Our  refinement  of  phase  arrival  times  for  these  events,  as  recorded  in 
Apatity  and  at  ARCESS,  have  given  quite  remarkable  results.  Using  the  autoregressive 
likelihood  technique  of  Pisarenko  et  al  ( 1987)  we  have  been  able  to  estimate  onset  times 
automatically  to  an  accuracy  (standard  deviation)  of  about  0.05  s  for  P-phases  and  0.15- 
0.20  s  for  S-phases.  These  accuracies  are  as  good  as  for  analyst  picks,  and  are  consider¬ 
ably  better  than  in  the  current  SigPro  analysis. 

Section  7.3  describes  the  third  step  in  the  post-processing  approach.  This  third  step  in  the 
post-processing  is  a  recomputation  of  the  location  estimate,  using  the  refined  arrival  times 
and  broad-band  azimuths  associated  with  a  fixed  frequency  band  (tied  to  the  initial  IMS 
location).  Again  using  the  set  of  58  known  Khibiny  Massif  events  as  a  reference,  we  have 
investigated  the  accuracy  of  the  LocSat  location  procedure  using  varying  amounts  of 
regional  corrections.  Using  a  depth  constrained  to  0,  we  find  that  the  epicenter  can  be  esti¬ 
mated  to  2-3  km  accuracy  (median  value),  with  “worst  case”  error  less  than  10  km,  even 
without  any  regional  corrections.  It  is  noted  that  the  method  needs  to  be  tested  for  addi¬ 
tional  target  sites  before  firm  conclusions  can  be  drawn. 

Section  7.4  describes  a  four-month  monitoring  experiment,  using  the  Continuous  Thresh¬ 
old  Monitoring  technique  applied  to  the  Novaya  Zemlya  test  site.  Data  from  the  NORESS, 
ARCESS  and  FINES  A  regional  arrays  were  used.  Starting  1  December  1992,  we  have 
been  compiling  daily  statistics  of  all  peaks  on  the  threshold  diagram  exceeding  m^,  =  2.75, 
and  associating  these  peaks  to  regional  or  teleseismic  events  whenever  possible.  In  addi¬ 
tion,  we  are  analyzing  smaller  peaks  (below  mb  =  2.75)  that  can  possibly  be  associated 
with  Novaya  Zemlya  epicenters. 

We  have  conducted  extensive  analysis  of  the  m^  =  2.5  seismic  event  at  Novaya  Zemlya  at 
09.29.25  GMT  on  31  December  1992.  This  event  was  detected  by  ARCESS  (P  and  S), 
Spitsbergen  (P  and  S),  NORESS  (P)  and  Apatity  (S).  In  addition,  the  Kola  Regional  Seis- 
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mological  Centre  provided  readings  for  the  station  Amderma  (Pg  and  Sn)  made  from  ana¬ 
log  recordings.  Our  results  indicate  that  the  epicenter  was  slightly  to  the  north  of  die 
northern  Novaya  Zemlya  test  site.  However,  there  are  some  uncertainties  in  the  travel  time 
tables,  and  an  on-site  location  cannot  entirely  be  ruled  out 

Section  7.5  contains  some  initial  processing  results  from  the  Spitsbergen  array.  Prelimi¬ 
nary  analysis  of  data  from  this  new  small-aperture  (1km)  array  on  the  Arctic  island  of 
Spitsbergen  has  revealed  a  large  variation  in  signal  and  noise  amplitudes  across  the  array.. 
Further  investigations  are  needed  in  order  to  determine  the  reason  for  this  variation,  and  it 
may  prove  necessary  to  redeploy  one  or  more  of  the  array  sensors.  Site  B2  in  this  array 
exhibits  a  particularly  favorable  signal  focusing  for  nearly  all  arrival  azimuths.  Initial 
studies  of  the  noise  characteristics  show  that  the  noise  is  quite  stable  over  time,  especially 
at  high  frequencies.  Also,  the  noise  amplitudes  at  high  frequencies  seem  to  be  similar 
across  the  array,  in  contrast  to  the  large  variation  seen  at  low  frequencies. 

Section  7.6  contains  results  from  an  evaluation  of  the  most  recent  version  of  the  Intelligent 
Monitoring  System  (IMS).  In  order  to  assess  the  performance  of  IMS,  the  automatic 
results  for  a  one-week  test  period  were  carefully  and  thoroughly  reviewed.  The  conclusion 
is  that  the  overall  performance  of  IMS  has  improved  considerably  from  previous  versions 
and  can  now  be  rated  as  very  satisfactory,  but  also  that  there  is  still  room  for  some 
improvements.  Nearly  60%  of  the  events  automatically  declared  by  IMS  were  judged  to 
be  acceptable  in  the  sense  that  they  correspond  to  real  seismic  events  for  which  only  rela¬ 
tively  minor  modifications  were  made  during  the  analyst  review  process.  The  majority  of 
the  remaining  20%  were  judged  to  be  false  events,  and  some  modifications  to  IMS  are  pro¬ 
posed  that  may  help  in  reducing  the  number  of  such  events.  Using  various  reference  bulle¬ 
tins,  it  was  found  that  some  events  were  missed  by  IMS,  even  if  data  were  available  to 
IMS  that  should  allow  it  to  associate  phases  and  form  events.  An  evaluation  of  the  manual 
corrections  made  to  the  automatically  determined  onset  times  shows  that  there  is  a  poten¬ 
tial  for  improving  these  by,  e.g.,  implementing  in  IMS  new  procedures  for  onset  time  esti¬ 
mation  like  those  described  in  section  7.2. 
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1  Summary 

This  Semiannual  Technical  Summary  describes  the  operation,  maintenance  and  research 
activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Norwegian  Regional  Seismic 
Array  (NORESS,  the  Arctic  Regional  Seismic  Array  (ARCESS)  and  the  experimental 
Spitsbergen  regional  array  for  the  period  1  October  1992  -31  March  1993.  Statistics  are 
also  presented  for  additional  seismic  stations,  which  through  cooperative  agreements  with 
institutions  in  the  host  countries  provide  continuous  data  to  the  NORSAR  Data  Processing 
Center  (NPDC).  These  stations  comprise  the  Finnish  Experimental  Seismic  Array 
(FINESA),  the  German  Experimental  Seismic  Array  (GERESS),  and  an  experimental 
regional  seismic  array  in  Apatity.  Russia. 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode,  and  the  perfor¬ 
mance  has  been  very  satisfactory.  Since  October  1991,  a  new  version  of  the  IMS  that 
accepts  data  from  an  arbitrary  number  of  arrays  and  single  3-component  stations  has  been 
operated. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  97.8%  as  compared  to  96.7%  for  the  previous  reporting  period.  A 
total  of  1733  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  The  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  tested  successfully  for  acquiring  such  data 
by  AFTAC.  Processing  of  requests  for  full  NORSAR  and  regional  array  data  on  magnetic 
tapes  has  progressed  according  to  established  schedules.  There  have  been  no  modifica¬ 
tions  made  to  the  NORSAR  data  acquisition  system. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS,  FINESA  and  GERESS  data  have  been  conducted 
throughout  the  period.  Data  from  two  experimental  small-aperture  arrays  at  sites  in  Spits¬ 
bergen  and  Apatity,  Kola  Peninsula,  have  been  recorded  and  processed  in  an  experimental 
mode.  Monthly  processing  statistics  for  the  arrays  as  well  as  results  of  the  IMS  analysis 
for  the  reporting  period  are  given. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  In  addition,  the  mainte¬ 
nance  center  has  been  involved  with  occasional  maintenance  of  equipment  for  FINESA. 
Other  activities  have  involved  testing  of  the  NORSAR  communications  systems  and  work 
in  connection  with  the  experimental  small-aperture  arrays  in  Spitsbergen  and  the  Kola 
Peninsula. 

Starting  1  October  1991,  an  effort  began  to  carry  out  a  complete  technical  refurbishment 
of  the  NORSAR  array.  This  project  is  funded  jointly  by  AFTAC,  ARPA  and  NTNF.  Dur¬ 
ing  the  repeating  period,  efforts  have  focused  upon  evaluation  and  laboratory  testing  of 
technical  options  for  field  instrumentation,  in  particular  state-of-the-art  A/D  converters, 
data  acquisition  and  synchronization  devices.  During  the  reporting  period,  we  have  also 
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been  testing  several  such  systems  under  realistic  operating  conditions  in  the  field.  When 
these  studies  have  been  completed,  a  recommendation  for  a  system  to  be  installed  will  be 
presented  to  the  funding  agencies. 

Summaries  of  six  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

Section  7. 1-7.3  describe  the  basic  principles  and  give  examples  of  application  of  a  promis¬ 
ing  new  approach  to  automatic  analysis  of  events  from  sites  with  recurring  seismic  activ¬ 
ity.  This  method,  which  we  have  called  “intelligent  post-processing  of  seismic  events”  is 
particularly  well  suited  to  implement  as  a  post-event  processing  technique  in  the  Intelli¬ 
gent  Monitoring  System. 

Section  7.1  describes  the  basic  principles  of  this  approach.  From  experience  with  analyst 
review  of  events  automatically  defined  by  the  Intelligent  Monitoring  System  (IMS),  we 
have  realized  that  the  quality  of  the  automatic  event  locations  can  be  significantly 
improved  if  the  event  intervals  are  reprocessed  with  signal  processing  parameters  tuned  to 
phases  from  events  in  the  given  region.  The  tuned  processing  parameters  are  obtained 
from  off-line  analysis  of  events  located  in  the  region  of  interest.  The  primary  goal  of  such 
intelligent  post-processing  is  to  provide  event  definitions  of  a  quality  that  minimizes  the 
need  for  subsequent  analysis  and  review. 

The  first  step  in  this  post-processing  is  to  subdivide  the  area  to  be  monitored  in  order  to 
identify  sites  of  interest.  Clearly,  calibration  will  be  easiest  and  potential  savings  in  man¬ 
power  are  largest  for  areas  with  high,  recurring  seismic  activity.  We  have  identified  8  min¬ 
ing  sites  in  Fennoscandia/W.  Russia/Estonia  and  noted  that  65.6  per  cent  of  the  events  of 
Ml  >  2.0  in  this  region  can  be  associated  with  one  of  these  sites.  This  result  is  based  on  1 
1/2  years  of  data. 

Section  7.2  describes  the  second  step  in  the  post-processing  approach.  This  second  step  is 
to  refine  the  phase  arrival  time  and  azimuth  estimates  using  frequency  filters  and  process¬ 
ing  parameters  that  are  tuned  to  the  initial  event  location  provided  by  the  IMS.  We  have 
studied,  as  a  first  example,  a  set  of  58  events  from  the  Khibiny  Massif  in  the  Kola  penin¬ 
sula.  Very  accurate  locations  of  these  events  has  been  provided  by  the  Kola  Regional  Seis- 
mological  Centre.  Our  refinement  of  phase  arrival  times  for  these  events,  as  recorded  in 
Apatity  and  at  ARCESS,  have  given  quite  remarkable  results.  Using  the  autoregressive 
likelihood  technique  of  Pisarenko  et  al  (1987)  we  have  been  able  to  estimate  onset  times 
automatically  to  an  accuracy  (standard  deviation)  of  about  0.05  s  for  P-phases  and  0. 15- 
0.20  s  for  S-phases.  These  accuracies  are  as  good  as  for  analyst  picks,  and  are  consider¬ 
ably  better  than  in  the  current  SigPro  analysis. 

Section  7.3  describes  the  third  step  in  the  post-processing  approach.  This  third  step  in  the 
post-processing  is  a  recomputation  of  the  location  estimate,  using  the  refined  arrival  times 
and  broad-band  azimuths  associated  with  a  fixed  frequency  band  (tied  to  the  initial  IMS 
location).  Again  using  the  set  of  58  known  Khibiny  Massif  events  as  a  reference,  we  have 
investigated  the  accuracy  of  the  LocSat  location  procedure  using  varying  amounts  of 
regional  corrections.  Using  a  depth  constrained  to  0,  we  find  that  the  epicenter  can  be  esti¬ 
mated  to  2-3  km  accuracy  (median  value),  with  “worst  case”  error  less  than  10  km,  even 
without  any  regional  corrections.  It  is  noted  that  the  method  needs  to  be  tested  for  addi¬ 
tional  target  sites  before  firm  conclusions  can  be  drawn. 
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Section  7.4  describes  a  four-month  monitoring  experiment,  using  the  Continuous  Thresh¬ 
old  Monitoring  technique  applied  to  the  Novaya  Zemlya  test  site.  Data  from  the  NORESS, 
ARCESS  and  FINES  A  regional  arrays  were  used.  Starting  1  December  1992,  we  have 
been  compiling  daily  statistics  of  all  peaks  on  the  threshold  diagram  exceeding  mb  =  2.75, 
and  associating  these  peaks  to  regional  or  teleseismic  events  whenever  possible.  In  addi¬ 
tion,  we  are  analyzing  smaller  peaks  (below  mb  =  2.75)  that  can  possibly  be  associated 
with  Novaya  Zemlya  epicenters. 

We  have  conducted  extensive  analysis  of  the  mb  =  2.5  seismic  event  at  Novaya  Zemlya  at 
09.29.25  GMT  on  31  December  1992.  This  event  was  detected  by  ARCESS  (P  and  S), 
Spitsbergen  (P  and  S),  NORESS  (P)  and  Apatity  (S).  In  addition,  the  Kola  Regional  Seis- 
mologicai  Centre  provided  readings  for  the  station  Amderma  (Pg  and  Sn)  made  from  ana¬ 
log  recordings.  Our  results  indicate  that  the  epicenter  was  slightly  to  the  north  of  the 
northern  Novaya  Zemlya  test  site.  However,  there  are  some  uncertainties  in  the  travel  time 
tables,  and  an  on-site  location  cannot  entirely  be  ruled  out. 

Section  7.5  contains  some  initial  processing  results  from  the  Spitsbergen  array.  Prelimi¬ 
nary  analysis  of  data  from  this  new  small-aperture  (1km)  array  on  the  Arctic  island  of 
Spitsbergen  has  revealed  a  large  variation  in  signal  and  noise  amplitudes  across  the  array. 
Further  investigations  are  needed  in  order  to  determine  the  reason  for  this  variation,  and  it 
may  prove  necessary  to  redeploy  one  or  more  of  the  array  sensors.  Site  B2  in  this  array 
exhibits  a  particularly  favorable  signal  focusing  for  nearly  all  arrival  azimuths.  Initial 
studies  of  the  noise  characteristics  show  that  the  noise  is  quite  stable  over  time,  especially 
at  high  frequencies.  Also,  the  noise  amplitudes  at  high  frequencies  seem  to  be  similar 
across  the  array,  in  contrast  to  the  large  variation  seen  at  low  frequencies. 

Section  7.6  contains  results  from  an  evaluation  of  the  most  recent  version  of  the  Intelligent 
Monitoring  System  (IMS).  In  order  to  assess  the  performance  of  IMS,  the  automatic 
results  for  a  one-week  test  period  were  carefully  and  thoroughly  reviewed.  The  conclusion 
is  that  the  overall  performance  of  IMS  has  improved  considerably  from  previous  versions 
and  can  now  be  rated  as  very  satisfactory,  but  also  that  there  is  still  room  for  some 
improvements.  Nearly  80%  of  the  events  automatically  declared  by  IMS  were  judged  to 
be  acceptable  in  the  sense  that  they  correspond  to  real  seismic  events  for  which  only  rela¬ 
tively  minor  modifications  were  made  during  the  analyst  review  process.  The  majority  of 
the  remaining  20%  were  judged  to  be  false  events,  and  some  modifications  to  IMS  are  pro¬ 
posed  that  may  help  in  reducing  the  number  of  such  events.  Using  various  reference  bulle¬ 
tins,  it  was  found  that  some  events  were  missed  by  IMS,  even  if  data  were  available  to 
IMS  that  should  allow  it  to  associate  phases  and  form  events.  An  evaluation  of  the  manual 
corrections  made  to  the  automatically  determined  onset  times  shows  that  there  is  a  poten¬ 
tial  for  improving  these  by,  e.g.,  implementing  in  IMS  new  procedures  for  onset  time  esti¬ 
mation  like  those  described  in  section  7.2. 
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2  NORSAR  Operation 

2.1  Detection  Processor  (DP)  operation 

There  have  been  58  breaks  in  the  otherwise  continuous  operation  of  the  NORSAR  online 
system  within  the  current  6-month  reporting  interval.  The  uptime  percentage  for  the 
period  is  97.8  as  compared  to  96.7  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  downtime  for  the 
days  between  1  October  1992  and  31  March  1993.  The  monthly  recording  times  and  per¬ 
centages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  Hardware  failure  3 

b)  Stops  related  to  program  work  or  error  2 

c)  Hardware  maintenance  stops  0 

d)  Power  jumps  and  breaks  0 

e)  TOD  error  correction  5 

f)  Communication  lines  48 


The  total  downtime  for  the  period  was  1 19  hours  and  3  minutes.  The  mean-time-between- 
failures  (MTBF)  was  3.6  days,  as  compared  to  2.0  for  the  previous  period. 
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Fig.  2.1.1.  Detection  Processor  uptime  for  October  (top),  November  (middle)  and 
December  (bottom)  1992. 
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Date 

Time 

Cause 

12  Oct 

0939 

-  1140 

Software  work 

21  Oct 

0720 

-  1315 

Line  failure 

02  Nov 

03  Nov 

1605 

-  0555 

Line  failure 

07  Dec 

08  Dec 

1801 

-  0606 

Line  failure 

21  Dec 

0423 

-0802 

Line  failure 

01  Jan 

1007 

-  1131 

Software  woric 

31  Jan 

0906 

-  1448 

Hardware  failure 

31  Jan 

1808 

- 

Hardware  failure 

01  Feb 

- 

Hardware  failure 

02  Feb 

-  1455 

Hardware  failure 

03  Mar 
04  Mar 

0759 

-0747 

Line  failure 

15  Mar 

1248 

-  1314 

Hardware  failure 

Table  2.1.1.  The  major  downtimes  in  the  period  1  October  1992  -  31  March  1993. 


Month 

DP  Uptime 
Hours 

DP  Uptime 
% 

No.  of 

DP  Breaks 

No.  of 
Days  with 
BrcaKS 

DPMTBF* 

(days) 

Oct  92 

735.24 

98.84 

11 

9 

2.6 

Nov  92 

705.09 

97.94 

14 

11 

2.0 

Dec  92 

728.10 

97.87 

4 

4 

6.1 

Jan  93 

723.57 

97.31 

8 

7 

3.4 

Feb  93 

724.51 

97.43 

4 

4 

6.0 

Mar  93 

700.45 

97.34 

17 

13 

1.6 

97.79 

58 

48 

3.6 

♦Mean-time-between-failures  -  total  uptime/no.  of  up  intervals. 


Table  2.1.2.  Online  system  performance,  1  October  1992  -  31  March  1993. 
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2.2  Array  communications 

General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  throughout  the  report¬ 
ing  period. 

The  most  prominent  events  which  have  affected  the  different  systems  have  been:  loss  of 
synch,  bad  cable,  spikes  on  data,  PCM  irregularities.  Digital  Access  Cross  Connection 
System  (DACCS),  SLEM,  line  outage  and  2701  data  adapter. 

Detailed  Summary 

October  (weeks  40-44),  28.9-1.11.92 

NTA/Lillestrom  replaced  the  carrier  frequency  system  pertaining  to  01 A-06C  (-02C),  with 
pulse  code  modulation  equipment  (PCM)  between  Lilies trom  and  Hamar  1  Oct  93.  An 
attempt  to  replace  the  02C  system  was  also  made,  but  without  success. 

01A,  which  had  been  down  since  18  June,  was  bade  in  operation  23  Oct  with  satisfactory 
performance. 

8  Oct  01B  started  with  deteriorated  data  (spikes)  and  excessive  error  rate,  after  the  original 
SLEM  was  reinstalled.  Between  9  and  12  Oct  the  subarray  was  out  of  operation.  During  a 
subarray  visit  13  Oct  a  SLEM  failure  occurred  and  after  this  the  system  was  down 
throughout  October. 

02C  continued  to  cause  us  trouble  after  NTA’s  first  attempt  to  replace  the  carrier  system 
with  PCM  equipment  (1  Oct).  Although  it  was  (according  to  NTA/Lillehammer)  changed 
back  to  the  original  system  the  same  day.  Between  2  and  7  Oct  the  subarray  was  affected 
by  different  irregularities  related  to  the  communication  path,  and  was  more  or  less  out  of 
operation  until  12  Oct,  when  it  was  demasked.  On  28  Oct  02C  was  again  affected  in  con¬ 
nection  with  NTA  work  in  Oslo. 

Other  subarrays  which  were  affected  in  October  were  03C  and  06C;  03C  between  4  and  6 
Oct  and  06C  between  22  and  26  Oct. 

All  systems  were  affected  21  Oct  between  0720  and  1316  GMT. 

Average  outages  for  October,  individual  weeks  40-44 


Week  40  (-01A.02C)  0.003% 

Week  41  (-01  A,02C)  0.001% 

Week  42  (-01  A.01B.02C)  0.002% 

Week  43  (N/A) 

Week  44-01B,02C,06C )  0.003% 
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November  (weeks  45-48),  2-29.11.92 

01B  resumed  operation  11  Nov  (week  46)  after  having  been  down  since  8  Oct,  due  to  a 
SLEM  failure. 

In  the  same  period  02C  and  04C  were  also  affected.  02C  had  synch  problems,  and  the  04C 
irregularity  was  caused  by  NTA/Hamar. 

02B  was  affected  between  12  and  16  Nov  and  from  17  through  30  Nov  due  to  a  bad  com¬ 
munications  cable.  Also  the  02C  synch  problem  continued  week  46  (39.5%  down),  and 
week  47  (17.4%  down). 

06C  was  affected  week  48  (38.1%  down),  probably  by  a  line  outage.  A  Modcomp  restart 
caused  06C  reinitialization. 

Average  outages  in  November,  individual  weeks: 


Week  45  (-01B,02C,04C1)  0.002% 

Week  46  (-01B,02B,02C1)  0.0006% 

Week  47  (-02B,02Q  0.003% 

Week  48(-02B,06C)  0.002% 


December  (weeks  49-53),  30.11.92-3.1.93 

02B  went  down  week  46  due  to  a  bad  communications  cable  located  in  a  swampy  area, 
with  snow  depths  exceeding  1.5  meters.  A  cable  repair  was  therefore  difficult  according  to 
NTA/Hamar.  They  suggested,  however,  a  start  date  mid-January  93.  21  December,  after  1 
a.m.  all  systems  were  affected  4  hours  caused  by  an  NTA  system  irregularity. 

Average  outage  in  December,  individual  weeks: 


Week  49  (-02B)  0.0008% 

Week  50  (-02B)  0.0003% 

Week  51  (-02B)  0.001% 

Week  52  (-02B)  :  2.380% 

Week  53  (-02B)  0.0008% 


January  (weeks  1-4),  4-31.1.93 

NTA/Hamar  finally  finished  the  repair  of  the  02B  communications  cable  on  29  Jan  93.  The 
subarray  had  been  out  ■  f  operation  since  week  46/92. 

After  a  few  days’  observation  of  the  performance,  the  subarray  communications  were 
accepted.  Communications  errors  were  not  discovered  between  29  Jan  0945  GMT  and  31 
Jan  approx  1800  GMT,  which  was  the  approximate  hour  that  the  communication  adapter 
2701  stopped,  and  the  data  transfer  between  the  Modcomp  and  the  IBM  4381  machine 
was  suspended.  The  2701  remained  down  throughout  31  Jan. 

Average  outages  in  January,  individual  weeks: 
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Week  1  (-Q2B)  0.002% 

Week  2  (-02B  0.002% 

Week  3  (-02B)  0.001% 

Week  4  (-,02B)  0.001% 


February  (weeks  5-8),  1-28193 

The  2701  communications  adapter  which  failed  31  Jan  93  approx.  1800  GMT  was 
repaired  3  Feb,  and  the  NORS AR  system  started  again  1 433  GMT.  After  restart  all  sys¬ 
tems  performed  satisfactorily. 

Average  outages  in  February,  individual  weeks: 

Week  3  0.002% 

Week  6  0.002% 

Week  7  0.002% 

Week  8  0.002% 

March  (weeks  9-12),  1-283.92 

3  March  a  DACCS  (Data  Access  Cross  Connection  System)  failed  at  the  NTA/Lillestrem 
premises  and  affected  all  NORSAR  subarray  communications  systems  (01 A-06C). 
According  to  NTA/Lillestrom  all  the  systems  (-02C)  were  in  operation  approx  2100  hrs 
GMT.  On  4  March  at  0544  GMT,  an  attempt  to  start  NORSAR  was  made,  but  without  suc¬ 
cess.  The  TOD  (Time  Code  Generator),  which  was  replaced  3  March,  was  reinstalled. 
0747  GMT  a  new  ONLINE  start  was  successful,  but  02C  was  still  inoperative. 

During  the  period  misalignment  of  received  data  and  TOD  generator  clock  pulses  have 
resulted  in  several  Modcomp  restarts.  4  March  02C,  which  remained  down  after  the  NTA7 
DACCS  failure  in  Lillestrom,  was  loop-tested  between  Kjeller  and  Sjusjoen.  No  errors 
were  observed. 

15  March  a  NORSAR  ONLINE  stop  occurred  at  1248  hrs  GMT,  up  25  min  42  sec  later. 

16  March  a  new  stop  occurred  at  1256  GMT,  and  the  system  was  down  for  about  5  min  42 
seconds. 

18  March  NTA/Hamar  carried  out  scheduled  work  on  the  02C  PCM  system  between  the 
Nes  peninsula  and  Sjusjoen. 

22  March  deteriorated  data  (spikes)  were  observed  on  01  A,  02B  and  02C.  At  the  same 
time  a  2-second  misalignment  of  received  NORSAR  data  and  TOD  clock  pulses  were 
detected.  In  order  to  find  out  if  the  seriously  deteriorated  data  had  created  a  situation  in  the 
Modcomp  which  resulted  in  the  misalignment,  the  three  subarrays  were  omitted.  In 
response  to  a  request,  NTA/Hamar  informed  us  that  an  attenuator  had  been  removed  in 
connection  with  changing  the  path  between  Loten  and  06C  from  PCM  cable  to  radiolink. 
After  having  reinstalled  the  attenuator,  NORSAR  performance  was  significantly 
improved. 
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Avenge  outages  in  March,  individual  weeks: 


Week  9  (-020  0.002% 

Week  10  (-020  0.001% 

Week  11  (all)  0.022% 

Week  12  (all)  0.522% 


OA.  Hansen 


Subarrays 

Oct  (5) 

Nov  (4) 

Dec  (5) 

Jan  (4) 

Feb  (4) 

Mar  (4 

Average 

28.9- 

1.11.93 

2-29.11.93 

30.11.92- 

3.1.93 

4-31.1.93 

1-28.2.93 

1-28.3.93 

1/2  year 

01A 

ooio1) 

0.004 

0.001 13) 

0.008 

0.001 

0.225 

0.047 

01B 

0.0302* 

0.00058) 

0.0009 M) 

0.0007 

0.0006 

0.002 

0.006 

02B 

0.00073) 

0.0019> 

100.0N/A 

92.755wa 

0.002 

0.329 

0.00320) 

02C 

0.0054) 

0.002,O) 

0.001 15) 

0.002 

0.002 

0.010,9) 

0.004 

03C 

0.0023) 

0.0004 

0.001,6) 

0.001 

0.0004 

0.011 

0.003 

04C 

0.0016) 

0.0001n) 

0.001 17) 

0.001 

0.004 

0.009 

0.003 

06C 

O.OtE7* 

0.00212) 

0.00318) 

0.003 

0.002 

0.369 

0.063 

AVER 

0.007 

0.001 

0.00121) 

0.01622* 

0.002 

0.136 

0.018 

Figures  representing  error  rate  (in  per  cent)  followed  by  number  1),  2),  etc.,  are  related  to  legend  below. 


Table  2.2.1.  Communications  performance.  The  numbers  represent  error  rates  in  per  cent 
based  on  total  transmitted  frames/week  (28  Sep  92  -  28  Mar  93). 


1).4) 

one  week  (43/44) 

9).  10) 

one  week  (45/48) 

2),  8),  19) 

average  2  weeks  (40,41/47,48)  (11,12/93) 

7),  11),12) 

average  3  weeks  (40,41,42/46,47,48/45,46,47) 

3),  5).  6) 

average  4  weeks  (40,41,42,44) 

13),  14),  15).  16) 

average  4  weeks  (49,50,51,53) 

17),  18) 

average  4  weeks  (49,50,51,53) 

20) 

average  4  months  (Oct,  Nov  92/  Feb,  Mar  93) 

21) 

average  6  subarrays  (01  A.02B,  02C-06Q 

22) 

average  6  subarrays  (01  A.02B,  02C-06Q 
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23  NORSAR  Event  Detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor  operation  are 
given.  The  table  lists  the  total  number  of  detections  (DPX)  triggered  by  the  on-line  detec¬ 
tor,  the  total  number  of  detections  processed  by  the  automatic  event  processor  (EPX)  and 
the  total  number  of  events  accepted  after  analyst  review  (teleseismic  phases,  core  phases 
and  total). 


Total 

DPX 

Total 

EPX 

Accepted  events 

Sum 

Daily 

P-phases 

Core 

Phases 

Oct  92 

13450 

1812 

241 

131 

372 

12.0 

Nov  92 
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1324 

189 

50 

239 

8.0 

Dec  92 
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1512 
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61 

250 

8.1 

Jan  93 

11324 

1265 

148 

74 

222 

myj 

Feb  93 

11450 

1241 

225 

61 

286 

10.2 

Mar  93 

13000 

1505 

257 

107 

364 

11.7 

1249 

484 

1733 

9.5 

Table  23.1.  Detection  and  Event  Processor  statistics,  1  October  1992  -  31  March  1993. 


NORSAR  Detections 

The  number  of  detections  (phases)  reported  by  the  NORSAR  detector  during  day  275, 
1992,  through  day  090, 1993,  was  74,237,  giving  an  average  of  412  detections  per  pro¬ 
cessed  day  (180  days  processed).  Table  2.3.2  shows  daily  and  hourly  distribution  of  detec¬ 
tions  for  NORSAR. 


T.  Schftyen 


. OPX  Hourly  distribution  of  detections 


NCMSAiaei  lap.  2-92/93 


May  1993 


rH09m*,w’»®P'-®ff‘O'"i0903*,m®r-«a'©rH09fn*,m®r-®0*orHrH0«r’>^u"»®r*“®tf*OrH09f'>*‘m®r>.a>0N©rH09oi'*>m 

oooooooooHHHriHHHHHF4(innnoinoin(N(innoooooooooHHHHHHHHHHr«noinr«n 

*i*)4i*i*»*i*>*j***i*>*j*)*»*j*j*j*»*>*>y*j*i***i*iii*)*<*)*)>>>>>>>>>>>>>>>>>>>>>>>>> 

VSJ?V5IS5Syj/25Siyi!V5li^yi!iS!S!^iji!S000o000000000000000000000 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOZZZzzzzzzzzzzzzzzzzzzzzzz 

H««n»®r*N«(,*«ff‘9,oo*OHa'Oin9,a'(7'ncDcoinof'«Nnff'nH»in»irtoff'9,a>9i®'oa'e(n(N#<r«oioo(»i 

nfi»<,fKnnin'Oin<,r«^<*><,M<,^nf,mnin^,noin«'nn<inN(»i*n^in/in^nfnn<«,r«(*3n(nnrt<"i,n»*i«r 

rH 

<ro*o*in®m*»®inHrm®®®r'>r--09o,3rH09e9©®®or-in^e'>r^o9®oa'ina»»>e'3'Oomr>09P'09o«®otns©®aM7sr»f^® 

HHnMrtHHH«nnHMHHHHHHrtH(NrtrtnHHr*WHnH  HHHINNCKNMNNNOHNrtHHHHtHHNM 

«oo9in9,03HinHO<J'Hfor»rt^'Oino'H(T<nfN03r'r'9,'OinHOO'Cn®(N«oo'a!)r-H<<ininonin'O^H9'a3HM 

nflHHd  MMnHHmHtlHHHHHHHHHNHHNHHNM  *-H,r,_,CNrHOI.-<OH.-H09r-«  (NHMHHNNHrtOI 

r-»no^«nr»mo^fn^,f,«ci^,^r*®^ioina*a\^-»o^iOf*>®r^»-«f^oo®<M®®^'«®rs<M®or>i®C9<^o<n«®®inflo<N 

HHHH  H  ♦  fl  H  H  rt  M  rlO  H  fi  Hf,  H  N  M  N  rt  H  HN(N  09  f>*  H  H  H  rt  H  <N 

*»M*r»r'*mrH*'0»®©090s®mr’3*,mo>®r-©®'©o®r**®«*>inr-©0«.©09®09®09ir»®w®<n**<*<7*r*®rHm®®®09<7>® 

HdHHI'IH  fl  H  «  Ct  W  «  09  rH  09  -H  rH  HHNHHHfVHHrt 

o«^r«r*®»r)Ovr9®fnno®^T9®^®r>*in®av'r«y'mmr'C9^r-tpn®OQor«<y«\o*oin<73^®r-®r~-«‘r«<yvo®r~m®mr*«N 
rH  .-h  09  rH  rH  rH  r*  c9  rH  09  c*»  09  rH  r*  rH  rH  «h  rH  nnn  ,h  rH  rH  oh  rH  rH  rH  HMHMHfiHnMHHHHH  09  ,h  rH  rH  .h  rH  rH 

H 

«irt»»ninrtn^<MftMoift'»(nnna>c»®nrtonnr-ci^inwH»or*HinoHH»iN^Boo»H®^in*in«in^n(*» 
rH  09  rH  «h  h  -4  «-»  r~  rH  -9  *-i  09  «-h  09  *-9  rH  rH  «-t  09  09  09  09  HMHHnn  HNHNiiNHXHHHHnn  -h  -»  rH  rH  -*  -h  rH 

^»»n®®<F»r,-<*f9m®*,O9rHO9ioo9r**®©*-9o-©in«ft<o®o*(*>*,t/30*o9©®®mr^®r'*O9<Mn®o>mo3»noinmvorHr*®a*<7» 
HCln  N  HNHMHNHHHHNHH  P(  m  Pt  h  h  hh  hn  hhhhhwnhh  hhhnhhh  hh  hh 

Hfir*-r®'OP*c*n'0®^ooa)0(N9,Ort'oa3r**r*rt9,M®nH«ir(*ip-Hnff'r«or**^r'riinf«nooi,'(*3n9,»-*HnH 
HHHH  H  H(1H(lHH  HrtHMHfl  HHH  HHHHHH  HPtHPIdCIHHHHHHHH  HHHHHHHN 

eHn®©*,mc9®p»<n*-9ir»®^\oi«©®®©o’>«no9®®©<no9*>o*®m^r--o909®^p-ri*ir»<'no9ino9®©»na>c909rHHr<*<nrH 

rH  +  net  HHI1  rH  rH  HmnHflHfNriHHHCldH  HHHHHnnnH  HHPtHn  09  -H  HHCtriHHCI 

>o®ffc*»rt*HO®p,»^c,9®oo9W«Hr,»o^oo,)«»)»H®®»H«H»no9®4r>o~®r9®m»no9o*a'<nooo90»o*o9<js<y»a'®^o9«o*or~ 

H  H  H  H  H  HHHH(iHf9f1HnHHf9HP1HHHHHH  rH09C9rH04rH09rH03rHrH09  03rH  HH  HHHHH  H  (9  H 

<r<nr«-®r-<*^,n®rt^Mnr*e»mor*^in^®®^#n®c9i*>fnor^»»^^Hr«»®®»nom®irKnminHfofn®r»OH»n*H«H 

HHH  H  HHHHHHfl  HH  (*3  H  H  H  H  <9  PI  HHH  (*)  H  H  H  HHdHHHPI  H  H  H  HHflHHHH 

HH^Hp«^inp‘ieH9,«OP*«®nff'«nnHH«'0(9H(7'asHHinn<T'r«rtnff*nnHino'nino(Mnr‘npiH®iD^9> 
(1HHH  HP»Hf9f9H  NHHHWHHHIOMMHHPIH  PlHPIHHHHHHPIHHH  09  09  rH  NP9HMNH  rH 

09®»no»^®®®®*‘®rH^09mrH^09M’)inrH.-»rHO®rHrH®®o9®®r-rHe9«'»n®©®ei®®«<TvrH<n©<r»®«09rH®®r» 

H  HH  H  H  MHnnnn  rH  09  09  09  09  09  09  H  H  H  H  H  H  H  H  H  f9  H  HHHH  H  09  H  H 


HHnP9H«»H 


^r»m®09  09io»9®r-i^^^®n^,»H09»n®^09^'ir*09r9fn*-9r>r'^*»n^i»Hm®m^^»nr*»09® 
09  rH  rH  09  I*109HH  09  rH  rH  n  rH  rH  ^  09  09  #H  HHHH  rH  HCIH  H  (N 

mo«9'9ir,>Hor9Hoa«noooooDiniAinHOOioo90ni3iMTt®r'OOH®HeB®®«o®oin9ioiine9<oio«H 
rH  rH  09  ^  »H  »H  ®  09  H  H  (*■)  09  H  rH  09  09  m  rH  09  rH  rH  >H  09  HHHH09I9HH  HHH  rH  09 

H 

HHO»no900r>asa'iO'0090«o-P909r‘r‘HH«0'00r-H'»ono9fl'f'-onOHonoi'0'OcMJ'0'®Hp'09,090'Offvoi 
rH  rH  09  HflH  ®  rH  H09  H  l»IH  09  HHHH  HHO  rH  09  rH  09  rH  rH  rH  rH  09  rH  rH  rH  rH  09  09  rH 

rH 

©o9*|^©^o9®*>H©<n®HH^HO»o*irnnom®®o*-®«o».©ine»or«-*.O'0'rHO*,«-9O®H©*.u->mo9®®©mr*om 
rH  rH  09  09  rH  rH  09  09  rH  rH  rH  rH  rH  rH  rH  rH  09  rH  rH  HHHOIHHH  09  09  rH  rH  09  rH  09  rH  rH  09  rH  rH  rH 

rH  09 

C9©ininiftof9p«.o9P99>«9'r,'HioioooMnin*H09tn«®o9P*HHP*inoo^ff>«HOrtinp*inn9'0®®9>»r«-intf»oiin 
rH  rH  09  09  H  (*)  >H  9*  09  H  H  H  H  H  Ol  rH  rH  09  rH  rH  rHrH0909rHrHO9O9  09  09  rH  rH  rH  rH  09  rH  rH  rH  H  rt  rH  H  rH  09 

rH  rH 

«r09ino^H09Oini»>rH«®r9^tDn®'fiP9r»(noH,fl0OinHi^OH'OBnnO'Hp9r*n'0ininH,r9HH09p'©B09HO'O 
rH  P909HHHHH009HH0909HH09HHM  09  09  09  rH  rH  rH  09  09  09  rHrHrH0909rHrHrHrH0909rH09rH  HHHHHP9M 

rH 

®^®®09  09P>»^©rH09^e9fne*©rH<n®rHPno»m<ho9m»H09®®®©f^©©o9m©©®®«mm^<Mr»e\®o9®®rH®o90 
rHrHrHrHrHrHrHrH09rHrHrH09rH0909rHrHfr)rH  09rHrHrH090909rHrH»T|rH  (*3rHf9Hf9HtHHH09('9HHH  09  rH  rH  03  rH  09  09 

ffiu>on«)0'OcD'<MOfriMPi«oO'Oiop'<ynooi/,>9'no909HO»o-nor‘'Oi*ion«OHP‘tf'p'rHHoo'<JO'00nn<no 
rH  rH  rH  rH  rH  n09rtrtrlHH09HHHH  09rHrHrH09rH0909090909  09  09  rH  Ol  rH  09  rH  09  rH  rH  09  09  09  09  rH  09  rH  03  rH  09  09  rH 

<hw>^rH^pn^rH09®®©®®®<no9rH©o><nin«rrHmrH®«nifta«p»^©©o*©p-o9mtf»o-a»rH09mr>o9®rHO-®o)®^no9 
HO9Hf9r(r9HHHHrHHHHOlHf9HHHMHHrHn090l  HOIfl  09  09  P'1  H  09  rH  rH  09  09  09  rH  «r>  rH  rH  09  rH  09  09  rH  09  rH 

Mrt'fi®«r9H09®9>«0»OM^rH9'O«®®(n*n^O91inn9'«nrtOH«M*)n0»OO«P*H,(J'H(JSHa5l0OM3O(*1O 
rHO9rHO9O9rHrHP“rHrHO9O909rHrHO9rHrH  rHl*309HH09rHHH09nn  09rH09H*rH09rHrHrH09rH  09HH09rHHHr4n 

©®-»09rH®ir>rH0»®rHprrHmp^Hr«H^m09(*>Hr©®®09r»®<9,r»^®f*‘©^©®®in©09©®rHoe»^®or*-®®Hr®®in 
rH  H  WH  HnH^P9H09<9O9HHHHHHHrtnHHf9p9P9HHP9H  rH  09  rH  Ol  rH  P'3  09  09  rH  09  09  Hr9P9HHNHrHHH 

09 

in»r^®ooHf9n9,in<©p»«oorioiM9,in«f*»0'OHOi(»)^in*MDOOHPiM*rirtvOP‘a3a>OHno)9,in'or**®ffio 
f'.p^O*P'r«-®®®®aOaO®®®®ON(T>Cr©<7'OS<T>^>^©OOOOOOOOOOrHrHrHrHrHrHrHrHrHrHOt0909090HfH0909090«On 
09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09  09<n0>0\pn030'tfnnnn<,^r,>mo^n<nr»0,,»nm0»r^fnfnr^pn0»030,imr^ 


13 


Table  23.2  (Page  1  of  4) 


NAO  . DPX  Hourly  distribution  of  detections 


NOMA*  Sei  ft*.  242*3 


M*y  1993 


•  p 
•u  u 

•  tj 


>■ 

«  >. 

>^o  « 


>» 

«  >. 

>1-0  * 


>> 

4  >. 

>vo  m 


>.T3 
>.>,4  n 
•  ■Oil 


*8  4*1!  4  4  *3  3?  4  2  4  4*3  38  4*2  f  i3  38  4*2  S' 4  >3  «ll*6<«,D#«4l.4  4TJ®*«&<l«T3tli*&<l«’Oi 

3sB?13|M3??35S33??!^SS3?1S|323?13|S33?|!|S33??S|S23?^S5 

(IflNMnOOOOOOOOOHHHHrtHHHHHMdMdddftdNrtnnOOOOOOOOOHrtrHHrtHHHHHIN 

I  22222222222222222222222222222222222244444444444444444444 

m  Mn*«««(rt^«rt®H«)ino«‘«nci^n<,iO'OHirt^^ioo(n»i09»oin<,«nHO®»r.Mn«irtoo«tfMO»^oo 

a  ♦nmOHr»(,nr«wnoa'>er«onfiirt^o»ffMn(ninr«on«'0(noff<^inO'0»o*tnr'OMr»iniftrif^^H^HMOfNCbn 

m  4fl>rtHNio»N<n(»iooin*OH^ino'O^Hff>nnnri»®inHie»oi«iff''0«'Oi,i4inom»p*'OH«irnrp|)H<Moo 
<n  c* <n  <n  cn  <n  ^  .h  ^  <n  <n  rtdHNMNM  «-*  ^-t^fNfn  nnmhm  hhhhhnh  r«  -« -« .-i  -«  c* 


o  mo®p*»OH«*ff.oowMn«'»OHH'0«<roiftr-r«®^^  in  m  tMnmoo^r-t^cDcoojMnor^voos 

<N  W*  fN  rinrlHH  »-4  *H  HHHH  nNNHnflHflH  HNNHNH  rt«  ro 

<y>  octr»or*fn«no«B^i40fl»4*c«r*.r^4*r^r^rt<n4*^a*®iooooo^o»cosotf'0<MP^inooooosfnor^r.Hr^o*«^f».f>«srO'« 

s-l  HftnrfHHHNHHrt  rlHHrtHHHHHHHHH  MHriHnnHHflHrfHHHHHMHrtH  HMHrHH  fN  «-»  fN 

to  04-i'.c«o,*©*-«toiotN<N©ininr‘4*r'»^fnr'-r«‘r».-i<*>fofn©®in©©fnfNi©^^®«s-i®inr*r*^a»®in<N©<N*Him©t©©® 

-4  HC»HC»Hf«r»HHC»PI  HHHflNHMHHHHHHHH  H  H  f»  (1  HH  (*l  H  H  H  HMH  Hrt  <N  fN  fN  *H  ~4  »~4  »H  •— 4 

r»  n^r.ioor*oc«oos-*4'min»ft«n»o«n4'ioomoososcio»r-oosoasr«o4v«4,4,»Or-4inoins-tr*^os®4'^ir«in®fno»o 
.H  fN  fN  <N  »-l  fN  CN  <N  <N  <N  <N  HdflPIHHNM  HHH«  <1  H  CIHrt  rtf»  H  rt  P»m  ,-t  fN  .H  —4  fN  <N 


fN  fN  m  <N  fN  fN  HflrtHHHnHMdHMM^rtNrtfnHHHCt  — 1  fN  —4  -4  fN  <N  Hl*| 

oiino«fi*»inn<n«oo«<NP»40HOnnN,a>4nH«o*«40»f,(oiHn4innHr*HHOO'is*H»B4o<ninff>^co 
HnmrtriHnnfiHHfn  p*(NH<N^<NH»HHHHHH«(*>fNH(i(,iHrtnnf,«H(*iHHHfN  •”>  •-*  .-1  fN  ni*i 


O 

CM 

V 

00 

Cu 


HP»nN,n4p.4»OHf»(»»4in«p.®oiOHf>in^(niOpv»ffiOHf«<nN,in'tfH«(*i4in'OP‘®0'OHf«n4in'Ops®0'0 
(*)nnrmnf<lH*m44<r4'N,4,4r«ri«i'ri/MninininioimniniA<6io4iO4ioi0  HHHHHHHH-irtN 


N 

n 

N 

*! 


NAO  .DPX  Hourly  distribution  of  detections 


NORSAftSciR*.  2-92*3 


May  1993 


4) 

Ai 

« 

Q 


C 

3 


W 


ggSgggSSSSS||2|||||||||||||||||||||||||S|m|SSSfeSSSSSfefe 


®04r'><na*®*''0n®04©c*©o®o'>o-®f»04®®®®f^m<'>©<©.-904r'-©infn®in.-4inr*-®\©04©®op'>©o>©c"©©'-'© 

OH-»p*HHMh-r«a)r«  ♦  n  a»  »  04  ®  ®  a*  ®  in  vo  ®  ♦  a*  ®  ®  o  o-  04  ®  o  r-  c*  r-  tr  ®  m  so  o  \o  +•  oi  on  m  ®  o'  04  o>  04  «-i  o> 


tn  m*,^OMno4®<r*o**,©oo4^oinf">o>r**,®r-4©®soo*mo»040»®o>ina\®in»H©r-*,®o***e9''’»p-.-HOO®^«n®otn® 

ON  *H  rH  rN  rH  W  •— 9  «— 9  0404  •-«  (NHHtNHH  fN  tN  (N  fNH  P»HNH(MHNHHH  Hrt^ON  .-9  04  04  OI  ON  -H 

04  o4e-4Wwno4®09f,-®^,©©^*-904®*‘®.-9ino>®r'-o4©in®r-®oi04®in*>4n<o.-9©r4mo">oin9rm®©©©9nina»r'-o-H® 

04  HnnHriMHHCin  HHClnfNNHHHMHMHPtHHrtnrJNH  HH-4MN  -rtNH  0104 

.h  <r»r-04r>.**«in049'"ivo©oo4«.-9in*'«ininf-OMnf-c**'CTN.-9in©®9''i»ncnr-o%©®4no'®oin^oJoi.HOo®\r-»n«'Oin** 

cs  HHH«HHH(Nnn  in  04  -*  a  -t  .h  r-9  -h  ,h  oi  os  rs  os  *h  -i  os  os  inhnmhnhh  nnomn  hnwcn  rs  .h 

o  .H®o4rH©®®omr-©©r>®'4,®04r-®a"Oi^04o-N*fn-i»<Mno4r-o*,0'>-*o-.-9r~»no'©*H\c®>n*-oo©r'-‘D9*io®9''i 

04  rH  04  04  rS  04  -I  O")  04  c-l  H  (N  H  H  H(NNNN(nHHC«MHH(NHP»HNHNH  IN  rs  04  09  .H  04  •-*  .N  «H  •— 1  r-< 

as  o®^>nin^HHoooHinooof4'in®o®iooc"or9inic(Nn^^o«^oin®inooo*«o'HOOfi®in{Noort 
.H  09  MHHHH  04.H.H09— 4r-90904--9.-«04*H  M  (N  f9  N  *H  »-l  H  H  (\  #H  H  «  <H  rH  09  rH  HCJ 

®  rt4-«4-ooHO'OinHOtf>«*HOWHrtm^«M«Mn^O'<ono®inHffMoo04'ino«M*iifl»oon»HO'09'n 
T— I  HHrlHHHH  m  rH  rH  HHNHH  M«N(*tPlrtHflHHHMHHnHNnr»HH  r-9  rH  rH  04  09  HHMH  04-9 

o»  <^<’90fia34'«4,OHOOof>»oooDioOHHOomn«iJ'»ff'noifiMff'ff>Hoi^Mino*o^«ff'OOinoNinoino 
r-9  Ol  *H  *-4  *H  HHHr9r9«<HC4M«HH«r9HH«r«<HHNMfSHH  HOlHHH  .-9  rH  04  04  r-9  CN 

®  9r>«rfomromo'<7'®»HO>o-H<y>*Hrofna>®^,{y'®®o<''tf>r-'oo-oor4^inrno4oo4^ooo'04m®04^-toor'r4r-vc>o— 

rH  HHC9HHHH  HCI  04  HHHHHH  HOC9(NMH09HH0909H0«nH09HHHHH(Nf-(H09r9  H  rs  IN  09  HH 

in  ioin®nP‘0090HOH04T'<}'«>OMno9noMe(}MNosrt4ioin(Mo^(Nnff'Ho-(Oina50(Ni*ir9<7w»ioo'OH4'<NO'fi» 

*-H  04  04  r-9  04  ®  04  04  r-9  —4  rH  rH  04  rH  «-9  rH  04  rH  04  04  rH  04  04  r>  rH  rH  r-9  04  rH  rH  »-9  rH  rH  «-9  rH  04  pH  04  rH  rH  *H^09(N  rH  rH 

■*•  f44'09®HnP'fninoooo9)(OH(,9CD4'Oino90or'nn®*omo9«inoHa3H(noioooH(j'OOOOH<j\09»on'0 

»-9  04  04  rH  rH  rH  04  rH  04  .-4  HH(NHr9r9HH(N(NO«HHHHHHHHnr9  04  04  m  09  rH  IN  04  09  04  09  rH  rH 

m  m»Hinfi09Hin«Hoooo4,NoopirH04'Hinnoff'ffNManooifiinooooo90oo40'CinoooHn4,0(»83 
r-9  04  O4O4rHO4O4rHrH0S  HH04HHH  H09H04H04H  HH  r-1  r-9  rH  m  04  04  **  04  fH.-HrHrH.-H  H  H  09  09  »H  rH 

04  m®®c*r*‘^or-9»s^40000404H040*®®e*®<nin®®p**r-^>o»n«04®o>®^'0®^oo»Hmp»«H®oo®04r*-^o«r>-> 

rH  f-4  04  rH  04  .HrH  rH  rH  rH  04  rH  rH  04  04  <•»  04  09  *■  rH  rH  04  rH  04  rH  rH  ^  *H  04  04  04  rH  Ol  rH  04  09  H  H  OS  04  rH  rH  04  rH  rH 

<h  inss«(n'fi«04o«ooonna'Hos\oo9(Difia'snsO(MnHMrMnoeHO'soir)0«inMnsoooHSioon'C 
H  H  04  04  hHS>HH  H  HOIHHnHHH(«lH04Hfl  «*>  HHHHHHHHH  04  04  04  rH  rH  rH  rH  04  04  H 

o  o4mHP»®.-9*,o4O4m©o©.-4®©O4OMnc**H®0'O4®in©cs*'9r<no-orN“r>-min4n*,o49n©oo9®<no>o©»*o4o»«©®® 

*— 9  1-9  on  rH  «n<-4«H04  04  HHHHHH  rH  H  09  n  H  09  04  H  H  09  H  H  HHH0I09  HOIHHH  rH  09  04  rH  rH 

tr>  om<no>r>090<*i(nfno40o«H®in®o><-9®o««>04®o404«>oo4®9s®09<NOkn®®inr»or>®<H«Hmoo®o>o>'®o9i® 

O  rH  n  04  H  rH  04  04  rH  rH  rH  04  rH  rH  rH  Ol  rH  04  rH  rH  «H  «n  «-»  «*4  04  04  04  ^-1  *H  04  HHH(SH  »H  OS  *H 

®  «mo4  04  04o*^H®®oo©^o'0'nnno404HCs®©o*H9r)«04fo^mm^®fn©'tfr*OHr9rsff»r^p^©oo4Hinwo®® 
O  H04HHH  H  H  04  HfnHHHH04HS104HH  •H0404rH«H*-9.H^H04-'«  H  09  H  04  H  09  H  04  <-9  <-9  04  -9 

o404*-9f»»®r-in®a*©®©©<*>r--a9^o’>in®^o<,ntno4®<no4rne'r,>a»04®f-**sp-«4T'><'4mo4©©®«-9m©©^r>r*‘^©®r> 
O  *H  riHHrtH  0S«H  HH09H  HrtHf  HH04H  04  09  --I  -4  04  »H  04  04  H  09  H  IN  H  HHC9C9 

®  o4»nmmo)ff*^-90«*«-ip«-©©©®04mp»p«-^9r-o,>®®oi®®©in®®04®o^9rjo*in®®^,®®09in^«»H©©a'04«n©©p^ 

o  C4«-»^S^S04  H  H  04  HHH04H  04  04  ^-9  ^-9  «H  09  r-9  .-4  H  H  H  H  H  04  HHHHIN  H  H  H  (N  »H 

in  ns'©«Oieoooi/,MOoof'a)'os‘irr^«oa'«oin«0'0»ininH43Hoo«ino'0'®HOnHOO'(No©ono«r*in 

O  *H  (NHHHHHHfNfO  «H  *H  04  *H  m  H>  r-9  04  H  H  04  H  04  04  H  04  04  H  04  H  04  H  IN  09  ■— 9  r— »  04  HrtfnOI  »H  OJ 

•s-  ®or»®®o*r-®m®ino©^,m®po^in®©09®^9®H^©o-®^^909®©p»>^,®4n^9®<,©o©»no4^9©p«-»n*rr-^9®-^ 

©  H09H09HHHHH04  04  H  04  09  H09HH04n09H04H040HrtH09  04  04  04  09HHN  04  04  O)  m  04  ^  04  04  ON  HJ  —I  04 

m  vOHor'OsHOinr909004,cDnnioo9'cosno9ooinoifl09(7<oo«inH09iNooinnHOnn'OiNioossHO®oin 
O  H09  09HH04  04  H09H  04  04  »H  »H  04  <-9  <H  <H  04  04  r>  04  04  «H  04  <H  «H  HHOI04  04HC4H  04  04  04  04  .-9  0>  04  HHOOINININ 

04  p*«-9in®®m*H®*H©©©©04o*04®040*04in®04a»04»n0404®inpoo4©»H<T*04  04r,>o,)m©©©®®9rff'<9r©^'»n*H^o>C'in 

O  •H040fltH*H«HrH040404  ^9^-9fS^9^904r-9^-9  09r-*0>0404nr>0l09  04  04  «-9  04  04  04  w4  *4  04  <-9  r-9  04  Ol  04  HH4>O9  09hCN 

«h  ooHr-(NHns,0909090oifl(NOH©sn^n(N<Doos,ofl3oocDOf-ofl3onB'flS'soovoa"9)ins©inH(Dr'Of9s 
©  *-9  04  04  04  04  04  *H  rH  04  04  HHHI*|0!04HHHHH09  09  04HHH04  O4HHH04HHmHlN  H  IN  H  IN  IN  r-9  04  09  04  H  «-9  04 

©  op‘innOHOf*‘aa«»ooinoooo90H09inso9Hsoi(j'Mnoinsa3<r!HfNns,n>ooon(MDnf,'00'OOi»«oin 
O  rH  04  04  OS  04  04  *H  r-9  OS  r-9  (NH04r9f9Hr9HHHni*1HrmHfNn  m  H  04  04  m  IN  04  OS  04  rH  09H09  04H  nOOHOIN 


>,  r-404oai«tn®o-®<T90rH04(n^in®r-®9>o«H04(n«in®o-®^Or-9049n<s>m®o-'®^Or-f(Nm9rm®r^©^OrH04iin^in® 

K  iNr9fNfNiN04090909nnnoooonoosss*ssss,ssininio^inimnininin'0N>i0i0«)^'0'0'0'fir'NOM^r-r- 
Q 


15 


Table  23.2.  (Page  3  of  4) 


. dpx  Hourly  distribution  of  detections 


NOftSAft  Sri.  ft*  2-92*3 


May  1993 


« 

u 


s 


a 

i 


«  «  >Vt>  *  « 

‘SS'EfrMSSS'li 

- - — 

illli 


nsmsmssn 


-ir3O4««*>O'O*r9©«rr**>0»-t 

r>o*»or-^,«>inf^«30r-*oo»H 


•  HrtH  O 


S>»  p*>  oo  &>  ©  «h  ra  m  ®  r-  oo  o*  o  O 

«  r- r*»  oo  »  ao  oo  «  ©  oo  ®  ®  oo  O'  < 

z  O  z 


n  r* 
**  r* 
v  tn 


o  O' 

c*  «h 

o  03 
(N 


rl 

O' 


o* 

f*^ 

OJ 

04 

04 

m 

o* 


o» 

o» 

CM 

O* 

03 


Oi 

o* 


Sfi 

34s 

«  u»  o 

£-a-s 

“M  8 

!§f 

11 1 

21* 

111 

CC  z>  s> 

o|  & 
z-s  « 

'S  o'5 

w  VM 


c 

o 

•3 

3 

,£* 

•fi 


00 

si 


a -8 

ftsS 

^  ri\ 


3^ 
■°  o 

s* 

•3 


SP 


8 


■§1 

"I  e 

*■3  - 

tsa 


_  a. 

3* 

<*-  S* 

O  13 

s-g 

3  % 

C/5 


£2 

■°  S 

a-s 


‘S 

Q 

3 

N  <*, 

st  2 

■o  3  3 

g  jl 


G. 

1 

I 

£ 

I 


16 


3  Operation  of  regional  arrays 

3.1  Recording  of  NORESS  data  at  NDPC,  Kjeller 

Table  3.1.1  lists  the  main  outage  times  and  reasons. 

The  average  recording  time  was  99.69%  as  compared  to  98.55%  during  the  previous 
reporting  period. 


# 


*) 


Date 

Time 

Cause 

11  Nov 

1925  -  1941 

Transmission  line  failure 

31  Dec 

2310  -  0000 

Software  work 

08  Mar 

1200  -  1216 

Hub  maintenance 

15  Mar 

0934  -  0952 

Transmission  line  failure 

Table  3.1.1.  Interruptions  in  recording  of  NORESS  data  at  NDPC,  1  October  1992  -  31 
March  1993. 


Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


October 

99.97 

November 

99.96 

December 

99.89 

January 

100.00 

February 

100.00 

March 

99.91 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J.  Torstveit 
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Fig.  3.1.1.  NORESS  data  recording  uptime  for  October  (top),  November  (middle)  and 
December  (bottom)  1992. 
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3 2  Recording  of  ARCESS  data  at  NDPC,  Kjeller 

Table  3.2.1  lists  the  main  outage  times  and  reasons. 

The  average  recording  time  was  99.61%  as  compared  to  99.25%  for  the  previous  reporting 


period. 

Date 

Time 

Cause 

27  Nov 

2123  -  2214 

Hardware  failure  NDPC 

28  Nov 

2153  -  2230 

Hardware  failure  NDPC 

16  Dec 

1057  -  1147 

Hardware  failure  NDPC 

18  Dec 

0641  -0912 

Software  work 

31  Dec 

2311  -  0000 

Software  work 

01  Feb 

1821  -  1927 

Satellite  link  failure 

01  Feb 

2054  - 

Satellite  link  failure 

02  Feb 

-  0516 

Table  3.2.1.  The  main  interruptions  in  recording  of  ARCESS  data  at  NDPC,  1  October 
1992  -  31  March  1993. 

Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 

October  99.99% 

November  :  99.78% 

December  :  99.41% 

January  :  100.00% 

February  :  98.49% 

March  :  99.98% 


Fig.  3.2.1.  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 


J.  Torstveit 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  October  (top),  November  (middle)  and 
December  (bottom)  1992. 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  January  (top),  February  (middle)  and 
March  (bottom)  1993. 
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3 3  Recording  of  FIN£SA  data  at  NDPC,  Kjeller 

The  average  recording  time  was  98.78%  as  compared  to  93.75%  for  the  previous  period. 
As  can  be  seen  from  Table  3.3.1  below,  the  main  reason  for  the  downtime  is  transmission 
line  failure. 


Date 

Time 

Cause 

23  Oct 

0434 

-0935 

Transmission  line  failure 

10  Dec 

0547 

-  1314 

Transmission  line  failure 

18  Dec 

1337 

-  1422 

Transmission  line  failure 

18  Dec 

1427 

-  1447 

Transmission  line  failure 

21  Dec 

1021 

-  1035 

Transmission  line  failure 

26  Jan 

1636 

-  1651 

Transmission  line  failure 

26  Jan 

2104 

-  2137 

Transmission  line  failure 

27  Jan 

1255 

- 

Software  work  in  Helsinki 

28  Jan 

-  1440 

29  Jan 

2316 

- 

Hardware  failure  NDPC 

30  Jan 

-  0030 

31  Jan 

2213 

- 

Transmission  line  failure 

01  Feb 

-  0712 

Transmission  line  failure 

10  Mar 

1043 

-  1057 

Transmission  line  failure 

10  Mar 

1334 

-  1349 

Transmission  line  failure 

10  Mar 

1525 

-  1544 

Transmission  line  failure 

Table  33.1.  The  main  interruptions  in  recording  of  FINESA  data  at  NDPC,  1  October 
1992  -  31  March  1993. 

Monthly  uptimes  for  the  FINESA  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmission  lines,  data  center  operation)  affecting  this  task  were 
as  follows: 


October 

99.13% 

November 

99.98% 

December 

98.79% 

January 

95.99% 

February 

98.91% 

March 

99.87% 

Fig.  3.3.1  shows  the  uptime  for  the  data  recor  ding  task,  or  equivalently,  the  availability  of 
FINESA  data  in  our  tape  archive,  on  a  day-by-o  t-y  basis,  for  the  reporting  period. 

J.  Torstveit 
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Fig.  33.1.  FINESA  data  recording  uptime  for  October  (top),  November  (middle)  and 
December  (bottom)  1992. 
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Fig.  33.1.  FINESA  data  recording  uptime  for  January  (top),  February  (middle)  and 
March  (bottom)  1993. 
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34  Event  detection  operation 

This  section  reports  results  from  one-array  automatic  processing  using  signal  processing 
recipes  and  “ronapp”  recipes  for  the  ep  program  (NORSAR  Sci.  Rep.  No  2-8839). 

Three  systems  are  in  parallel  operation  to  associate  detected  phases  and  locate  events: 

1.  The  ep  program  with  “ronapp”  recipes  is  operated  independently  on  each  array  to 
obtain  simple  one-array  automatic  solutions. 

2.  The  Generalized  Beamforming  method  (GBF)  (see  F.  Ringdal  and  T.  Kvasma 
(1989),  A  mulitchannel  processing  approach  to  real  time  network  detection,  phase 
association  and  threshold  monitoring,  BSSA  Vol  79,  no  6, 1927-1940)  processes  the 
four  arrays  jointly  and  presents  locations  of  regional  events. 

3.  The  IMS  system  is  operated  on  the  same  set  of  arrivals  as  ep  and  GBF  and  reports 
also  teleseismic  events  in  addition  to  regional  ones. 

IMS  results  are  reported  in  section  3.S  and  GBF  results  in  section  3.6. 

In  addition  to  these  three  event  association  processes,  we  are  running  test  versions  of  the 
so-called  Threshold  Monitoring  (TM)  process.  This  is  a  process  that  monitors  the  seismic 
amplitude  level  at  the  four  regional  arrays  continuously  in  time  to  estimate  the  upper  mag¬ 
nitude  limit  of  an  event  that  might  go  undetected  by  the  network.  The  current  TM  process 
is  beamed  to  several  sites  of  interest,  including  the  Novaya  Zemlya  test  site.  Simple  dis¬ 
plays  of  so-called  threshold  curves  reveal  instants  of  particular  interest;  i.e.,  instants  when 
events  above  a  certain  magnitude  threshold  may  have  occurred  in  the  target  region. 
Results  from  the  three  processes  described  above  are  used  to  help  resolve  what  actually 
happened  during  these  instances. 

NORESS  detections 

The  number  of  detections  (phases)  reported  from  day  275, 1992,  through  day  090, 1993, 
was  36,783  giving  an  average  of  202  detections  per  processed  day  (182  days  processed). 

Table  3.4.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS. 

Events  automatically  located  by  NORESS 

During  days  275, 1992,  through  090, 1993, 2074  local  and  regional  events  were  located  by 
NORESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  11.4  events  per  processed  day  (183  days  processed).  63%  of  these  events  are  within 
300  km,  and  88%  of  these  events  are  within  1000  km. 

ARCESS  detections 

The  number  of  detections  (phases)  reported  during  day  275, 1992,  through  day  090, 1993, 
was  91 ,474,  giving  an  average  of  503  detections  per  processed  day  ( 1 82  days  processed). 


♦) 


•  • 
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Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  275, 1992,  through  090, 1993, 3359  local  and  regional  events  were  located  by 
ARCESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
18 .5  events  per  processed  day  (182  days  processed).  43%  of  these  events  are  within 
300  km,  and  84%  of  these  events  are  within  1000  km. 

FINESA  detections 

The  number  of  detections  (phases)  reported  during  day  275, 1992,  through  day  090, 1993, 
was  58,399,  giving  an  average  of  321  detections  per  processed  day  (182  days  processed). 

Table  3.4.3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 

Events  automatically  located  by  FINESA 

During  days  275, 1992,  through  090, 1993, 2835  local  and  regional  events  were  located  by 
FINESA,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  15.6  events  per  processed  day  (182  days  processed).  65%  of  these  events  are  within 
300  km,  and  87%  of  these  events  are  within  1000  km. 

GERESS  detections 

The  number  of  detections  (phases)  reported  from  day  275, 1992,  through  day  090, 1993, 
was  7,456,  giving  an  average  of  182  detections  per  processed  day  (41  days  processed). 

Table  3.4.4  shows  daily  and  hourly  distribution  of  detections  for  GERESS. 

Events  automatically  located  by  GERESS 

During  days  275, 1992,  through  090, 1993, 717  local  and  regional  events  were  located  by 
GERESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  17.5  events  per  processed  day  (41  days  processed).  72%  of  these  events  are  within 
300  km,  and  90%  of  these  events  are  within  1000  km. 

Apatity  detections 

The  number  of  detections  (phases)  reported  from  day  275, 1992,  through  day  090, 1993, 
was  90,712,  giving  an  average  of  501  detections  per  processed  day  (181  days  processed). 

As  described  in  earlier  reports,  the  data  from  Apatity  is  transferred  by  one-way  (simplex) 
radio  links  to  Apatity  city.  The  transmission  suffers  from  radio  disturbances  that  result  in  a 


large  number  of  small  data  gaps  and  spikes  in  the  data.  Although  the  communication  pro¬ 
tocol  may  correct  such  errors  by  requesting  retransmission  of  data,  this  cannot  be  done  at 
Apatity.  For  such  error  corrections,  a  two-way  radio  link  is  needed  (duplex  radio).  Some 
of  the  detections  are  consequently  due  to  bad  data  spikes.  However,  it  should  be  noted  that 
noise  from  cultural  activities  and  from  the  nearby  lakes  cause  most  of  the  unwanted  detec¬ 
tions.  These  unwanted  detections  are  “filtered”  in  the  signal  procesing,  as  they  give  seis¬ 
mic  velocities  that  are  outside  accepted  limits  for  regional  and  teleseismic  phase 
velocities. 

Table  3.4.5  shows  daily  and  hourly  distribution  of  detections  for  Apatity. 

Events  automatically  located  by  Apatity 

During  days  275, 1992,  through  090, 1993, 2691  local  and  regional  events  were  located  by 
Apatity,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of 
15.0  events  per  processed  day  (180  days  processed).  31%  of  these  events  are  within 
300  km,  and  66%  of  these  events  are  within  1000  km. 

Spitsbergen  detections 

The  number  of  detections  (phases)  reported  from  day  275, 1992,  through  day  090, 1993, 
was  19,359,  giving  an  average  of  150  detections  per  processed  day  (129  days  processed). 

Table  3.5.6  shows  daily  and  hourly  distribution  of  detections  for  Spitsbergen. 

Events  automatically  located  by  Spitsbergen 

The  installation  at  Spitsbergen  is  not  completed,  and  automatic  processing  for  location  of 
events  has  only  been  done  for  experimental  purposes.  See  also  section  7.5. 
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Figs.  7. 1.2-7. 1.4  show  the  event  distribution  in  Fennoscandia  for  all  magnitudes,  >  1 
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3.5  IMS  operation 

« 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  December  1989 
and  was  operated  at  NORSAR  from  1  January  1990  for  automatic  processing  of  data  from 
ARCESS  and  NORESS.  A  second  version  of  IMS  that  accepts  data  from  an  arbitrary 
number  of  arrays  and  single  3-component  stations  was  installed  at  NORSAR  in  October 
<  1991,  and  regular  operation  of  the  system  comprising  analysis  of  data  from  the  4  arrays 

ARCESS,  NORESS,  FINESA  and  GERESS  started  on  15  October  1991.  As  opposed  to 
the  first  version  of  IMS,  the  one  in  current  operation  also  locates  events  at  teleseismic  dis¬ 
tance. 

,  During  this  reporting  period,  two  new  data  sources  were  made  available  to  IMS:  the  array 

near  Apatity  on  the  Kola  peninsula  of  Russia,  and  the  array  on  Spitsbergen.  The  analysts 
had  waveforms  available  from  these  two  arrays  from  1 1  December  1992.  On  14  December 
1992,  phase  detections  from  the  Apatity  array  were  included  in  the  automatic  phase  asso¬ 
ciation.  The  phase  detections  from  the  Spitsbergen  array  were  made  available  to  the  ana- 
,  lysts  on  5  February  1993.  These  detections  are  not  used  in  the  automatic  phase 

association,  but  can  be  added  manually  during  analysis. 

The  operational  stability  of  IMS  has  been  very  good  during  the  reporting  period.  In  fact 
the  IMS  event  processor  (pipeline)  has  had  no  downtime  of  its  own;  i.e.,  all  data  available 
^  to  IMS  have  been  processed  by  IMS. 

Phase  and  event  statistics 

Table  3.5.1  gives  a  sumary  of  phase  detections  and  events  declared  by  IMS.  From  top  to 
bottom  the  table  gives  the  total  number  of  detections  by  the  IMS,  the  number  of  detections 
that  are  associated  with  events  automatically  declared  by  the  IMS,  the  number  of  detec¬ 
tions  that  are  not  associated  with  any  events,  the  number  of  events  automatically  declared 
by  the  IMS,  the  total  number  of  events  defined  by  the  analyst,  and  finally  the  number  of 
events  accepted  by  the  analyst  without  any  changes  (i.e.,  from  the  set  of  events  automati¬ 
cally  declared  by  the  IMS). 
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Oct  92 

Nov  92 

Dec  92 

Jan  93 

Feb  93 

Mar  93 

Total 

Phase  detections 

38026 

26150 

39885 

40148 

49504 

50710 

206397 

•Associated  phases 

5835 

5439 

5936 

5838 

6127 

8923 

32263 

-  Unassociated  phases 

32191 

20711 

33949 

34310 

43377 

41787 

174134 

Events  automatically 
declared  by  IMS 

1742 

1671 

1785 

1755 

1765 

2554 

9530 

No.  of  events  defined  by 
the  analyst 

1275 

1359 

1136 

1141 

1430 

2125 

7191 

No.  of  events  accepted 
without  modifications 

1018 

880 

951 

954 

1105 

1623 

5513 

I 

Table  3.5.1.  IMS  phase  detections  and  event  summary 


U.  Baadshaug 
B.  Ferstad 
B.Kr.  Hokland 
L.B.  Loughran 


3.6  GBF  operation 

Events  automatically  located  by  GBF 

During  days  275, 1992,  through  090, 1993, 9280  local  and  regional  events  were  located  by 
GBF.  This  gives  an  average  of  51.0  events  per  processed  day  (182  days  processed).  70% 
of  these  events  are  within  300  km  of  nearest  station,  and  85%  of  these  events  are  within 
1000  km  of  the  nearest  station. 

70.8%  of  these  events  were  defined  by  2  regional  phases.  Teleseismic  phases  are  currently 
not  used  by  GBF.  86.7%  of  all  events  had  3  defining  phases  or  less. 

T.  Kvsema 
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4  Improvements  and  Modifications 

4.1  NORSAR 

NORSAR  data  acquisition 

No  modification  has  been  made  to  the  NORSAR  data  acquisition  system. 

The  data  are  recorded  on  a  30-hour  circular  disk  buffer  on  the  IBM  system,  and  archived 
onto  1/2  inch  magnetic  tapes.  In  addition  to  this,  the  data  are  now  regularly  transmitted  to 
a  SUN  system  for  recording  on  a  48-hour  circular  disk  buffer. 

NORSAR  detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the  IBM  4381  com¬ 
puter  during  this  reporting  period. 

Detection  statistics  are  given  in  section  2.  In  addition  to  the  detection  processing  done  on 
IBM,  the  dp  program  is  doing  regular  detection  processing  on  a  SUN  system,  using  the 
unix-based  circular  disk  buffer  (see  below).  A  detection  SNR  threshold  of  20.0  triggers 
automatic  saving  of  waveforms  into  CSS  3.0  data  files. 

NORSAR  event  processing 

There  have  been  no  changes  in  the  routine  processing  of  NORSAR  events,  using  the  IBM 
system. 

In  parallel  with  the  IBM  processing,  routine  event  processing  is  also  done  on  a  SUN  com¬ 
puter  using  the  “old”  IBM  time  delay  correction  data  base  that  has  been  converted  to 
SUN/UNIX.  The  automatic  solutions  produced  are  equal  to  or  better  than  the  old  system 
with  a  lower  false  alarm  rate.  Alert  messages  are  sent  to  USGS  for  events  above  magni¬ 
tude  5.5. 

NORSAR  refurbishment 

As  reported  earlier,  the  main  problem  in  this  refurbishment  is  to  find  24-bit  A/D  convert¬ 
ers  with  low  power  consumption.  We  have  received  evaluation  units  from  Refraction 
Technology,  Teledyne  Geotech  and  Science  Horizons.  During  the  reporting  period,  the 
available  units  have  been  extensively  tested,  and  they  appear  to  perform  satisfactorily. 
Details  on  the  refurbishment  effort  will  be  reported  on  separately. 

42  Regional  Arrays 

Detection  processing 

The  routine  detection  processing  of  the  arrays  is  running  satisfactorily  on  each  of  the 
array’s  SUN-3/280  acquisition  systems.  The  same  program  is  used  for  NORSAR,  NOR- 
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ESS,  ARCESS,  FINESA,  GERESS,  Apatity  and  Spitsbeigen,  but  with  different  “recipes”. 
The  beam  table  for  NORESS  and  ARCESS  is  found  in  NORSAR  Sci.  Rep.  No.  1-89/90. 
The  beam  table  for  FINESA  and  GERESS  is  found  in  NORSAR  Sci.  Rep.  No.  1-90/91. 
The  beam  table  for  Apatity  is  found  in  NORSAR  Sci.  Rep.  No.  1-92/93,  and  that  for  Spits¬ 
bergen  is  given  in  the  present  report  (section  7.3). 

Detection  statistics  are  summarized  in  section  3. 

Signal  processing.  Phase  estimation 

This  process  performs  f-k  and  polarization  analysis  for  each  detection  to  determine  phase 
velocity,  azimuth  and  type  of  phase,  and  the  results  are  put  into  the  ORACLE  detection 
data  base  for  use  by  the  IMS. 

Event  Processing.  Plot  and  epicenter  determination 

A  description  of  single-array  event  processing  is  found  in  NORSAR  Sci.  Rep.  No.  2-88/ 
89,  and  NORSAR  Sci.  Rep.  No.  2-89/90. 


J.  Fyen 
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5  Maintenance  Activities 

5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  at  the  Maintenance  Center  (NMC)  Hamar,  and 
NDPC  activities  related  to  monitoring  and  control  of  NORSAR,  including  monitoring  of 
NORESS,  ARCESS,  FINESA,  GERESS  and  Apatity. 

Activities  involve  preventive  and  corrective  maintenance,  planning  and  installation  of  the 
Spitsbergen  array  (28  Oct  -  7  Nov  92)  and  completion  of  the  Apatity  installation  (22  Sep  - 
3  Oct  92). 

NORSAR 

Visits  to  subarrays  in  connection  with: 

•  Test  of  SLEM  equipment  in  order  to  locate  a  failure 

•  Reinstallation  of  SLEM  used  for  test  purposes 

•  Adjustment  of  gain  SP/LP  channels 

•  Adjustment  of  MP/FP 

•  Adjustment  of  DC  offset  SP  channels 

•  Replacement  of  MP/FP  motors  incl.  adjustments 

•  SLEM  rest 

•  Replacement  of  BE-card 

•  Comm,  cable  test 

•  Modem  and  line  test 


NMC 

•  Preparations  for  the  Spitsbergen  installation 

•  NORSAR  refurbishment  preparations 

•  Other  tasks  related  to  repair  of  SLEM  and  related  subarray  equipment 
NORESS 

•  Repair  of  C6,  D6 
Spitsbergen 

•  Installation  of  the  Spitsbergen  array  (Oct/Nov  92) 

•  Installation  of  a  second  windmill  and  extra  batteries  (Feb  93) 

Apatity 

•  Installation  continued  through  30  Sep  92  and  was  completed  3  Oct  92 
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Subarray/ 


NORSAR 

01B 

Installation  of  OlA  analog  unit  in  01B  SLEM  for  test 
purposes  (possible  because  OlA  was  still  down) 

6  Oct 

01B 

Spikes  in  data  still  required  CTV  engagement  A  test  with  a 
data  scope  at  die  NDPC  and  the  CTV  modem  (in  C-loop) 
proved  that  the  communications  system  was  OK.  A  Digital 
Unit  from  the  01B  SLEM  was  installed  in  the  OlA  SLEM  for 
continued  testing 

12  Oct 

01B 

Visit  to  the  SA  in  connection  with  installation  of  original 
SLEM  with  a  new  power  unit  Power  unit  failed  (lack  of +15 
Volt).  No  spare  unit  was  available. 

15  Oct 

OlA 

Reinstalled  the  SLEM  used  for  test  purposes  on  01B. 

Attempts  to  adjust  the  gain  SP/LP  resulted  in  improper  data. 

15  Oct 

OlA 

A  new  visit  to  the  subarray  resulted  in  proper  data  and 
demasking  of  the  array. 

22  Oct 

NPDC 

Daily  checks  of  the  following  arrays  have  been  carried  out 
i.e„  NORSAR,  NORESS,  ARCESS,  FINESA  and  Apatity. 
SP/LP  instruments  have  been  calibrated.  Free  Period  (FP)  and 
Mass  Pos.  (MP)  were  measured.  Those  outside  specifications 
adjusted  (when  feasible  from  NDPC). 

Oct 

NORSAR 

01B 

SLEM  brought  to  NMC  for  repair 

9  Nov 

01B 

SLEM  installed 

11  Nov 

02B 

Found  communication  inoperative 

20  Nov 

04C 

Replaced  BE-card  channel  1.  Adjusted  gain  all  SP/LP  chan¬ 
nels.  Adjusted  FP/MP  all  LP  instruments 

30  Nov 

NDPC 

Daily  checks  of  the  following  arrays  have  been  carried  out 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINESA,  Apatity  and 
Spitsbergen  (partly). 

SP/LP  instruments  were  calibrated  weeks  45, 47  and  48.  FP/ 
MP  were  measured.  Those  outside  specifications  adjusted 
(when  feasible  from  NDPC). 

Nov 
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Subarray/ 

area 

Task 

Date 

NMC 

The  staff  have  been  engaged  in  preparations  in  connection 
with  the  refurbishment  of  the  NORSAR  array.  In  addition, 
installation  of  the  Spitsbergen  array  (28  Oct  -  7  Nov  has  been 
carried  out 

Nov 

NORSAR 

03C 

Adjusted  channel  gain  all  SP/LP  channels.  Adjusted  MP/FP 

V,  NS,  EW  seismometers.  Adjusted  DC  offset  all  SP  chan¬ 
nels. 

2  Dec 

04 C 

LP  seismometer  NS  repaired 

4  Dec 

02B(tel) 

Replaced  the  telemetry  receiver  station  batteries. 

4  Dec 

NMC 

The  NMC  staff  continued  with  preparations  for  the  refurbish¬ 
ment  of  the  NORSAR  array. 

Dec 

NDPC 

Daily  checks  of  the  following  arrays  have  been  carried  out 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINESA,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  weeks  45, 49  -  52.  FP/  MP 
were  measured.  Those  outside  specifications  adjusted  (when 
feasible  from  NDPC). 

Dec 

NORSAR 

1993 

02B 

A  SLEM  reset  carried  out 

11  Jan 

NMC 

NORSAR  refurbishment  preparations  continued. 

Jan 

NDPC 

Daily  checks  of  the  following  arrays  have  been  carried  out 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINESA,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR).  FP /  MP  were 
measured.  Those  outside  specifications  adjusted  (when  feasi¬ 
ble  from  NDPC). 

Jan 

Spitsbergen  A  second  windmill  and  extra  batteries  were  mounted. 
NMC  NORSAR  refurbishment  preparations  continued. 


19  Feb 
Feb 
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Subarray/ 

area 

Task 

Date 

NDPC 

Daily  checks  of  the  following  arrays  have  been  carried  out, 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINESA,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR).  FP/  MP  were 
measured.  Those  outside  specifications  adjusted  (when  feasi¬ 
ble  from  NDPC). 

NORSAR 

02C 

Replaced  NS  MP  motor.  Adjusted  EW.NS  Mass  Position  and 
Free  Period  (MP/FP). 

11  Mar 

02C 

Replaced  EW  MP/FP  motor  and  adjusted. 

29  Mar 

NORESS 

Replaced  C6. 

23  Mar 

Final  repair  of  D6.  NORSAR  refurbishment 

31  Mar 

NPDC 

Daily  checks  of  the  following  arrays  have  been  carried  out, 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINESA,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR).  FP/  MP  were 
measured.  Those  outside  specifications  adjusted  (when  feasi¬ 
ble  from  NDPC). 

Table  5.1.  Activities  in  the  field  and  the  NORSAR  Maintenance  Center,  including  NDPC 
activities  related  to  NORSAR,  NORESS,  ARCESS,  FINESA,  GERESS,  Apatity  and 
Spitsbergen  1  Oct  92-31  Mar  93. 


5.2  Array  status 

As  of  31  March  1992  the  following  NORSAR  channels  deviated  from  tolerances: 


01A01 

8  Hz  filter 

02 

8  Hz  filter 

04 

30  dB  attenuation 

02B01 

08 

09 

02C01 

08 

03C04 

O.A.  Hansen 
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7  Summary  of  Technical  Reports  /  Papers  Published 

7.1  Intelligent  post- processing  of  seismic  events  -•  Part  1 : 

Bade  approach 

Introduction 

This  is  the  first  in  a  series  of  three  contributions  in  this  report  addressing  the  topic  of  intel¬ 
ligent  post-processing  of  seismic  events.  In  this  first  contribution  we  discuss  how  to  subdi¬ 
vide  the  area  to  be  monitored  in  order  to  identify  sites  of  particularly  high  seismic  activity. 
We  further  introduce  the  basic  idea  behind  this  post-processing  technique,  which  is  to  use 
as  a  starting  point  the  initial  event  location  provided  by  the  Intelligent  Monitoring  System 
(IMS)  and  then  use  region-specific  information  to  refine  the  solution.  By  applying  this 
technique  to  areas  with  significant  recurring  seismic  activity,  such  as  mining  sites,  a  con¬ 
siderable  part  of  the  analyst  work  can  be  eliminated. 

Since  15  October  1991,  the  Intelligent  Monitoring  System  (IMS,  Bache  et  al,  1993)  has 
been  processing  seismic  data  from  four  high-frequency  arrays  in  northern  Europe.  These 
are  NORESS  and  ARCESS  in  Norway,  FINESA  in  Finland  and  GERESS  in  Germany. 
During  October  1992  a  small-aperture  array  was  installed  near  Apatity  on  the  Kola  penin¬ 
sula,  and  during  late  October/early  November  1992  another  small-aperture  array  on  the 
island  of  Spitsbergen  became  operational.  The  data  from  these  installations  are  now 
included  in  the  IMS  processing  for  the  production  of  the  event  bulletin. 

Since  four  of  the  arrays  providing  data  to  the  IMS  are  located  in  Fennoscandia,  see  Fig. 
7.1.1,  the  IMS  event  bulletin  shows  an  excellent  event  detection  capability  for  this  region. 
Ringdal  (1991 )  found  that  for  Fennoscandia/NW  Russia,  a  network  consisting  of  NOR¬ 
ESS,  ARCESS  and  FINESA  has  a  90%  detection  capability  close  to  ML  2.0.  Near  the 
individual  arrays,  the  detection  capability  is  considerably  better,  and  consequently  a  large 
number  of  events  less  than  ML  1.0  are  detected. 

IMS  event  statistics 

The  basic  principle  of  the  post-processing  method  is  to  start  by  subdividing  the  area  to  be 
monitored  into  smaller  areas,  and  subsequently  apply  region-specific  analysis  to  each  such 
area.  As  an  example,  we  will  consider  in  some  detail  the  statistics  oi  events  in  Fennoscan¬ 
dia  and  NW  Russia  for  the  18-month  time  period  10/15/91  -  04/15/93.  We  will  only  con¬ 
sider  “well-defined”  events;  thus  we  ignore  events  with  author  identification  “yes/no”  and 
“ESAL/Poor_Loc”  in  the  origin  table. 

For  the  time  period  10/15/91  to  04/15/93  the  IMS  bulletin  contains  19503  well-defined 
events.  65.6%  (12799)  of  these  events  are  located  in  the  Fennoscandian/NW  Russian 
region  defined  by  the  map  of  Fig.  7.1.2, 15.8%  (3089)  are  located  within  5  degrees  of  the 
GERESS  array,  and  the  remaining  18.5%  (3615)  are  distributed  around  the  rest  of  the 
world,  mostly  at  teleseismic  distances  from  the  regional  array  network. 
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Figs.  7. 1.2-7. 1.4  show  the  event  distribution  in  Fennoscandia  for  all  magnitudes,  ML  >  1 
and  Ml  >  2,  respectively.  In  each  figure,  we  have  marked  the  approximate  geographical 
extent  of  8  main  mining  areas.  Table  7.1.1  lists  these  mining  sites  and  gives  details  on  the 
number  and  percentages  of  events  associated  to  the  sites  at  various  magnitudes. 

From  the  three  figures  and  Table  7.1.1  we  can  make  the  following  general  observations: 

•  Out  of  the  total  12799  events,  6317  (49.4  per  cent)  are  above  ML  =  1 ,  and  only  1131 
(8.8  per  cent)  are  above  ML  2. 

•  The  total  percentage  of  events  associated  with  the  8  mining  sites  is  47.88%  (all  mag¬ 
nitudes),  56.66%  (Ml  >  1.0)  and  65.61%  (Ml  >2.0).  Thus,  these  sites  become  more 
dominant  for  the  largest  events,  in  terms  of  relative  number  of  events  reported. 

•  Some  mining  sites  have  a  relatively  high  proportion  of  latge  events  (Ml  >2.0).  This 
is  particularly  noticeable  for  the  mining  areas  in  Western  Russia/Estonia.  On  the 
other  hand,  the  Kiruna  mine  has  the  largest  number  of  events  altogether,  but  almost 
none  of  these  are  above  Ml  2. 

Being  based  on  about  1  1/2  year  of  data,  the  statistics  discussed  here  should  be  reasonably 
representative  for  the  situation  in  the  Fennoscandian/NW  Russia  region.  Thus,  analysis  of 
recurring  events  from  these  mining  areas  is  a  significant  workload  for  the  analyst.  An 
automatic  method  to  improve  the  automatic  analysis  so  as  to  obtain  location  precisions 
comparable  with  the  analysts’  results  would  be  a  significant  development.  In  this  and  the 
next  two  sections  of  this  report,  we  will  show  that  such  an  improvement  is  possible  for  a 
well-calibrated  mining  area  (the  Khibiny  Massif). 

General  outline  of  the  method 

Most  automatic  detection  processor  algorithms  work  without  any  a  priori  assumptions  as 
to  when  and  where  a  seismic  event  occurred.  This  is,  of  course,  quite  reasonable,  and  to 
some  extent  inevitable.  The  detector  (SigPro)  associated  with  the  IMS  works  in  this  way. 
As  a  result,  some  of  the  SigPro  output  parameters,  which  later  will  be  used  by  the  IMS 
ESAL  system,  are  less  than  optimum. 

However,  once  an  initial  event  location  is  given  by  the  IMS,  it  is  possible  to  use  this  initial 
location  successively  in  an  automatic  iteration  scheme.  Each  iteration  gives  a  more  precise 
location,  which  in  turn  allows  the  automatic  program  to  place  successively  stronger  con¬ 
straints  on  the  processing  parameters. 

As  a  first  example  (see  Section  7.2),  we  can  consider  the  estimation  of  signal  arrival  time. 
Given  that  an  event  has  occurred  in  a  certain  area,  the  automatic  program  can  select  a  set 
of  optimum  filter  bands  and  beam  parameters  for  this  area,  prior  to  reassessing  the  arrival 
time  estimate.  As  shown  in  Section  7.2,  this  can  lead  to  a  remarkable  improvement  in  tim¬ 
ing  precision.  The  examples  given  in  Section  7.2  make  use  of  an  autoregressive  likelihood 
technique  (Pisarenko  et  a],  1987;  Kushnir  et  al,  1990).  It  is  noteworthy  that  this  method 
seems  to  require  that  the  search  be  limited  to  a  relatively  short  time  window  in  order  to 
work  well.  If  an  initial  location  and  origin  time  is  known,  we  can  obtain  the  required  short 
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data  in  a  low  passband,  e.g.  1  -  2  Hz,  where  S  has  a  low  SNR,  and  the  first  S-onset  is  in 
manv  rjwt  missed  With  this  nrohlem  in  mind  we  have  found  it  difficult  to  iustifv  a  com- 
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time  window  for  the  search.  The  method  is  therefore  well  suited  to  a  post-processing 
application. 

Another  example  is  the  estimation  of  azimuth  from  either  arrays  or  three-component  sta¬ 
tions.  The  advantage  of  using  a  fixed  frequency  interval  for  broadband  F-k  analysis  was 
convincingly  demonstrated  by  Kvaema  and  Ringdal  (1986).  Again,  a  prerequisite  was  the 
knowledge  that  the  event  in  question  was  located  in  a  certain  known  area. 

Section  7.3  demonstrates  that  the  approach  of  doing  post-processing  based  on  IMS  initial 
solutions  has  the  potential  of  providing  an  order-of-magnitude  improvement  in  location 
precision,  at  least  in  certain  cases  such  as  the  Khibiny  Massif  near  the  Apatity  array.  The 
improvement  may  be  less  if  no  network  station  is  located  close  to  the  source,  but  it  should 
still  be  significant.  For  example,  the  data  from  Kvaema  and  Ringdal  (1986)  indicate  that  a 
single  array  (NORESS)  would  be  capable  of  locating  the  Blisjo  explosions  to  within  an 
accuracy  of  10  km  or  better  at  a  distance  of  300  km.  This  is  compared  to  the  typical  uncer¬ 
tainty  of  about  30  km  in  traditional  single-array  location  estimates  at  this  distance 
(Mykkeltveit  and  Ringdal,  1981). 

In  general,  it  is  true  that  regional  corrections  are  required  in  order  to  compute  an  optimum 
location.  Again,  the  post-processing  analysis  is  well  suited  toward  this  end,  because  the 
corrections  can  be  tied  to  the  general  area,  to  which  the  initial  IMS  processing  assigns  the 
event. 

In  this  context,  it  is  important  to  note  that  no  regional  travel-time  tables  need  to  be 
involved  as  long  as  an  adequate  set  of  calibration  events  for  the  general  area  are  available. 
The  corrections  for  systematic  bias  may  be  made  both  to  the  phase  arrival  times  and  to  the 
estimated  azimuths.  Again,  this  subject  will  be  discussed  in  detail  in  Section  7.3,  in  con¬ 
nection  with  an  applicaiton  of  the  method  to  the  Khibiny  Massif  area  in  the  Kola  penin¬ 
sula. 

T.  Kvserna 

U.  Baadshaug 
F.  Ringdal 
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Region 

All  Mag 

mmm 

Number 

Percent 

Number 

Percent 

Number 

Percent 

Fte»»sc«»dia/NW 

12799 

100.00 

6317 

49.36 

mum 

8.84 

Estonia 

1487 

11.62 

1159 

18.36 

225 

19.89 

Karelia 

379 

2.96 

212 

3.36 

70 

6.19 

Khibiny 

1374 

10.74 

1106 

17.51 

233 

20.60 

Kiruna 

1953 

15.26 

634 

10.04 

11 

0.97 

Kostomuksha 

69 

0.54 

69 

1.09 

47 

4.16 

Kovdor 

112 

0.88 

99 

1.57 

34 

3.01 

Nikel 

620 

4.84 

181 

2.87 

104 

9.20 

Siilinjaervi 

134 

1.05 

119 

1.88 

18 

1.59 

Tbtalfbr8mines 

6128 

47.88 

3579 

56.66 

742 

65.61 

Table  7.1.1.  Distribution  of  events  in  mining  regions  of  Fennoscandia  and  NW  Russia. 
Events  with  author  identificantion  “yes/no”  and  “ESAL/Poor_Loc”  in  the  origin 
table  are  not  included  in  the  statistics. 
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Fig.  7.1.1.  Map  showing  the  location  of  the  regional  arrays  currently  used  by  the  Intelli¬ 
gent  Monitoring  System  in  operation  at  the  NORSAR  processing  center. 


12799  IMS  events 


Fennoscandia/NW  Russia 
10/15/91  -  4/15/93 

0*  5*  10*  15*  20*  25*  30*  35*  40* 


0*  5“  10*  15”  20’  25’  30”  35”  40’ 


Fig.  7.1.2.  Location  of  12,799  events  (all  magnitudes)  processed  by  the  IMS  for  an  18- 
month  period.  Only  event  solutions  of  satisfactory  quality  have  been  included  (see 
text  for  details).  Note  the  concentration  of  events  in  selected  mining  areas. 


6317  IMS  events 

Fennoscandia/NW  Russia  (ml  >1.0) 
10/15/91  -  4/15/93 


O'  5*  1(T  15*  20‘  25’  30*  35*  40* 


Fig.  7.1.3.  Same  as  Fig.  7.1.2,  but  showing  only  events  of  ML  >  1.0.  The  total  number  of 
events  is  6317  for  the  18-month  period. 
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1 131  IMS  events 

Fennoscandia/NW  Russia  (ml  >  2.0) 
10/15/91  -  4/15/93 


Fig.  7.1.4.  Same  as  Fig.  7.1.2,  but  showing  only  events  of  Ml  >  2.0.  The  total  number  of 
events  is  1131  for  the  18-month  period. 
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7 2  Intelligent  post-processing  of  seismic  events  -  Part  2: 

Accurate  determination  of  phase  arrival  times  using  autoregessive 
likelihood  estimation 

Introduction 

A  precise  estimation  of  the  onset  time  of  seismic  phases  is  necessary  to  obtain  an  accurate 
event  location.  In  the  context  of  automatic  processing  of  the  regional  arrays  operated  by 
NORSAR,  a  two-step  procedure  has  been  in  use  since  the  first  regional  array,  NORESS, 
was  established  in  1985  (Mykkeltveit  and  Bungum,  1984).  This  procedure  consist  of  first 
applying  a  series  of  STA/LTA  detectors  in  parallel  to  a  set  of  filtered  beams.  When  one  or 
more  of  the  STA/LTA  detectors  exceed  a  predefined  threshold,  a  phase  detection  is 
declared  and  a  detection  time  is  found.  Subsequently,  a  time  domain  phase  timing  algo¬ 
rithm  is  applied  to  the  filtered  beam  with  the  highest  SNR,  using  the  detection  time  as  the 
starting  value.  A  detailed  description  of  this  algorithm  is  found  in  Mykkeltveit  and  Bun- 
gum  (1984). 

These  estimates  of  the  onset  times  are  used  by  the  automatic  phase  association  and  event 
location  procedure  (ESAL)  of  the  Intelligent  Monitoring  System  (IMS)  (Bache  et  al, 
1993),  and  the  fully  automatic  processing  results  are  finally  reviewed  and  corrected  by  the 
analyst  using  the  Analyst  Review  Station  (ARS)  of  the  IMS.  Through  the  analyst  review 
we  have  experienced  that  in  most  cases  the  phase  onsets  need  to  be  adjusted.  In  section  7.6 
of  this  report,  comparisons  between  the  automatic  and  manual  onset  time  estimates  of  var¬ 
ious  seismic  phases  are  presented.  For  P-type  phases  the  standard  deviations  of  automatic 
-  manual  onset  time  generally  varies  between  0.5  and  0.7  seconds.  It  should  be  noticed  that 
these  estimates  represent  averages  for  all  seismic  phases  for  a  one-week  period,  and  thus 
include  phases  with  very  different  SNR’s,  frequency  characteristics  and  signatures.  For 
phases  with  an  high  SNR  and  an  instantaneous  signature,  the  performance  of  the  phase 
timing  algorithm  is  significantly  better.  Nevertheless,  one  of  the  conclusions  from  the 
study  presented  in  section  7.6  is  that  in  order  to  improve  the  precision  of  the  automatic 
event  locations  provided  by  the  IMS  and  in  order  to  reduce  the  analyst’s  workload,  there  is 
a  strong  need  to  improve  the  precision  of  the  automatic  onset  time  estimates. 

In  this  study  we  will  investigate  the  potential  automatic  use  of  an  onset  picker  based  on 
autoregressive  likelihood  estimation  (Pisarenko  et  al,  1987,  Kushnir  et  al,  1990).  Both  a 
single-component  version  (ESTON1)  and  a  three-component  version  (ESTON3)  of  this 
method  will  be  tested  on  data  from  events  located  in  the  Khibiny  Massif  on  the  Kola  Pen¬ 
insula  of  Russia,  recorded  at  the  Apatity  array,  the  Apatity  three-component  station  and 
the  ARCESS  array,  see  Fig.  7.2.1. 

The  Khibiny  Massif  events 

Six  apatite  mines  are  located  within  an  area  of  about  10  km2  in  the  Khibiny  Massif  (see 
Fig.  7.2.1). A  detailed  description  of  these  mines  and  the  mining  activity  is  found  in 
Mykkeltveit  (1992).  The  coordinates  of  the  mines  are  given  in  Table  7.2.1.  Notice  that 
Mine  I  consists  of  both  an  underground  mine  and  an  open-pit  mine.  Although  we  have  no 
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explicit  information  on  the  exact  sizes  of  these  mines,  interpretation  of  various  maps  sug¬ 
gest  that  the  typical  size  is  about  1  km2.  The  Kola  Regional  Seismological  Centre  has 
since  the  beginning  of  1991  provided  NORSAR  with  information  on  mining  blasts  in  the 
six  Khibiny  mines.  The  information  provided  contains  an  assignment  of  the  relevant  mine 
(I-VI),  P  (and  normally  also  S)  arrival  times  at  the  analog  APA  (Apatity)  station  (co¬ 
located  with  the  three-component  station),  the  amplitude  and  period  of  the  signal,  and  the 
total  charge  size  in  tons.  Detailed  information  on  the  58  events  used  in  this  study  is  given 
in  Table  12.2. 

At  the  Apatity  array  (APAO)  and  the  Apatity  three-component  station  (APZ9)  which  are 
located  within  a  distance  range  of  18  -  49  km  from  the  different  Khibiny  mines,  the  seis¬ 
mograms  show  clear  P  (Pg),  S  (Sg  or  Lg)  and  Rg  arrivals  (see  Figs.  7.2.2  and  7.2.3).  At 
ARCESS  which  is  located  about  400  km  from  the  mines,  a  clear  Pn  and  emergent  Pg,  Sn 
and  Lg  phases  are  observed  (see  Fig.  7.2.4).  Table  7.2.3  gives  detailed  information  on  the 
distances  and  azimuths  from  APZ9,  APAO  and  ARCESS  to  the  six  Khibiny  mines. 

Estimation  of  P  arrival  time 

We  will  not  go  into  any  detail  concerning  the  theory  of  the  autoregressive  likelihood  tech¬ 
nique  for  onset  time  estimation,  but  refer  to  Pisarenko  et  al  (1987)  and  Kushnir  et  al 
(1990)  for  details.  We  will,  however,  concentrate  on  the  practical  aspects  of  implementing 
the  method  as  part  of  an  automatic  processing  sequence.  The  onset  time  procedure 
requires  that  we  have  available  an  approximate  timing  of  the  phase  arrival,  and  the  search 
for  the  exact  onset  is  limited  to  an  interval  around  the  initial  arrival  time.  The  initial  arrival 
time  can  be  obtained  in  several  ways: 

•  Pisarenko  et  al  (1987)  and  Kushnir  et  al  ( 1990)  suggest  that  an  optimal  event  detec¬ 
tor  based  on  an  Bayesian  approach  using  autoregressive  modelling  of  the  data 
should  be  used  for  phase  detection  and  approximate  timing. 

•  Another  alternative,  which  we  will  use  for  P-phases  at  APAO  and  ARCESS,  is  to 
take  the  onset  time  provided  by  the  IMS  as  the  initial  estimate.  This  onset  time  is 
calculated  using  the  algorithm  described  in  the  Introduction. 

•  A  third  alternative  is  to  use  the  predicted  phase  arrival  time,  e.g.,  derived  from  an 
initial  event  location  and  origin  time  or  from  the  expected  pattern  of  phase  arrivals 
from  events  in  a  given  region.  This  is  the  approach  to  be  used  for  both  P-  and  S- 
phases  at  APZ9  and  for  S-phases  at  APAO. 

The  autoregressive  likelihood  estimator  is  based  on  regarding  the  signal  onset  time  as  the 
moment  in  time  when  the  statistical  features  of  the  observed  time  series  (single-compo¬ 
nent  or  three-component)  are  abruptly  changed.  The  single-component  method  thus  takes 
into  account  changes  in  both  power  and  frequency  content,  and  it  is  therefore  important 
that  the  broadband  signal  waveforms  are  retained  (no  narrow  bandpass  filtering).  In  addi¬ 
tion  to  taking  into  account  changes  in  power  and  frequency  content,  the  three-component 
method  is  also  sensitive  to  changes  in  the  polarization  characteristics  of  the  three-compo¬ 
nent  data.  From  experimenting  with  both  ESTON1  and  ESTON3  we  have  found  that  in 
order  to  obtain  stable  estimates,  the  data  should  first  be  lowpass  filtered  and  decimated 
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according  to  the  highest  frequencies  of  the  signals.  For  P-phases  from  the  Khibiny  events 
recorded  at  the  Apatity  stations  and  at  ARCESS,  signal  frequencies  up  to  the  Nyquist  fre¬ 
quency  of  20  Hz  are  in  almost  all  cases  observed.  Consequently,  no  lowpass  filtering  or 
decimation  is  applied  to  the  data.  However,  a  low-order  prediction  error  filter,  designed 
from  a  25  second  noise  sample  preceding  the  P-phase,  is  applied. 

For  the  arrays  APAO  and  ARCESS,  the  slowness  and  azimuth  of  the  P-phases  are  com¬ 
puted  using  broad-band  f-k  analysis  (Kvaema  and  Doombos,  1986).  Kvaema  and  Ringdal 
(1986)  showed  that  the  key  to  obtain  stable  slowness  and  azimuth  estimates  of  phases 
from  a  given  region  is  to  process  the  data  in  a  fixed  frequency  band,  using  a  fixed  time 
window  positioning.  The  time  window  positioning  for  the  f-k  analysis  and  the  fact  that  the 
onset  time  is  computed  from  a  steered  beam  (with  steering  delays  found  from  f-k  analysis) 
make  these  two  processes  to  be  closely  integrated.  Through  extensive  testing  we  have 
found  that  for  P-phases  from  the  Khibiny  mines  recorded  at  APAO,  the  most  stable  slow¬ 
ness  and  azimuth  estimates  are  obtained  if  the  f-k  spectrum  is  computed  in  the  frequency 
band  4.0  - 10.0  Hz  using  the  9  SPZ-component  sensors  of  the  array.  The  time  window  has 
a  length  of  2  seconds  and  starts  0.3  seconds  ahead  of  the  P-arrival.  For  ARCESS  the  most 
stable  estimates  are  obtained  if  all  25  SPZ-component  sensors  are  processed  in  the  fre¬ 
quency  band  3.0  -  6.0  Hz  using  a  3  second  time  window  starting  0.3  seconds  ahead  of  the 
P-arrival.  Further  refinement  can  be  made  by  including  a  procedure  that  checks  the  SNR  in 
the  predefined  frequency  band,  and  adjusts  the  frequency  cutoffs  in  the  case  of  low  SNR. 

For  APAO  and  ARCESS,  the  automatic  algorithm  used  for  estimating  the  P-onsets  from 
events  from  the  Khibiny  Massif  can  be  summarized  in  the  following  way: 

•  Use  the  P  onset  from  the  IMS  as  the  initial  arrival  time  estimate. 

•  Generate  a  beam  with  steering  delays  corresponding  to  the  IMS  f-k  results  and  pre¬ 
whiten  the  beam  with  a  low  order  prediction  error  filter. 

•  Apply  ESTON1  to  the  prewhitened  beam  in  an  interval  of  ±2  seconds  around  the 
initial  onset  estimate.  A  sliding  window  length  of  1  second  is  used  by  ESTON1,  and 
the  autoregressive  modelling  is  of  order  3. 

•  Position  the  time  window  according  to  the  onset  time  from  ESTON1  and  run  the  f-k 
analysis  as  oudined  above. 

•  Generate  a  new  beam  with  steering  delays  corresponding  to  the  new  f-k  results  and 
prewhiten  the  beam  with  a  low  order  prediction  error  filter. 

•  Apply  ESTON1  to  the  new  prewhitened  beam  in  an  interval  of  ±2  seconds  around 
the  previous  ESTON1  onset  estimate. 

•  Position  the  time  window  according  to  the  last  onset  time  estimate  from  ESTON1 
and  run  the  f-k  analysis  once  more  as  outlined  above. 

For  the  three-component  station  APZ9  the  situation  is  somewhat  different  as  this  station  is 
not  part  of  the  automatic  IMS  processing.  APZ9  is,  however,  located  only  18  km  from 
APAO,  and  we  therefore  choose  to  use  the  onset  time  at  APAO  as  the  initial  arrival  time 
estimate.  The  processing  can  be  summarized  as  follows: 

•  Use  the  P  onset  at  APAO  as  the  initial  arrival  time  estimate. 
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•  Prewhiten  the  Z-component  APZ9  with  a  low  order  prediction  error  filter. 

•  Apply  ESTON 1  to  the  prewhitened  Z-component  APZ9  in  an  interval  of  ±5  seconds 
around  the  initial  onset  estimate  from  APAO.  The  sliding  window  length  used  by 
ESTON  1  is  1  second  and  the  autoregressive  modelling  is  of  order  3. 

•  Reestimate  the  onset  time  applying  ESTON  1  to  the  prewhitened  z-component  APZ9 
in  an  interval  of  ±2  seconds  around  the  previous  ESTON  1  onset  estimate. 

Precision  of  automatic  P-onsets  using  the  single-component  estimator  (ESTON  1) 

The  precision  of  the  automatic  P-onsets  using  ESTON  1  has  been  assessed  by  two  different 
methods.  The  first  method  is  to  compare  the  ESTON  1  estimates  to  the  best  manual  pick. 
The  purpose  of  this  approach  is  to  obtain  information  on  any  bias  in  ESTON  1  estimates, 
and  also  to  check  the  consistency  between  the  automatic  and  manual  onset  estimates.  In 
the  second  method  we  estimate  the  standard  deviation  of  the  phase  picks  by  looking  at  the 
consistency  of  the  arrival  time  differences  between  phases  from  events  located  in  the  same 
mine,  see  Sereno  (1990).  In  this  way,  an  unbiased  estimate  of  the  standard  deviation  is 
found  for  both  the  automatic  and  manual  picks. 

In  Fig.  7.2.5.a  the  time  differences  between  the  automatic  onsets  from  IMS  and  the  manual 
pick  at  APAO  are  presented  as  a  function  of  signal-to-noise  ratio  (SNR)  on  the  prewhit¬ 
ened  beam.  The  standard  deviation  of  these  differences  is  as  high  as  0.43  seconds.  In  com¬ 
parison,  the  time  difference  between  the  automatic  ESTON  1  onsets  and  the  manually 
picked  onsets,  given  in  Fig.  7.2.5.b  has  a  standard  deviation  as  low  as  0.02  seconds,  and 
has  a  systematic  positive  bias  of  0.05  seconds,  i.e.  the  ESTON  1  onsets  are  consistently 
picked  a  bit  late. 

In  Fig.  7.2.6  similar  plots  are  presented  for  Pn  phases  recorded  at  ARCESS.  The  automatic 
onsets  from  IMS  shows  significantly  less  scatter  (o  *  0.13  s)  than  for  APAO.  On  the  other 
hand,  the  automatic  ESTON  1  onsets  has  a  somewhat  larger  scatter  (o  =  0.03  s)  and  a 
somewhat  larger  positive  bias  (0.08  s)  than  at  APAO.  Nevertheless,  the  accuracy  of  the 
ESTON  1  time  picks  appears  to  be  significantly  better  than  those  of  the  current  processing 
system. 

To  address  possible  dependency  of  ESTON  1  on  the  signal-to-noise  ratio,  we  have  in  Fig. 
7.2.7  plotted  the  time  difference  between  the  automatic  ESTON  1  onsets  and  the  manually 
picked  onsets  for  both  APAO  and  APZ9.  In  this  way,  we  get  an  overview  of  the  scatter  for 
SNR’s  ranging  from  about  3  to  about  500.  We  find  from  this  figure  that  the  time  differ¬ 
ences  seem  to  be  almost  independent  of  SNR,  which  again  indicate  that  the  ESTON1 
onset  estimator  works  well  for  quite  low  signal-to-noise  ratios.  For  direct  comparison,  we 
have  in  Fig.  7.2.8  plotted  the  ARCESS  data  on  the  same  scale,  but  for  ARCESS  no  low 
SNR  phases  are  observed,  as  all  SNR’s  exceed  6.  By  SNR,  we  mean  the  maximum  of  the 
linear  ratio  STA/LTA  (i.e.,  short  term  average  divided  by  long  term  average.) 

The  results  presented  above  show  that  an  automatic  onset  time  procedure  for  P-phases 
using  ESTON  1  can  give  a  remarkable  improvement  in  precision  compared  to  the  current 
algorithm  used  in  the  IMS.  The  ESTON  1  onset  estimator  has  a  bias  that  is  dependent  on 
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the  dominant  frequency  of  the  signal.  At  the  Apatity  stations  the  average  dominant  P  fre¬ 
quency  is  1 3  Hz  and  the  average  bias  is  0.05  s,  whereas  at  ARCESS  the  average  dominant 
P  frequency  is  6  Hz  and  the  average  bias  is  0.08  s.  The  difference  in  bias  can  not  be  due  to 
differences  in  sampling  rate,  since  both  APAO,  APZ9  and  ARCESS  have  a  sampling  fre¬ 
quency  of  40  Hz.  For  the  high  frequency  signals  at  the  Apatity  stations,  the  bias  shows  no 
clear  dependency  on  SNR,  at  least  not  for  SNR  above  3.  As  expected,  the  results  also  sug¬ 
gest  that  the  precision  of  ESTON1  increases  with  increasing  signal  frequencies.  This  can 
be  inferred  from  the  fact  that  the  automatic  P-onset  estimates  of  the  high-frequency  sig¬ 
nals  at  the  Apatity  stations  are  more  consistent  with  the  manual  picks  than  the  P-onsets  at 
ARCESS,  which  generally  have  a  lower  dominant  frequency.  According  to  Pisarenko  et  al 
(1987)  there  are  no  analytical  expressions  for  the  theoretical  biases  and  variances  of 
ESTON1  and  ESTON3,  and  it  is  therefore  necessary  to  obtain  empirical  values.  In  any 
case,  the  biases  and  variances  are  less  than  0.1  s  for  the  P-phases  considered  in  this  study. 

An  unbiased  estimate  of  the  measurement  variance  is  determined  from  the  arrival  time  dif¬ 
ference  between  two  phase  observations  for  repeated  events  in  the  same  mine. 
Specifically: 


*\ypick^u2,pick 


I  I  <aw,i: 

k  =  1  i  =  1 _ 

W obs  ~  ^ mines ) 


where  of  and  o\  are  the  picking  variances  of  each  phase,  A  7^  is  the  ith  observation  of 
the  arrival-time  difference  for  the  kth  mine.  <a Tobl)i  is  the  mean  arrival  time  difference  for 
the  kth  mine.  Noii  is  the  total  number  of  observations  (at  all  mines),  and  NMi„„  is  the  num¬ 
ber  of  mines. 

By  computing  the  arrival-time  differences  between  the  P  observations  at  the  three  stations 
APA0(1),  APZ9(2)  and  ARCESS(3),  we  get  three  equations  of  the  type  above,  with  alto¬ 
gether  three  unknowns.  These  equations  can  easily  be  solved  to  obtain  estimates  of  each 
individual  variance  value.  For  example,  for  the  automatic  picks  we  obtain: 

o*  +  o|  =  (0.065) 2 

o]  +  o\  =  (0.065) 2 

o\  +  o\  =  (0.074)2 

which  gives  o,  =  0.04,  o2  =  o.os  and  o3  =  o.os  (see  Table  7.2.4). 
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The  standard  deviations  of  both  the  manual  picks  and  the  automatic  picks  from  ESTON1 
are  given  in  Table  7.2.4.  All  standard  deviations  are  less  than  or  equal  to  0.06  s,  and  a  part 
of  this  variability  is  likely  due  to  the  fact  that  the  events  of  each  mine  are  not  located  at  the 
same  spot,  but  are  distributed  within  the  mine.  Under  the  assumption  that  each  mine  has 
an  extent  of  1  x  1  km  (which  is  reasonable  from  interpretation  from  various  maps),  we 
have  computed  the  maximum  theoretical  P  travel-time  difference  between  APAO  and 
APZ9  for  the  six  Khibiny  mines.  The  largest  value  is  obtained  for  Mine  I,  where  a  0.09  s 
travel-time  difference  is  possible.  The  smallest  values  are  found  for  Mine  V  and  Mine  VI, 
where  a  0.04  s  travel-time  difference  is  possible.  This  clearly  suggest  that  location  vari¬ 
ability  within  each  mine  can  have  a  significant  impact  on  the  estimates  of  the  picking  pre¬ 
cision.  If  we,  however,  assume  that  the  distribution  of  the  events  within  each  mine 
increases  all  picking  error  estimates  with  a  similar  amount,  a  likely  interpretation  of  the 
differences  in  picking  precision  between  the  three  stations  is  as  follows: 

•  The  P  picks  at  APZ9  are  less  precise  than  at  APAO  due  to  generally  lower  SNR. 

•  The  P  picks  at  ARCESS  are  less  precise  than  at  APAO  due  to  generally  lower  domi¬ 
nant  frequency  of  the  signals. 

The  results  of  Table  7.2.4  show  that  the  automatic  ESTON1  method  matches  the  human 
precision,  and  that  for  the  relatively  high  SNR  P  arrivals  analyzed  in  this  study,  the  stan¬ 
dard  deviation  of  the  automatic  picks  is  well  below  0.1  s. 

Estimation  of  S  arrival  time  at  the  Apatity  stations 

As  seen  from  Fig.  7.2.2  and  7.2.3,  the  seismograms  of  the  Khibiny  events  show  clear  P,  S 
and  Rg  arrivals  at  the  Apatity  stations  APAO  and  APZ9.  The  S-onsets  do,  however, 
become  more  emergent  with  increasing  source-receiver  distance.  The  S-phases  typically 
have  the  largest  SNR  on  the  transverse  component,  but  they  are  also  clearly  observed  on 
the  radial  and  vertical  components.  For  automatic  estimation  of  S-onsets,  we  have  experi¬ 
mented  with  both  the  single-component  ESTON1  onset  estimator  applied  to  the  transverse 
component,  and  the  three-component  ESTON3  onset  estimator  applied  to  the  three-com¬ 
ponent  data,  and  found  that  the  ESTON3  method  gave  the  best  results.  The  most  stable 
results  were  obtained  if  the  data  were  first  filtered  in  a  relatively  wide  band  between  2.0 
and  8.0  Hz  and  then  decimated  to  a  sampling  rate  of  20  Hz. 

The  Rg  arrivals  occur  very  close  in  time  to  the  S-onset  Moreover,  the  Rg  phase  is  disper¬ 
sive.  We  therefore  did  not  succeed  in  estimating  the  Rg-onset  in  a  reliable  way  with  the 
autoregressive  likelihood  estimation  technique.  In  order  to  design  an  automatic  processing 
sequence  for  the  Khibiny  events  at  the  Apatity  stations  APAO  and  APZ9,  we  have  utilized 
the  expected  pattern  of  phase  observations  from  events  in  this  region,  and  the  procedure  is 
as  follows: 

•  Estimate  P-onset  as  previously  outlined 

•  Identify  the  peak  of  the  Rg-phase  from  an  STA  envelope  created  from  the  filtered  z- 
component  (0.8  -  2.0  Hz).  The  search  interval  for  the  Rg  maximum  is  currently  lim¬ 
ited  to  20  seconds  after  the  P-onset. 
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•  For  the  APAO  array,  the  slowness  and  azimuth  of  the  Rg  phase  is  estimated  using 
data  in  a  5  s  window  starting  3  s  ahead  of  the  Rg  peak.  A  frequency  band  of  0.8  -  2.0 
Hz  is  used  in  the  f-k  analysis. 

•  For  estimation  of  S-onsets,  the  three-component  data  are  first  filtered  and  decimated 
as  outlined  above,  and  the  ESTON3  estimator  is  then  applied  within  a  time  interval 
that  starts  2  s  after  the  P-onset  and  stops  at  the  time  of  the  Rg  peak.  The  sliding  win¬ 
dow  length  used  by  ESTON3  is  2  seconds  and  the  autoregressive  modelling  is  of 
order  3. 

•  For  the  APAO  array,  the  slowness  and  azimuth  of  the  S  phase  is  estimated  using  data 
in  a  2  s  window  starting  0.3  s  ahead  of  the  S-onset.  A  frequency  band  of  2-  3  Hz  is 
used  in  the  f-k  analysis. 

Illustrations  of  the  automatic  processing  sequence  are  shown  if  Figs.  7.2.9  and  7.2.10. 
Notice  the  clear  peak  of  the  ESTON3  likelihood  function  at  the  S-onset  for  both  events. 

Precision  of  automatic  S-onsets  using  the  three-component  estimator  ( ESTO ’  3) 

When  comparing  the  automatic  ESTON3  onsets  to  the  manually  picked  S-onsets  at  APAO, 
we  find  that  the  standard  deviation  of  the  time  differences  is  0.18  s  and  that  the  bias  is  0.01 
s.  For  APZ9  the  standard  deviation  is  0.12  s  and  the  bias  is  0.0  ''  s.  Due  to  the  relatively 
large  scatter  in  the  observations,  this  bias  of  0.04  is  not  significantly  different  from  0. 

To  compute  the  measurement  variance  of  the  S-onset  estimates,  we  have  again  investi¬ 
gated  the  consistency  of  the  arrival  time  difference  between  two  phase  observations  for 
repeated  events  in  the  same  mine,  according  to  equation  7.2.1.  By  using  the  P-onsets  at 
both  APAO,  APZ9  and  ARCESS  as  references,  for  which  the  measurement  variances  are 
known  (see.  Table  7.2.4),  we  can  reliably  estimate  both  the  manual  and  the  automatic  S- 
onset  variances  at  both  APAO  and  APZ9.  The  results  are  summarized  in  Table  7.4.5 

We  see  from  Table  7.4.5  that  the  picking  uncertainty  estimates  based  on  all  three  reference 
P-phases  are  very  consistent,  indicating  that  the  method  for  estimating  uncertainty  works 
well.  The  fact  that  the  precision  of  the  manual  and  automatic  picks  are  about  equal,  indi¬ 
cates  that  the  automatic  S-onset  algorithm  using  ESTON3  matches  the  human  precision. 
Table  7.4.5  also  show  that  the  S-onsets  at  APZ9  (o  =  0.13  s)  are  generally  more  precise 
than  at  APAO  (o  =0.19  s).  We  believe  that  this  is  due  to  the  fact  that  the  S-phases  become 
more  emergent  with  increasing  source-receiver  distance,  as  illustrated  in  Figs.  7.2.2  and 
7.2.3.  The  six  Khibiny  mines  are  located  within  a  distance  range  of  18  -  33  km  from 
APZ9,  and  within  a  distance  range  of  32  -  49  km  from  APAO.  see  Table  7.2.3. 

It  would  have  been  interesting  to  compare  the  precision  of  the  automatic  S-onsets  from 
ESTON3  to  the  precision  of  the  automatic  S-onsets  used  by  the  IMS.  We  have,  however, 
experienced  that  the  continuous  processing  of  the  Apatity  array  data  (APAO)  has  problems 
in  detecting  and  estimating  the  onsets,  slowness  and  azimuths  of  secondary  phases  that 
have  little  time  separation  and  large  differences  in  frequency  content.  The  first  S-detection 
of  the  Khibiny  events  recorded  at  APAO  is  in  almost  all  cases  declared  in  a  low  frequency 
band,  which  is  typical  for  detection  of  Rg.  Consequently  the  onset  routine  prefilters  the 
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data  in  a  low  passband,  e.g.  1  -  2  Hz,  where  S  has  a  low  SNR,  and  the  first  S-onset  is  in 
many  cases  missed.  With  this  problem  in  mind  we  have  found  it  difficult  to  justify  a  com¬ 
parison  between  the  automatic  S-onsets  from  ESTON3  and  the  automatic  S-onsets  from 
the  procedure  providing  data  to  the  IMS. 

Precision  of  azimuth  estimates  from  broad-band  f-k  analysis 

As  described  above,  the  estimation  of  onset  time  and  the  estimation  of  slowness  and  azi¬ 
muth  by  f-k  analysis  are  closely  integrated  processes.  In  one-array  location  of  seismic 
events,  the  azimuth  estimates  are  necessary  to  be  able  to  compute  an  event  location,  and  in 
the  event  location  procedure  of  the  IMS  (Bratt  and  Bache,  1988)  the  phase  azimuth  esti¬ 
mates  are  required  to  be  accompanied  with  an  uncertainty  estimate. 

In  Table  7.5.6,  the  mean  and  the  standard  deviation  of  the  azimuth  residuals  relative  to  the 
Khibiny  mine  locations  of  the  phases  recorded  at  the  arrays  APAO  and  ARCESS  are  pre¬ 
sented.  We  see  from  the  table  that  the  P  and  Rg  azimuths  at  APAO  have  about  the  same 
standard  deviation  (3.9  and  3.4  degrees,  respectively),  but  that  the  P  azimuths  have  a  sys¬ 
tematic  bias  of  8.2  degrees.  The  Rg  azimuth  bias  is  as  low  as  -1 .8  degrees.  The  ARCESS  P 
azimuths  have  a  very  low  standard  deviation  (0.9  degrees),  but  a  systematic  bias  of  4.6  is 
consistently  observed.  The  S-azimuths  at  APAO  show  a  very  large  scatter.  This  show  that 
if  the  Khibiny  events  are  located  without  introducing  corrections  for  the  azimuth  biases, 
the  Rg  azimuths  should  be  given  the  smallest  a  priori  uncertainty.  If  the  systematic  biases 
are  removed,  the  a  priori  uncertainty  of  P  and  Rg  at  APAO  become  comparable.  With  the 
systematic  bias  removed,  the  a  priori  uncertainty  of  ARCESS  P  is  very  small,  but  it  should 
be  noticed  that  in  the  event  location  procedure  (Bratt  and  Bache,  1988)  the  a  priori  azi¬ 
muth  uncertainty  is  scaled  by  the  source-receiver  distance,  such  that  an  azimuth  observa¬ 
tion  at  400  km  distance  with  an  uncertainty  of  0.9  degrees  (ARCESS  P)  is  given  less 
weight  than  an  azimuth  observation  at  40  km  distance  with  an  uncertainty  of  3.9  degrees 
(APAO  P). 


Conclusions 

The  results  presented  in  this  study  show  that  very  precise  automatic  estimates  of  phase 
onsets  from  events  in  the  Khibiny  Massif  can  be  obtained  with  the  autoregressive  likeli¬ 
hood  estimation  method.  Implementation  of  the  method  requires  that  we  have  available 
approximate  estimates  of  the  phase  arrival,  and  we  have  shown  that  such  approximate  esti¬ 
mates  can  oe  obtained  from  the  IMS  event  definitions  (phase  association  and  event  loca¬ 
tion)  or  from  analysis  of  previous  events  in  the  region.  In  this  way,  the  autoregressive 
likelihood  estimation  method  can  provide  phase  onsets  that  match  the  human  precision. 
The  uncertainties  and  biases  of  the  automatic  onset  estimates  of  various  phases  at  the  Apa¬ 
tity  stations  and  at  ARCESS  have  been  quantified,  and  the  precision  of  the  automatic 
phase  picks  shows  very  large  improvement  in  comparison  to  the  automatic  phase  onsets 
from  the  continuous  processing  providing  input  to  the  IMS. 

In  addition  to  the  automatic  onset  estimation,  we  have  estimated  the  phase  azimuths  and 
slowness  using  oroad-band  f-k  analysis.  This  has  been  done  using  data  in  a  fixed  fre¬ 
quency  band,  using  a  fixed  window  positioning,  as  suggested  by  Kvaema  and  Ringdal 


< 


75 


(1986)  for  the  purpose  of  obtaining  increased  stability.  ForP-phases  recorded  at  ARCESS, 
the  results  are  in  accordance  with  those  of  Kvaema  and  Ringdal  ( 1986),  where  a  scatter  of 
only  ±1  degree  were  observed  at  NORESS  for  Pn  phases  from  repeated  events  at  a  dis¬ 
tance  of  300  km.  Comparing  to  the  overall  azimuth  uncertainty  of  phases  recorded  at 
ARCESS  and  NORESS  (Sereno,  1990),  we  find  that  if  the  systematic  biases  are  removed 
from  our  azimuth  estimates,  the  event  location  precision  can  be  significantly  improved. 
Without  introducing  azimuth  corrections,  the  Rg  azimuths  at  APAO  show  to  be  quite  reli¬ 
able. 

We  realize  that  in  order  to  obtain  accurate  event  locations,  precise  onset  time  and  azimuth 
estimates  are  necessary,  but  not  sufficient  If  the  theoretical  travel-time  model  used  in  the 
event  location  deviates  from  the  true  travel-times,  the  precision  of  the  event  locations  will 
be  reduced.  To  overcome  this,  we  will  in  Section  7.3  of  this  report  discuss  the  introduction 
of  travel-time  corrections. 

A  natural  extension  to  this  study  will  be  to  investigate  the  performance  of  ESTON1  and 
ESTON3  applied  to  phases  at  the  other  stations  of  the  network  recording  the  Khibiny 
events.  At  NORESS  and  FINESA  the  detected  Pn  phases  have  quite  low  SNR,  and  it 
would  be  interesting  to  quantify  the  precision  of  the  automatic  onset  estimates  of  these 
phases.  It  also  remains  to  test  and  implement  the  automatic  picking  of  Sn  and  Lg  phases  at 
ARCESS,  NORESS  and  FINESA. 

During  the  work  with  the  autoregressive  likelihood  estimation  method,  we  have  experi¬ 
enced  that  the  display  of  the  likelihood  functions,  as  illustrated  in  Figs.  7.2.9  and  7.2.10 
can  assist  the  analyst  in  picking  correct  phase  onsets.  In  the  context  of  interactive  analysis 
of  seismic  data,  we  believe  that  the  idea  of  making  such  likelihood  functions  available  to 
the  analyst  should  be  pursued. 


T.  Kvaerna 
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Mine 

Location 

I,  Underground 

67.6702°N 

33.7285°E 

I,  Open-pit 

67.665°N 

33.744°E 

II 

67.647°N 

33.761°E 

III 

67.63 1°N 

33.835°E 

IV 

67.624°N 

33.896°E 

V 

67.632°N 

34.01 1°E 

VI 

67.665°N 

34.146°E 

Table  7.2,1.  The  table  gives  the  location  of  the  six  mines  in  the  Khibiny  Massif  shown  in 
Fig.  7.2.1. 
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Event 

Mine 

ESTONl  P-arrival 
(APAO) 

Size 

l 

I(U) 

1992-312:06.31.49.425 

185  tons 

2 

I(U) 

1992-355:07.10.51.931 

190  tons 

3 

1(0) 

1992-366:04.35.22.825 

14  tons 

4 

1(0) 

1993-017:06.51.15.831 

4  tons 

5 

I(U) 

1993-031:04.14.11.263 

140  tons 

6 

KO) 

1993-038:07.13.16.881 

7  tons 

7 

I(U) 

1993-045:07.39.45.200 

185  tons,  double 

8 

1(0) 

1993-052:07.24.49.706 

14  tons,  double 

9 

1(0) 

1993-059:04.24.01.850 

30  tons 

10 

1(0) 

1993-066:06.58.54.356 

12  tons 

11 

1(0) 

1993-073:11.28.09.731 

10  tons 

12 

1(0) 

1993-080:11.26.36.106 

10  tons 

13 

II 

1992-334:04.08.51.725 

130  tons 

14 

II 

1992-366:09.41.50.156 

100  tons 

15 

II 

1993-024:08.04.48.656 

24  tons 

16 

II 

1993-066:04.14.37.475 

60  tons 

17 

III 

1992-354:08.38.32.806 

15  tons 

18 

III 

1992-361:09.53.07.550 

320  tons 

19 

in 

1993-016:09.20.47.756 

20  tons 

20 

m 

1993-030:10.16.00.581 

20  tons 

21 

hi 

1 993-030: 1 3.07. 1 5.650 

265  terns 

22 

hi 

1993-044:12.05.41.725 

300  tons 

23 

hi 

1993-051:07.28.29.931 

16  tons 

24 

iii 

1993-058:12.01.17.356 

17  tons 

25 

in 

1993-065:08.51.42.356 

18  tons 

26 

hi 

1993-065:11.00.36.931 

1 30  tons 
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Event 

Mine 

ESTON1  P-arrival 
(APAO) 

Size 

27 

II 

1993-065:11.29.46.150 

Induced  earth¬ 
quake,  Ml  2.3 

28 

III 

1993-086:10.35.05.925 

14  tons 

29 

IV 

1992-353:12.29.32.681 

340  tons 

30 

IV 

1992-360:12.37.40.556 

380  tons 

31 

IV 

1992-365:11.39.53.300 

214  tons 

32 

IV 

1993-006:11.01.24.331 

310  teats 

33 

IV 

1993-012:12.04.00.406 

130  tons 

34 

IV 

1993-022:12.50.56.806 

310  tons 

35 

IV 

1993-036:12.35.08.831 

230  tons 

36 

IV 

1993-043:13.37.27.181 

350  tons 

37 

IV 

1993-057:12.31.59.475 

380  tons 

38 

IV 

1993-064:12.32.40.931 

300  tons,  double 

39 

IV 

1993-075:12.29.49.650 

70  tons 

40 

IV 

1993-082:15.08.34.250 

95  tons 

41 

IV 

1993-089:14.44.12.256 

100  tons 

42 

V 

1992-346:08.56.58.906 

183  tons 

43 

V 

1992-360:08.03.37.581 

138  tons 

44 

V 

1993-006:08.18.41.606 

101  tons 

45 

V 

1993-015:14.08.51.075 

188  tons 

46 

V 

1993-029:10.07.17.831 

195  tons 

47 

V 

1993-036:08.05.42.006 

unknown  size 

48 

V 

1993-050:12.34.53.506 

203  tons 

49 

EM 

1993-057:10.33.31.006 

149  tons 

50 

V 

1993-078:09.45.13.906 

260  tons,  general 

51 

V 

1993-085:08.59.22.900 

146  tons 

52 

VI 

1992-318:06.41.16.250 

60  tons 

53 


VI 


1992-332:07.05.26.950 


unknown  size 
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Event 

Mine 

EST0N1  P-arrival 
(APAO) 

Size 

54 

VI 

1992-353:07.35.26.956 

73  tons 

55 

VI 

1992-365:07.23.22.631 

80  tons 

56 

VI 

1993-022:08.20.09.231 

285  tons 

57 

VI 

1993-043:08.13.29.381 

140  tons 

58 

VI 

1993-071:07.33.22.906 

140  tons,  general 

Table  7.2.2.  This  table  contains  information  on  the  58  Khibiny  Massif  events  analyzed  in 
this  study.  For  each  event,  the  assigned  mine,  the  automatic  P-onset  at  the  Apatity 
array  and  the  event  size  are  given.  I(U)  means  the  underground  Mine  I,  whereas  1(0) 
means  the  open-pit  mine.  When  the  event  size  is  given  as  “double”,  the  reported 
charge  is  distributed  between  two  explosions  that  are  closely  separated  in  time.  If  the 
term  “general”  is  used,  the  reported  charge  is  distributed  among  several  explosions 
that  are  closely  separated  in  time. 


Mine 

APAO 

APZ9 

ARCESS 

Delta  (km) 

Az  (deg) 

Delta  (km) 

Az  (deg) 

Delta  (km) 

Az  (deg) 

I(U) 

32.0 

76.8 

17.8 

50.4 

393.8 

118.0 

1(0) 

32.6 

78.0 

18.0 

53.1 

394.7 

118.0 

n 

33.0 

81.7 

17.5 

59.8 

396.4 

118.2 

hi 

35.9 

85.2 

19.6 

69.0 

400.0 

118.1 

IV 

38.4 

86.6 

21.8 

73.3 

402.6 

118.0 

m 

43.4 

85.7 

26.7 

74.4 

406.0 

117.5 

VI 

49.4 

81.8 

33.3 

70.8 

408.6 

116.6 

Table  7.2 J.  Distance  and  azimuths  from  the  three  stations  APAO,  APZ9  and  ARCESS  to 
the  mines  considered  in  this  study. 
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^automatic 

P,  Apatity  array 

0.04  s 

0.04  s 

P,  Apatity  3 -comp. 

0.06  s 

0.05  s 

Pn,  ARCESS 

0.06  s 

0.05  s 

Table  7.2.4.  Estimated  standard  deviations  of  manual  and  automatic  (ESTON1)  P-onsets 
at  the  Apatity  array,  the  Apatity  three-component  station  and  at  ARCESS. 


S-phase 

Reference  phase 

°  manual 

-automatic 

APAO 

Pat  APAO 

0.20 

0.19 

P  at  APZ9 

0.20 

0.19 

Pn  at  ARCESS 

0.18 

0.20 

AMO  • 

IW 

H1B1! 

APZ9 

Pat  APAO 

0.13 

0.15 

P  at  APZ9 

0.12 

0.14 

Pn  at  ARCESS 

0.13 

0.16 

APZ9 

' 

Afr6t*gr  "  ■-  ' 

0.13 

0.15 

Table  7.2.5.  Estimated  standard  deviations  of  manual  and  automatic  (ESTON3)  S-onsets 
at  APAO  and  APZ9.  In  addition  to  the  average  uncertainty,  we  give  for  each  phase 
the  uncertainty  estimates  computed  using  the  three  different  reference  phases. 
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Mean  (°) 


P,  Apatity  array 

8  2 

Pn,  ARCESS 

4.6 

S,  Apatity  array 

-1.0 

Rg,  Apatity  array 

-1.8 

Table  7.2.6.  Mean  and  standard  deviation  of  azimuth  residuals  relative  to  Khibiny  mine 
locations. 
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33*  00'  33*  30'  34*  00' 


Fig.  7.2.1.  In  the  upper  part ,  a  large  reference  area  is  shown.  The  location  of  the 
ARCESS  array  is  given  by  a  filled  circle,  and  the  location  of  the  Khibiny  Massif 
region  is  shown.  The  lower  part  shows  a  detailed  picture  of  the  Khibiny  Massif 
region.  The  locations  of  the  six  mining  sites  are  given  by  laige  numbers  1-6.  The 
Apatity  array  is  shown  as  a  filled  circle  and  the  three-component  station  in  the  town 
of  Apatity  is  shown  as  a  large  triangle. 
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Mine  1.  140  tons,  distance  18  km 
P:  Prewhitened,  S:  Sondooss  2.0  -  g.U  Hz.  Ra:  Sondooss  0.8  -  2.0  Hz 
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12459. 


2646.2 
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4880.7 
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Fig.  7.2.2.  Seismograms  from  an  event  at  the  underground  Mine  I  recorded  at  the  Apatity 
three-component  station.  The  source-receiver  distance  is  18  km.  The  upper  trace, 
emphasizing  the  P-phase,  is  the  z-component  prewhitened  with  a  low-order  predic¬ 
tion  error  filter  designed  from  a  25  s  noise  sample  preceding  the  P-phase.  The  three- 
component  data  rotated  with  the  azimuth  to  the  mine  can  be  seen  at  traces  2-4.  The 
data  are  filtered  in  a  relatively  wide  passband  of  2  -  8  Hz,  and  these  three  traces 
clearly  show  an  instantaneous  S-phase.  The  lower  trace  is  bandpass  filtered  between 
0.8  and  2.0  Hz,  and  clearly  illustrate  the  low-frequency  Rg  phase. 
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Fig.  7.23.  Seismograms  from  an  event  at  mine  VI  recorded  at  the  Apatity  array.  The 
source-receiver  distance  is  49  km.  The  only  difference  from  the  traces  of  Fig.  7.2.2  is 
that  the  upper  trace  is  the  array  beam  steered  with  steering  delays  corresponding  to 
the  slowness  and  azimuth  of  the  P-arrival.  Notice  that  in  this  case  the  S-phase  is 
much  more  emergent  than  the  S-phase  of  Fig.  7.2.2. 
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Fig.  7.2.4.  ARCESS  recording  of  the  event  also  shown  in  Fig.  7.2.2.  The  source-receiver 
distance  is  394  km.  The  upper  trace  is  the  array  beam  steered  with  the  slowness  and 
azimuth  of  the  P-phase.  The  beam  is  prewhitened  with  a  low-order  prediction  error 
filter.  The  lower  three  traces  are  the  rotated  three-component  data  at  the  central  ele¬ 
ment  (ARAO)  of  the  ARCESS  array.  The  data  are  bandpass  filtered  between  2  and  5 
Hz.  Notice  that  both  Pn,  Pg,  Sn  and  Lg  are  clearly  seen,  and  that  Sn  and  Lg  are  most 
prominent  on  the  transverse  component. 
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a) 


C  Apatity  array 


Mean:  0.0590 
St.  Dev:  0.4329 


b) 


SigPro  -  Manual  P-onset  (s) 


Mean:  0.0493 
St.  Dev:  0.0192 


Estonl  -  Manual  P-onset  (s) 

Fig.  7.2J5.  These  two  plots  show  the  time  differences  between  the  automatic  P-onsets  and  • 

the  manual  pick  at  APA0  as  a  function  of  signal-to-noise  ratio  (SNR)  on  the  prewhit¬ 
ened  beam.  All  58  P-observations  are  included.  Figure  a)  shows  the  time  differences 
between  the  automatic  P-onsets  from  the  IMS  and  the  manual  picks,  whereas  figure 
b)  shows  the  time  differences  between  the  automatic  P-onsets  from  ESTON 1  and  the 
manual  picks.  The  two  vertical  dotted  lines  represent  ±o2.  9 
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a) 


b) 


O  ARCESS  array 


Mean:  -0.132 
St.  Dev:  0.1287 

-2-10  1  2 
SigPro  -  Manual  P-onset  (s) 


Mean:  0.0755 
St.  Dev:  0.0317 


-2-10  1  2 


Estonl  -  Manual  P-onset  (s) 


Fig.  7.2.6.  Plots  similar  to  those  of  Fig.  7.2.5,  but  in  this  case  for  Pn-onsets  at  ARCESS. 
For  one  of  the  Khibiny  Massif  events,  the  ARCESS  array  was  not  operational,  and 
consequently  only  57  P-phases  were  analyzed. 
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-0.1  0.0  0.1 


Estonl  -  Manual  P-onset  (s) 

« 


4 

Fig.  7.2.7.  This  figure  show  the  time  differences  between  the  automatic  P-onsets  from 
ESTON1  and  the  manual  pick  at  both  Apatity  stations  (APAO  and  APZ9)  as  a  func¬ 
tion  of  signal-to-noise  ratio  (SNR).  Compared  to  Fig.  7.2.5  the  x-axis  is  strongly 
4  expanded.  Notice  that  P-phases  at  APZ9  (triangles)  generally  have  lower  SNR.  The 

average  dominant  frequency  of  these  P-phases  is  13  Hz. 
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Fig.  12.9.  The  seismograms  of  this  figure  are  Apatity  array  (APAO)  recordings  of  the 
Mine  I  event  also  given  in  Fig.  7.2.2.  The  source-receiver  distance  is  32  km.  Trace 
no.  2  is  the  prewhitened  P-beam  used  for  onset-time  estimation  using  ESTON1.  The 
upper  trace  gives  the  likelihood  function  from  ESTON1  after  processing  an  interval 
of  ±2  s  around  the  initial  P-onset  estimate,  and  the  peak  of  this  likelihood  function 
correspond  to  the  estimated  onset  The  lower  trace  is  the  vertical  component 
APA0_sz  filtered  in  a  low  passband  (0.8  -  2.0  Hz)  to  enhance  the  Rg  phase,  and  the 
trace  above  it  is  the  STA  envelope.  The  peak  of  this  envelope  is  defined  as  the  peak 
of  the  Rg  phase.  After  the  P-onset  and  the  Rg  maximum  is  found,  we  are  searching 
for  the  S-onset  using  the  three-component  ESTON3  estimator.  The  search  interval 
starts  2  s  after  the  P-onset  and  stops  at  the  Rg  peak,  as  seen  from  the  ESTON3  likeli¬ 
hood  function  of  trace  no.  3.  The  three-component  data  processed  by  ESTON3  are 
given  in  traces  4-7.  Notice  that  there  is  no  need  to  rotate  the  three-component  data 
before  using  ESTON3,  but  to  visualize  that  the  S-phase  has  the  laigest  SNR  on  the 
transverse  component,  we  have  in  this  figure  rotated  the  data.  Notice  that  both  the 
ESTON1  and  the  ESTON3  likelihood  functions  show  clear  peaks  at  the  P-  and  S- 
onsets. 
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Ffc.  7.2.10  The  traces  of  this  plot  are  similar  to  those  of  Fig.  7.2.9.  The  seismograms  are 
Apatity  array  (APAO)  recordings  of  the  Mine  VI  event  also  shown  in  Fig.  7.2.3.  The 
source-receiver  distance  is  49  km.  Notice  that  the  S-phase  is  much  more  emeigent 
than  the  S-phase  of  Fig.  7.2.9,  but  the  ESTON3  likelihood  function  still  shows  a 
clear  peak  at  the  S-onset. 
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13  Intelligent  post- processing  of  seismic  events  -  Part  3: 

•  Precise  relocation  of  events  in  a  known  target  region 

Introduction 

From  experience  with  analyst  review  of  events  automatically  defined  by  the  Intelligent 

•  Monitoring  System  (IMS)  (Bache  et  al,  1993),  we  have  realized  that  the  quality  of  the 
automatic  event  locations  can  be  significantly  improved  if  the  event  intervals  are  repro¬ 
cessed  with  signal  processing  parameters  tuned  to  phases  from  events  in  the  given  region. 
The  tuned  processing  parameters  are  obtained  from  off-line  analysis  of  events  located  in 
the  region  of  interest.  The  primary  goal  of  such  intelligent  post-processing  is  to  provide 

•  event  definitions  of  a  quality  that  minimizes  the  need  for  subsequent  analysis  and  review. 

In  this  third  paper  on  intelligent  post-processing  of  seismic  events  we  test  an  event  loca¬ 
tion  method  that  makes  use  of  accurate  retiming  of  seismic  phases  and  accurate  re-esti- 
mates  of  azimuth  using  broad-band  analysis  of  array  data.  The  key  to  obtaining  such 
9  accurate  estimates  is  the  knowledge,  provided  by  the  IMS,  of  an  approximate  event  loca¬ 

tion  in  a  known  area  with  good  calibration  information.  Various  aspects  of  generalizing 
the  method  and  incorporating  it  into  the  IMS  are  also  discussed. 

Location  of  regional  events 

®  Seismic  event  location  has  traditionally  relied  upon  P-wave  arrival  times  from  a  large 

number  of  stations.  For  small  events  observed  only  at  regional  distances,  it  has  been  nec¬ 
essary  to  include  arrival  times  of  later  phases  as  well  as  arrival  azimuths  in  order  to  obtain 
acceptable  location  accuracy  (Bratt  and  Bache,  1988;  Thurber  et  al,  1989).  With  a  sparse 
network,  the  number  of  arrival  time  observations  may  often  be  so  low  that  azimuth  infor- 

•  mation  is  given  considerable  weight  in  the  location  scheme.  However,  as  the  number  of 
observations  increases,  the  relative  importance  of  the  azimuths  decreases  sharply,  in  view 
of  the  generally  large  uncertainties  (10  degrees  or  more)  that  are  associated  with  azimuth 
estimates  (e.g.,  Suteau-Henson,  1990;  Bame  et  al,  1990). 

•  Kvaema  and  Ringdal  (1986)  demonstrated  that  the  uncertainty  of  azimuth  estimates  from 
regional  arrays  can  be  greatly  reduced  by  applying  a  broad-band  estimation  scheme  in 
which  the  frequency  band  is  kept  constant  for  a  given  epicentral  region.  Thus,  for  NOR- 
ESS  observations  from  B15sj0  explosions  (distance  300  km)  with  high  SNR,  they  found  an 
azimuth  standard  deviation  of  only  1.0  degrees  for  Pn  phases  and  1 .5  degrees  for  Sn  and 

•  Lg  phases. 

This  observation  suggests  that  given  an  initial  (approximate)  event  epicenter,  it  should  be 
feasible  to  determine  an  optimum  frequency  band  (based  upon  previous  observations)  and 
apply  broad-band  F-k  analysis  in  order  to  obtain  very  accurate  azimuth  estimates,  which 
0  can  then  be  used  to  refine  the  event  location.  The  basic  idea  of  this  paper  is  to  establish  and 

test  a  method  for  location  refinement  which  takes  into  account  both  the  stable  azimuth 
estimates  obtained  by  fixed  frequency  broad-band  analysis  and  the  accurate  timing  obser¬ 
vations  described  in  Section  7.2.  We  use  as  an  example  IMS -reported  events  in  the  Khib- 
iny  Massif  recorded  by  the  nearby  Apatity  array  and  the  more  distant  ARCESS  array.  We 
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show  that  a  significant  improvement  can  be  obtained  relative  to  the  current  IMS  process¬ 
ing.  Moreover,  the  method  is  completely  automatic  and  can  easily  be  incorporated  into  the 
IMS,  as  illustrated  in  Fig.  7.3.1. 

As  shown  in  Section  7.1,  the  majority  of  seismic  events  in  Fennoscandia  can  be  associated 
with  mining  activity  at  a  few  sites.  The  Khibiny  Massif  is  one  of  the  mining  areas  with  the 
most  events.  Obtaining  accurate  automatic  solutions  for  events  from  this  and  other  mining 
areas  will  lead  to  a  considerably  lower  workload  for  the  analyst,  and  will  consequently 
allow  more  time  to  be  spent  analyzing  events  of  particular  interest. 

Method 

The  automatic  method  applied  in  this  study  consists  of  the  following  steps: 

1 .  The  IMS  detects  and  locates  a  seismic  event  within  a  given  distance  from  the  target 
region. 

2.  For  each  phase  detection  by  a  station  in  the  network,  the  estimated  arrival  time  is 
refined  using  the  method  of  Pisarenko  et  al  (1987)  (see  Section  7.2).  A  standard 
deviation  is  assigned  to  each  refined  estimate. 

3.  For  each  phase  detection  by  one  of  the  arrays,  an  "optimum"  frequency  band  for  F-k 
analysis  is  extracted  from  a  data  base  which  we  assume  has  been  previously  estab¬ 
lished  and  calibrated. 

4.  Broad-band  F-k  analysis  is  then  applied  to  each  phase.  The  resulting  values  are  cor¬ 
rected  for  systematic  bias  and  assigned  a  standard  deviation,  based  upon  SNR,  phase 
type  and  previously  observed  calibration  information. 

5.  The  LocSat  program  (Bratt  and  Bache,  1988)  is  then  applied  to  the  revised  data  set, 
and  a  new  event  location  estimate  is  obtained.  In  practice,  zero  depth  is  assumed  a 
priori  in  the  automatic  process. 

In  the  case  study  discussed  in  detail  in  this  paper,  we  will  analyze  events  in  the  Khibiny 
Massif  recorded  by  IMS,  for  which  we  have  independent  location  information.  We  will 
primarily  make  use  of  the  Apatity  array,  using  arrival  times  for  the  P  and  S  phases  and  azi¬ 
muths  derived  from  the  Rg  phase  in  order  to  refine  the  location  estimates  as  described 
above.  In  addition,  we  will  consider  possible  improvement  in  the  location  accuracy  when 
using  additional  available  stations. 

Data 

The  data  for  this  study  comprises  58  events  at  6  mines  in  the  Khibiny  Massif.  The  mines 
are  only  a  few  kilometers  apart  (Mykkeltveit,  1992),  and  they  each  have  a  dimension  rang¬ 
ing  from  a  few  hundred  meters  to  about  1  km.  These  dimensions  are  small  enough  so  that 
we  ignore  the  areal  extent  of  each  mine.  A  list  of  the  event  parameters  is  given  in  Section 
7.2. 

To  obtain  an  initial  location,  we  have  used  the  automatic  IMS  epicenter  solution  for  each 
event.  The  automatic  IMS  solution  relative  to  the  "true"  epicenter  for  each  of  the  events  is 
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displayed  graphically  in  Fig.  7.3.2.  Although  many  of  the  events  are  very  accurately 
located  (to  within  a  few  km),  there  is  a  fair  amount  of  scatter  in  the  location  estimates.  The 
median  "error"  is  10.6  km.  In  a  few  cases,  the  location  is  wrong  by  several  tens  of  km;  this 
is  Jue  to  occasional  erroneous  phase  identification  by  die  automatic  process. 

The  Apatity  array  has  been  described  by  Mykkeltveit  et  al  ( 1992),  and  initial  data  analysis 
results  have  been  presented  by  Ringdal  and  Fyen  (1992).  The  latter  paper  discusses  in  par¬ 
ticular  the  high  stability  of  azimuth  estimates  using  the  Rg  phase  for  shallow  events  at 
local  distances. 

In  Section  7.2,  some  examples  of  Apatity  array  recordings  of  events  in  the  Khibiny  Massif 
are  shown.  The  filtered  recordings  (0.8-2.0  Hz)  are  dominated  by  the  Rg  phase.  In  con¬ 
trast,  the  P  and  S  phases  are  far  more  high-frequent,  as  is  to  be  expected  at  such  close  dis¬ 
tances.  Because  of  the  low  frequency  of  the  Rg  phase,  this  phase  has  a  high  coherency 
across  the  array  and  thus  provides  more  stable  slowness  estimates  than  are  obtained  from 
the  other  phases. 

For  the  present  study,  we  have  used  Apatity  array  recordings  for  all  58  events.  For  all  but 
six  events,  high-frequency  three-component  data  from  the  array  site  have  been  available. 
Furthermore,  all  events  have  associated  3-component  broadband  recordings  from  the  sta¬ 
tion  installed  in  the  town  of  Apatity,  15  km  from  the  array  site.  ARCESS  array  data  have 
been  available  for  all  but  one  of  the  events. 

Results 

In  the  following,  we  present  results  from  analyzing  the  58  events  in  the  data  base  under 
various  scenarios.  Each  scenario  assumes  that  we  have  an  operational  situation  in  which  a 
certain  amount  of  calibration  information  has  been  assembled. 

a)  Use  of  Apatity  array  only;  P  and  S  arrival  times,  Rg  azimuth,  no  calibration. 

In  this  scenario,  we  assume  no  prior  knowledge  except  for  the  "optimum"  frequency 
bands  for  timing  and  azimuth  estimation.  The  travel-time  tables  are  those  that  the 
IMS  uses  for  Fennoscandia  in  general,  with  no  regional  corrections. 

The  results  are  displayed  in  Figs.  7.3.3.  We  note  that  there  is  a  small  systematic  bias, 
as  there  seems  to  be  a  trend  of  events  being  located  to  the  north  and  east  of  the  sites. 
This  indicates  that  some  improvement  could  be  obtained  by  introducing  a  regional 
velocity  model  for  the  Khibiny  area,  together  with  regional  azimuthal  corrections. 

Nevertheless,  the  scatter  relative  to  the  6  mines  in  this  plot  (using  uncalibrated  Apa¬ 
tity  array  data  only)  is  significantly  reduced  compared  to  the  IMS  results  plotted  in 
Fig.  7.3.2. 

The  improvement  can  be  quantified  by  considering  that  the  median  "error"  is  only 
3.1  km,  with  a  worst  case  "error"  of  6.7  km.  This  can  be  compared  to  the  IMS  results 
of  10.6  km  and  75.3  km,  respectively.  Thus  the  reprocessing,  even  without  regional 
correction,  could  significantly  improve  the  precision  in  the  epicenter  solutions  if 
incorporated  in  the  automatic  IMS  processing. 
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b)  Use  of  Apatity  array  data  only,  with  a  regional  P-S  travel-time  bias  correction  and  a 
correction  for  systematic  Rg  azimuth  bias. 

The  resulting  location  plot  is  shown  in  Fig.  7.3.4.  The  median  error  is  reduced  to  2.1 
km,  and  the  "worst  case"  error  is  6.1  km. 

c)  Use  of  Apatity  array  data  plus  P  and  S  times  from  the  3-component  station  in  Apat¬ 
ity.  Regional  corrections  have  been  applied  both  for  the  array  and  the  3-comp,  sta¬ 
tion. 

The  resulting  location  plot  (Fig.  7.3.3)  shows  that  the  median  error  is  1.4  km,  and 
the  maximum  error  is  5.4  km. 

d)  This  final  scenario  includes  both  Apatity  array  data,  Apatity  3-comp,  data  and 
ARCESS  P  arrival  times.  Regional  corrections  have  been  applied  for  all  the  data. 

The  resulting  location  plot  (Fig.  7.3.6)  shows  a  median  error  of  1.6  km,  and  a  maxi¬ 
mum  error  of  7.1  km.  It  is  interesting  to  note  that  this  is  not  quite  as  good  as  case  c). 
Thus,  it  is  not  necessarily  an  improvement  to  add  data,  even  if  (as  in  this  case)  the 
additional  data  are  extremely  accurate  (error  £  +  0.05  seconds  on  P-recordings). 

We  should  note  that  the  travel-time  corrections  used  in  this  paper  are  relative  and  based 
upon  the  same  data  set  that  we  have  evaluated.  In  view  of  the  large  number  of  events,  the 
possible  bias  introduced  by  this  procedure  should  be  negligible.  As  the  origin  time  of  the 
Khibiny  events  are  unknown,  we  have  fixed  the  travel-time  correction  of  P  at  the  Apatity 
three-component  station  APZ9  to  0.  The  P  phase  at  APZO  has  travelled  a  shorter  distance 
and  spent  less  time  in  the  earth  than  the  other  phases  used  in  the  event  locations.  Thus,  the 
influence  of  an  erroneous  travel-time  model  is  likely  to  be  the  least  for  this  phase. 

Travel-time  corrections: 


p 

APZ9 

0  s  (fixed) 

p 

APA0 

-0.10  s 

Pn 

ARCESS 

-0.31  s 

S 

APZ9 

0.22  s 

S 

APA0 

0.09  s 

Conclusions 

In  this  study  we  have  conducted  extensive  off-line  analysis  of  58  mining  explosions  in  the 
Khibiny  Massif  recorded  at  the  Apatity  array.  Independent  locations  of  these  explosions 
are  provided  by  seismologists  from  the  Kola  Science  Centre.  Most  of  these  events  show 
clear  P,  S  and  Rg  phases  at  the  nearby  Apatity  array  located  30-50  km  away  from  the  min¬ 
ing  areas  and  the  events  have  also  been  detected  by  the  ARCESS  array.  By  using  the 
onset-time  estimator  of  Pisarenko  et  al  (1987)  and  comparing  to  the  best  manual  pick,  we 
have  found  that  P-onsets  can  be  automatically  estimated  with  an  accuracy  of  better  than 
0.1  seconds,  and  S-onsets  with  an  accuracy  better  than  0.5  seconds.  In  addition,  the  azi¬ 
muth  estimates  from  F-k  analysis  of  the  P  and  Rg  phases  show  to  be  accurate  well  within 
±  5  degrees  after  removing  the  biases.  The  key  to  achieving  stable  azimuth  estimates  is  to 
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process  the  data  in  a  fixed  frequency  band,  using  a  fixed  time  window  positioning,  as  dem¬ 
onstrated  by  Kvsema  and  Ringdal  (1986)  for  BlSsja  recordings  at  NORESS. 

Our  observations  suggest  that  based  on  data  from  the  Apatity  array  alone,  we  are  able  to 
locate  these  events  (assuming  0  depth)  with  a  median  error  of  about  3  km  relative  to  the 
true  location.  Even  better  accuracy  can  be  achieved  using  calibration  information,  i.e.,  cor¬ 
recting  the  azimuth  and  arrival  time  observations  for  systematic  bias.  The  excellent  preci¬ 
sion  of  the  automatic  phase  onsets  and  azimuth  estimates  also  indicate  that  the  need  for 
subsequent  manual  analysis  of  these  events  may  be  eliminated. 

After  we  have  established  the  tuned  processing  parameters  for  events  in  this  region,  a  nat¬ 
ural  next  step  will  be  to  automatically  activate  such  intelligent  post-processing  every  time 
the  IMS  locates  an  event  in  the  Khibiny  region.  Extensions  of  the  method  to  other  mining 
areas  will  also  be  considered. 

T.  Kvseraa 
F.  Ringdal 
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Fig.  7 .3.1.  Schematic  view  of  the  principle  behind  intelligent  post-processing  of  seismic 
events. 
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Expert  system  location  errors 
Median:  10.6  km  Maximum:  85.3  km 
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Fig.  73.2.  Location  of  the  six  mining  sites  in  the  Khibiny  Massif  (large  numbers  1-6)  and 
the  locations  of  the  58  reference  events  (small  numbers  1-6)  as  given  by  the  auto¬ 
matic  IMS  processing.  In  the  upper  part,  a  large  reference  area  is  shown,  with  the 
mines  plotted  as  filled  squares.  The  lower  part  shows  a  detailed  picture  for  the  area 
near  the  mines.  The  small  number  (1-6)  associated  with  each  event  represents  the 
mine  in  which  die  event  actually  occurred.  The  Apatity  array  is  shown  as  a  filled  cir¬ 
cle  and  the  three-comp,  station  in  the  town  of  Apatity  is  shown  as  a  filled  triangle. 
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Apatity  array  location  errors  (uncalibrated) 
Median:  3.1  km  90%  quantile:  4.9  km  Maximum:  6.7 ! 


31*  32“  33“  34“  35“  31 


Fig.  7.3.3.  Same  as  Fig.  7.3.2,  but  the  event  locations  (small  numbers)  hav 
taken  from  the  post-processing  results  using  uncalibrated  Apatity  am 
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Apatity  stations  location  errors  (calibrated) 
Median:  1.4  km  90%  quantile:  2.7  km  Maximum:  5.4  km 
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Fig.  73.5.  Same  as  Fig.  7.3.2,  but  with  the  event  locations  resulting  from  combining  Apa 
tity  array  data  and  3-comp,  data  (both  calibrated)  in  the  post-processing. 
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7.4  Monitoring  a  moratorium:  An  experiment  in  continuous  seismic 
threshold  monitoring  of  the  northern  Novaya  Zemlya  test  site 

Introduction 

After  1  October  1992  a  moratorium  on  nuclear  testing  has  been  in  effect  for  the  US,  UK, 
Russia  and  France.  On  this  background,  we  have  applied  the  continuous  threshold  moni¬ 
toring  technique  to  the  northern  Novaya  Zemlya  test  site  for  a  full  six-month  period  (1 
October  -  31  March  1993),  using  the  NORESS,  ARCESS  and  FINESA  regional  arrays. 
Starting  1  December  1992,  we  have  compiled  daily  statistics  of  all  peaks  on  the  threshold 
diagram  exceeding  m b  =  2.73,  and  associated  these  peaks  to  regional  or  teleseismic  events 
whenever  possible.  In  addition,  we  have  analyzed  smaller  peaks  (below  -  2.75)  that 
can  possibly  be  associated  with  Novaya  Zemlya  epicenters. 

The  theoretical  background  for  and  applications  of  the  continuous  seismic  threshold  mon¬ 
itoring  metod  (CSTM)  have  been  described  in  several  articles.  The  approach  was  intro¬ 
duced  by  Ringdal  and  Kvaema  (1989),  who  showed  that  by  continuously  monitoring  the 
seismic  amplitude  level  at  several  seismic  stations  or  arrays,  one  can  at  any  time  obtain  an 
instant  network- based  magnitude  threshold  for  a  given  target  region.  The  magnitude 
threshold  can  be  interpreted  as  the  maximum  magnitude  of  a  possible  clandestine  explo¬ 
sion,  given  a  predefined  level  of  confidence.  In  the  context  of  a  comprehensive  or  thresh¬ 
old  test  ban  treaty,  the  continuous  assessment  of  the  magnitude  thresholds  makes  it 
possible  to  focus  attention  upon  those  specific  time  intervals  when  realistic  evasion  oppor¬ 
tunities  exist,  while  retaining  confidence  that  no  treaty  violation  has  occurred  at  other 
times. 

Previous  results  from  experimental  monitoring 

Kvaema  (1992)  presented  results  from  a  one-month  experiment  of  continuously  monitor¬ 
ing  the  northern  Novaya  Zemlya  test  site.  Data  from  the  Fennoscandian  regional  array  net¬ 
work  (ARCESS,  FINESA,  and  NORESS),  see  Fig.  7.4.1,  were  used  to  calculate  the 
magnitude  thresholds.  It  was  found  that  the  test  site  could  be  consistently  monitored  at  a 
very  low  magnitude  level  (typically  mb  =  2.5).  In  fact,  every  occurrence  of  the  threshold 
exceeding  mb  =  2.6  could  be  explained  as  resulting  from  an  identified  interfering  event 
signal  either  at  teleseismic  or  regional  distance,  except  for  three  instances  when  a  short 
gap  in  ARCESS  recording  caused  the  network  threshold  to  increase. 

The  excellent  capability  of  the  Fennoscandian  regional  array  network  to  monitor  the 
northern  Novaya  Zemlya  test  site  was  further  confirmed  by  an  experiment  where  record¬ 
ings  of  the  Novaya  Zemlya  nuclear  test  of  October  24, 1990  were  downscaled  to  mb  =  2.6 
and  superimposed  on  different  noise  intervals  (Kvaema,  1991). 

In  the  context  of  using  CSTM  as  a  tool  in  routine  monitoring,  it  is  important  to  determine 
how  the  method  will  work  under  different  conditions.  Variability  in  the  seismic  noise 
level,  occurrences  of  large  earthquakes  and  aftershock  sequences,  station  downtimes  and 
data  quality  problems  are  all  factors  that  will  influence  the  performance  of  CSTM. 


NORSAR  Sci.  Rep.  2-92/93 


May  1993 


Analysis  of  network  threshold  peaks 

Our  monitoring  experiment  was  conducted  in  the  same  way  and  with  the  same  parameter 
settings  as  used  by  Kvaema  (1992).  In  Kvaema  (1992)  the  monitoring  results  were  pre¬ 
sented  in  terms  of  plots  covering  one  data  day  each.  In  Figs.  A-l  to  A-29  of  the  Appendix 
of  that  report,  each  covering  one  day  of  February,  1992,  all  time  periods  where  the  net¬ 
work  magnitude  thresholds  at  the  90%  confidence  level  exceeded  m*,  =  2.6  were  identi¬ 
fied.  A  similar  approach  was  used  for  this  experiment 

For  the  remainder  of  this  paper,  the  term  magnitude  threshold  implies  the  magnitude 
threshold  at  the  90%  confidence  level. 

Figs.  7.4.2  shows  a  typical  example  of  a  one-day  plot  The  upper  three  traces  of  each  fig¬ 
ure  represent  the  magnitude  thresholds  obtained  from  the  three  indivitual  arrays,  whereas 
the  bottom  trace  illustrates  the  network  threshold.  Typically,  the  individual  array  traces 
have  a  number  of  significant  peaks  for  each  24-hour  period,  due  to  signals  from  interfering 
events  (regional  or  teleseismic).  On  the  network  trace,  the  number  and  sizes  of  these  peaks 
are  significantly  reduced,  because  an  interfering  event  usually  will  not  provide  matching 
signals  at  all  stations.  From  probabilistic  considerations,  it  can  in  such  cases  be  inferred 
that  the  actual  network  threshold  is  lower  than  these  individual  peaks  might  indicate. 

The  arrows  on  the  one-day  threshold  plots  indicate  peaks  with  network  magnitude  thresh¬ 
old  exceeding  mb  =  2.75.  A  T  at  the  arrow  indicates  that  the  peak  is  caused  by  signals 
from  a  teleseismic  event,  whereas  an  R  indicates  signals  from  a  regional  or  local  event. 
Fig.  7.4.3  shows  summary  statistics  for  one  data  day,  with  an  explanation  of  all  threshold 
peaks  exceeding  2.75.  Such  summary  statistics  were  generated  for  each  day  during  the 
four-month  period,  1  December  1992  -  31  March  1993. 

The  peaks  in  the  threshold  traces  are  caused  by  either  large  teleseismic  events  or  by 
regional  events.  The  regional  events  were  all  processed  and  located  by  the  Intelligent 
Monitoring  System  (IMS)  (Bache  et  al.,  1993),  and  die  teleseismic  events  were  located 
either  by  the  IMS  or  by  the  QED  service  of  the  USGS. 

Fig.  7.4.4  shows  a  histogram  of  the  number  of  peaks  exceeding  given  magnitude  thresh¬ 
olds.  During  the  entire  four-month  period,  there  were  only  40  peaks  exceeding  mb  =  3.0. 
Each  of  these  peaks  could  be  unambiguously  associated  with  either  a  regional  or  teleseis¬ 
mic  event  Consequently,  at  the  specified  confidence  level,  we  can  state  that  no  seismic 
event  of  mb  >  3.0  occurred  at  the  test  site  during  this  four-month  period. 

The  event  at  Novaya  Zemlya  on  31  December  1992 

Figs.  7.4.5  and  7.4.6  show  the  threshold  plots  and  peak  statistics  for  31  December  1992. 
At  0929  GMT  that  day,  a  peak  occurred  corresponding  to  an  event  at  Novaya  Zemlya, 
located  by  the  IMS.  This  peak  (which  is  the  only  one  associated  to  a  Novaya  Zemlya  event 
during  the  six-month  interval)  had  an  upper  magnitude  limit  of  2.6.  The  actual  event  mag¬ 
nitude,  according  to  Carter  et  al  (1993)  was  2.5. 
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This  event  was  detected  by  ARCESS  (P  and  S),  Spitsbergen  (P  and  S),  NORESS  (P)  and 
Apatity  (S).  In  addition,  the  Kola  Science  Centre  provided  readings  for  the  station 
Amderma  (Pg  and  Sn)  made  from  analog  recordings.  Tables  7.4.1  and  7.4.2  summarize 
the  available  observations  for  this  event  Table  7.4.3  gives  the  results  of  applying  the 
LOCSAT  program  (Bratt  and  Bache,  1988)  to  this  parameter  set. 

Our  results  indicate  that  the  epicenter  was  slightly  to  the  north  of  the  test  site.  However, 
there  are  some  uncertainties  in  the  travel  time  tables,  and  an  on-site  location  cannot 
entirely  be  ruled  out.  We  are  hesitant  to  introduce  travel-time  corrections  in  this  case,  since 
no  good  reference  event  is  available  for  some  of  the  key  stations. 

A  plot  of  the  IMS-processed  traces  is  shown  in  Fig.  7.4.7.  Notice  in  particular  the  high 
SNR  at  the  Spitsbergen  B2  single  sensor  (filter  band  8-16  Hz)  and  on  the  ARCESS  array 
beam.  Fig.  7.4.8  shows  the  IMS  solution  for  this  event  The  Amderma  station  is  not 
included  in  the  IMS  processing,  but  adding  its  data  does  not  cause  any  significant  change 
in  the  event  location. 

Conclusions 

This  work  has  documented  the  practical  capability  of  the  Continuous  Seismic  Threshold 
Monitoring  method  to  monitor  a  specific  nuclear  test  site  at  a  very  low  threshold  over  an 
extended  time  period. 

Specifically,  we  have  used  the  Fennoscandian  array  network  (NORESS,  ARCESS  and 
FINESA)  to  monitor  the  northern  Novaya  Zemlya  test  site  for  a  full  four-month  period  at  a 
threshold  of  mb  =  2.75.  We  have  identified  only  one  instance  where  an  event  close  to  this 
threshold  has  occurred  near  the  test  site.  In  fact,  the  event  magnitude  (2.5)  was  below  our 
target  threshold,  but  the  peak  was  still  easily  identified  on  the  threshold  trace. 

Recently,  additional  array  stations  have  been  installed  or  are  planned  for  installation  in  the 
Arctic  region.  These  stations  would  contribute  to  further  improving  the  CSTM  capability, 
both  for  Novaya  Zemlya  and  on  a  general  regional  basis.  This  will  be  the  subject  for  addi¬ 
tional  studies  in  the  future. 

T.  Kvaerna 
F.  Ringdal 
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Station 

Lat  (°N) 

Lon  (°E) 

Type  of  Station 

AMD  (Amderma) 

69.742 

61.655 

3-comp  analog 

APA  (Apatity) 

67.600 

33.000 

Array 

ARAO  (ARCESS) 

69.535 

25.506 

Array 

FIAO  (FINESA) 

61.444 

26.079 

Array 

NRAO  (NORESS) 

60.735 

11.541 

Array 

SVA  (Spitsbergen) 

78.180 

16.350 

Array 

Table  7.4.1.  Location  of  stations  used  in  this  study. 


Station 

Phase 

Type1* 

Arrival  time 

Azimuth 

St  dev.2) 

AMD 

Pg 

t 

09.30.43.7 

- 

2.0 

AMD 

Sn 

t 

09.31.20.7 

- 

2.0 

APA 

Sn 

t 

09.33.22.0 

- 

2.0 

ARC 

Pn 

t 

09.31.48.7 

- 

1.0 

ARC 

Pn 

a 

- 

62.5 

15.0 

ARC 

Sn 

t 

09.33.37.2 

- 

2.0 

ARC 

Sn 

a 

- 

58.4 

15.0 

NRS 

P 

t 

09.34.04.3 

- 

1.0 

NRS 

P 

a 

- 

24.0 

15.0 

SVA 

Pn 

t 

09.31.50.7 

- 

1.0 

SVA 

Sn 

t 

09.33.41.7 

- 

2.0 

Table  7A2.  Observed  arrivals  for  the  31  Dec  92  event 
’’  t  =  time,  a  *  azimuth 

2)  A  priori  standard  deviation  in  seconds  (time)  or  degrees  (azimuth) 
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Final  location  estimate  (+/-  S.D.): 


Utituda : 
Longitude 
Depth: 
Relative  O.T. : 
Absolute  O.T. : 


73.620  dog. 
SS.196  d eg. 
0.000  Ina. 
-79.301  sac. 
-79.301  sac. 


6.673  ka. 

7.177  Ion. 

0.000  ka.  (Fiaad) 
0.849  sac. 

0.849  sac. 


:  1969  12  31  23:58:40.70 


Confidence  region  at  0.90  level: 
Baal  an jor  axis:  18.1  ka. 

Saati -minor  axis:  10.6  ka. 

Major-axis  strika:  49.3  dag. 

Orig.  tiaa  arror:  1.4  sac. 


18.1  ka.  ■  0.16  dag. 

10.6  ka.  -  0.10  dag. 

49.3  dag.  clockwise  froa  North 
1.4  sac. 


Standard  arrors  (signs) : 

Prior:  1.00  (  9999  dag.  of  fraadoa) 
Posterior:  1.00  (  10007  dag.  of  fraadoa) 
Posterior:  0.65  (Normalised  sample  S.D.) 


Azimuthal  weighting:  1.00 
Effective  rank  of  matrix:  2.00 
Ms a Inna  azimuthal  GAP:  195  dag. 
-  No  daag>ing  required  ! 


Data 

Residuals 

Distance 

Azimuth 

Data 

Ariv  ID  Statn  Phase  Type  at 

True  No  realized 

(deg.) 

(deg.) 

Import  Err 

245771 

APA 

Sn 

245771 

AW 

*9 

245771 

AW 

Sn 

245782 

MRS 

P 

245782 

MRS 

P 

245782 

ARC 

Pn 

245782 

ARC 

Pn 

245782 

ARC 

Sn 

245782 

ARC 

Sn 

245782 

SVA 

Pn 

245782 

SVA 

Sn 

-0.087 

-0.044 

9.417 

241.62 

0.314 

0.900 

0.450 

4.371 

149.27 

0.319 

-1.495 

-0.748 

4.371 

149.27 

0.455 

1.068 

1.068 

20.495 

254.54 

0.016 

-9.782 

-0.652 

20.495 

254.54 

0.000 

-0.317 

-0.317 

10.096 

261.08 

0.055 

9.672 

0.645 

10.096 

261.08 

0.000 

-1.009 

-0.504 

10.096 

261.08 

0.220 

5.572 

0.371 

10.096 

261.08 

0.000 

-0.320 

-0.320 

10.244 

313.75 

0.408 

-0.033 

-0.016 

10.244 

313.75 

0.211 

Table  7.4  J.  Epicenter  solution  for  the  31  Dec  92  event  at  Novaya  Zemlya  using  the  LOC- 
SAT  program  (Bratt  and  Bache,  1988).  The  depth  has  been  constrained  to  0.  See 
Table  7.4.1  for  station  locations  and  Table  7.4.2  for  observed  arrival  data  and 
assumed  a  priori  standard  deviations. 
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Threshold  monitoring  of  Novaya  Zemlya 


Fig.  7.4.1.  Location  of  the  target  area  (the  northern  Novaya  Zemlya  test  site)  for  the  mon¬ 
itoring  experiment  The  locations  of  the  three  arrays  NORESS  (A  =  2280  km), 
ARCESS  (A  =  1100  km)  and  FINES  A  (A  =  1780  km)  are  indicated. 
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Network  maximum:  2.99 


u 

Network  m 

is 


H 

» 


T 

* 


.Lu 


10  tt  12  13  14  IS  IS  17  IS  It  30  21  22  23  24 


Hours  (GMT) 
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Fig.  74.2.  Example  of  continuous  threshold  monitoring  of  Novaya  Zemlya  for  one  day 
(13  February  1993).  A  “threshold  trace”  is  shown  for  each  of  the  3  arrays,  and  the 
combined  network  threshold  trace  is  shown  at  the  bottom.  Note  that  the  network 
threshold  is  well  below  magnitude  2.3  almost  all  the  time. 
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Histogram  of  TM  peaks 


Magnitude  230-2.75  2.75-3.00  3.00  -  3.25  3.25-330  330  -> 

No  of  peaks  554  67  24  10  6 


Fig.  7.4.4.  Histogram  showing  the  distribution  of  peaks  on  the  network  threshold  trace  for 
the  four-month  monitoring  experiment 
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Network  maximum:  3.89 


Fig.  7.4.5.  Threshold  plot  for  31  Dec  92  (see  Fig.  7.4.3  for  explanation).  The  small 
Novaya  Zemlya  event  at  09.29  GMT  is  indicated. 
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Waveforms  of  Novaya  Zemiya  Event  31  Dec  1992  •  Magnitude  2.5 


Apatity 

ARCESS 

Svalbard 

FINESA 

NORESS 


Fig.  7.4.7.  Waveform  plots  for  5  regional  arrays  for  the  Novaya  Zemiya  event  on  31  Dec  • 

92.  There  are  P-phase  detections  at  ARCESS,  Svalbard  (Spitsbergen)  and  NORESS, 
and  S-phase  detections  at  Apatity,  ARCESS  and  Svalbard. 
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Novaya  Zemlya  event  31  Dec  1992 


OritfnM: 

saaat 

Thee : 

31  Ok  IM  0333*4.4.- 12  v 

Latfhede : 

73J 

Longitude : 

SS2 

Depth : 

oo  ♦-0.0 

Body  wave  Magnitude : 

314 

Other  Magnitude  : 

227 

Locating  detections  : 

« 

Semi -Major  mds  : 

47J 

Semi-Minor  axis  : 

202 

Strike  of  semi-major  axis  : 

1112 

Standard  error  of  observation  : 

10000 

Estimated  type  : 

r 

Source/orgtnator : 

Locatof:bemt 

Fig.  7.4.8.  IMS  solution  and  associated  error  ellipse  for  the  Novaya  Zemlya  event  (m^  ~ 
2.5)  on  31  Dec  92. 
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7i  Initial  processing  results  from  the  Spitsbergen  small-aperture 
array 

During  the  fall  of  1992,  a  new  regional  array  was  installed  near  the  town  of  Longyearbyen 
on  the  Arctic  island  of  Spitsbergen  (see  Mykkeltveit  et  al,  1992).  Fig.  7.5.1  shows  the 
location  of  this  array  and  Fig.  7.5.2  shows  the  array  geometry.  As  of  early  1993,  six  of  the 
nine  array  sites  had  been  finalized,  and  this  study  is  based  on  these  six  instruments. 

In  the  following  we  report  on  some  initial  results  from  analyzing  data  from  the  new  array. 
It  is  emphasized  that  these  results  are  preliminary,  since  the  available  data  cover  only  a 
few  months,  and  only  a  partial  array  has  been  installed  until  now.  Thus  a  more  comprehen¬ 
sive  assessment  must  await  the  collection  of  more  complete  array  data  over  a  longer  time 
span. 

Noise  spectra 

Fig.  7.5.3a  shows  an  example  of  corrected  noise  spectra  for  the  6  vertical  elements  of  the 
Spitsbergen  array,  taken  at  00.00  hours  (GMT)  on  day  345, 1992.  For  comparison,  the 
ARCESS  spectra  for  the  same  number  of  channels  and  taken  at  the  same  time  are  shown 
in  Fig.  7.5.3b.  From  these  figures,  it  is  seen  that  the  Spitsbergen  array  has  practically  the 
same  noise  level  as  ARCESS,  especially  for  the  “best”  Spitsbergen  sensor  (B2).  However, 
the  Spitsbergen  array  has  a  much  larger  variation  in  noise  level  than  ARCESS,  especially 
at  low  frequencies. 

Fig.  7.5.4a  shows  41  corrected  noise  spectra,  taken  daily  at  00.00  GMT  during  the  period 
10  December  1992  -  18  January  1993  for  the  Spitsbergen  array  (instrument  B2).  These 
spectra  represent  typical  winter  conditions.  For  comparison,  a  similar  set  of  spectra  during 
summer  (23  May  -  7  June)  is  shown  in  Fig.  7.5.4.b.  The  range  is  similar  for  die  two  time 
periods,  although  the  highest  noise  levels  at  low  frequencies  are  observed  during  winter. 

Noise  suppression 

Previous  studies  have  shown  that  regional  arrays  are  very  effective  in  suppressing  seismic 
noise,  thus  providing  gains  that  are  often  of  the  order  of  Jn,  or  sometimes  even  in  excess 
of  this  value  (N  being  the  number  of  sensors).  Such  better  than  Jn  suppression  occurs 
when  particular  subgeometries  are  chosen,  enhancing  the  suppression  of  noise  at  certain 
frequencies.  As  a  first  check  on  the  capabilities  of  the  Spitsbeigen  array  in  this  regard,  we 
have  computed  noise  suppression  curves  for  the  initial  array  geometry. 

We  calculate  an  uncorrected  power  density  spectrum  first  by  prewhitening  60  seconds  of 
data.  Then  we  estimate  the  autocorrelation  for  6  partially  overlapping  windows  (window 
length  12  seconds),  and  compute  power  density  spectra  from  the  average  autocorrelation, 
with  compensation  for  the  prewhitening  filter. 
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An  average  spectrum  is  obtained  by  averaging  the  individual  channel  spectra  for  the  array. 
The  averaging  is  done  after  a  logarithmic  transform  of  the  spectra,  and  the  standard  devia¬ 
tion  for  each  frequency  point  is  calculated.  Each  spectrum  is  pointwise  compared  to  the 
average  spectrum.  If  a  value  is  outside  1 .5  standard  deviation  from  the  mean  value,  the 
point  is  considered  an  outlier.  If  a  single  channel  spectrum  has  more  than  60%  outliers,  die 
spectrum  is  excluded  and  a  new  average  spectrum  is  estimated.  A  beam  is  calculated  using 
only  those  channels  for  which  the  spectrum  was  accepted.  The  suppression  is  then  the 
beam  spectrum  divided  by  the  average  spectrum. 

Fig.  7.5.5  shows  noise  suppression  in  the  frequency  range  0-20  Hz  for  an  infinite-velocity 
beam  (no  time  delays)  for  the  Spitsbergen  array.  The  figure  shows  an  average  of  24  noise 
samples,  taken  hourly  on  day  345, 1992.  We  note  that  the  noise  suppression  is  between  5 
and  10  dB  above  2  Hz.  This  is  about  the  theoretical  value  (8  dB  for  6  sensors).  We  antici¬ 
pate  further  investigations  into  this  subject  after  the  full  array  is  deployed. 

Detection  processing 

Since  the  array  was  installed,  the  Spitsbergen  data  have  been  subjected  to  continuous  on¬ 
line  detection  processing  (DP)  at  NORSAR  using  the  standard  small-array  detection  algo¬ 
rithm.  The  initial  beam  deployment  is  shown  in  Table  7.5.1. 

As  an  example  of  a  regional  seismic  event  recorded  and  processed  using  Apatity  data,  we 
show  here  results  for  a  low-magnitude  earthquake  (mb  ~  2.0)  west  of  the  array  (on  the 
North  Atlantic  Ridge). 

Fig.  7.5.6  shows  individual  recordings  for  this  event.  Note  the  large  variations  in  signal 
amplitudes  across  the  array.  In  particular,  the  B2  instrument  shows  excellent  signal-to- 
noise  ratio  for  this  event  as  well  as  for  other  events  we  have  processed.  The  reason  for  this 
large  amplitude  variation  is  as  yet  unknown.  There  are  two  distinct  P-onsets  and  one  clear 
S -onset  on  the  waveforms. 

Frequency-wavenumber  analysis 

Frequency-wavenumber  solutions  for  P  and  S  are  shown  in  Fig.  7.5.7.  In  spite  of  the  small 
aperture  of  the  array  and  the  incomplete  deployment,  the  peaks  in  the  F-K  diagram  are 
well-defined.  However,  there  is  a  tendency  to  side  lobes  in  the  plots. 

From  Fig.  7.5.7  we  note  that  the  phase  velocities  of  each  of  the  two  phases  are  reasonably 
consistent  with  the  phase  type,  the  S  phase  having  the  slowest  velocity.  The  estimated  azi¬ 
muths  are  quite  consistent  (256  and  260  degrees,  respectively). 

In  conclusion,  our  preliminary  analysis  indicates  that  the  Spitsbergen  array  will  be  an 
important  supplement  to  the  seismic  array  network  in  Northern  Europe.  Its  inclusion  into 
the  Intelligent  Monitoring  System  (IMS)  will  in  particular  serve  to  improve  the  location 
precision  and  source  characterization  of  the  large  number  of  events  in  the  Barents  Sea  and 
adjacent  regions.  The  excellent  recordings  of  this  array  from  the  small  Novaya  Zemlya 
event  on  31  Dec  1992  (see  subsection  7.4)  is  a  further  illustration  of  this  point.  Further 
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evaluation  of  the  anay  capabilities  will  be  made  after  the  9-element  array  has  been  com¬ 
pleted. 

J.  Fyen 
F.  Ringdal 
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BEAM 

Velocity 

Azimuth 

Filter 

band 

Threshold 

N 

Configuration 

SOU 

99999.9 

0.0 

0.5 

_ 

1.5 

4.4 

6 

AOB 

S021 

99999.9 

0.0 

1.0 

- 

3.0 

4.4 

6 

AOB 

S031 

99999.9 

0.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S032 

11.0 

30.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S033 

11.0 

90.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S034 

11.0 

150.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S035 

11.0 

210.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S036 

11.0 

270.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S037 

11.0 

330.0 

1.5 

- 

3.5 

4.4 

6 

AOB 

S038 

15.0 

88.0 

1.5 

- 

3.5 

3.9 

6 

AOB 

S039 

10.0 

95.0 

1.5 

- 

3.5 

3.9 

6 

AOB 

S041 

99999.9 

0.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S042 

10.2 

30.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S043 

10.2 

90.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S044 

10.2 

150.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S045 

10.2 

210.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S046 

10.2 

270.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S047 

10.20 

330.0 

2.0 

- 

4.0 

4.4 

9 

AOAB 

S048 

15.0 

88.0 

2.0 

- 

4.0 

3.9 

9 

AOAB 

S049 

10.0 

95.0 

2.0 

- 

4.0 

3.9 

9 

AOAB 

S051 

99999.9 

0.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S052 

8.9 

30.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S053 

8.9 

90.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S054 

8.9 

150.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S055 

8.9 

210.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S056 

8.9 

270.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S057 

8.9 

330.0 

2.5 

- 

4.5 

4.4 

9 

AOAB 

S058 

15.0 

88.0 

2.5 

- 

4.5 

3.9 

9 

AOAB 

S059 

10.0 

95.0 

2.5 

- 

4.5 

3.9 

9 

AOAB 

S061 

99999.9 

0.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S062 

10.5 

30.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S063 

10.5 

90.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S064 

10.5 

150.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S065 

10.5 

210.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S066 

10.5 

270.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S067 

10.5 

330.0 

3.0 

- 

5.0 

4.4 

9 

AOAB 

S068 

15.0 

88.0 

3.0 

- 

5.0 

3.9 

9 

AOAB 

S069 

10.0 

95.0 

3.0 

- 

5.0 

3.9 

9 

AOAB 

S071 

99999.9 

0.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

S072 

11.1 

30.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

S073 

11.1 

90.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

S074 

11.1 

150.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

S075 

11.1 

210.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

S076 

11.1 

270.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

Table  7.5.1.  Spitsbergen  beam  table,  valid  from  1992/328  (23  November  1992).  The  table 
shows  the  name  of  the  beam,  velocity  (km/sec),  azimuth  (degrees),  filter  band  (Hz), 
STA/LTA  threshold,  and  configuration.  The  configuration  is  described  with  number 
of  sensors  and  a  configuration  code.  Here,  AOAB  means  center  AO  SPZ  plus  A-ring 
plus  B-ring,  and  AOB  means  AO  SPZ  plus  B-ring.  SI01  -  SI06  are  incoherent  beams 
using  SPZ  channels  only.  (Page  1  of  2) 
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S077 

11.1 

330.0 

3.5 

- 

5.5 

4.4 

9 

AOAB 

S081 

99999.9 

0.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

SO  82 

9.5 

30.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

S083 

9.5 

90.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

S084 

9.5 

150.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

S085 

9.5 

210.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

S086 

9.5 

270.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

S087 

9.5 

330.0 

4.0 

- 

8.0 

4.4 

9 

AOAB 

S091 

99999.9 

0.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S092 

10.5 

30.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S093 

10.5 

90.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S094 

10.5 

150.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S095 

10.5 

210.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S096 

10.5 

270.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S097 

10.5 

330.0 

5.0 

- 

10.0 

4.9 

9 

AOAB 

S101 

99999.9 

0.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S102 

9.9 

30.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S103 

9.9 

90.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S104 

9.9 

150.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S105 

9.9 

210.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S106 

9.9 

270.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S107 

9.9 

330.0 

8.0 

- 

9.0 

4.9 

9 

AOAB 

S201 

99999.9 

0.0 

1.0 

- 

3.0 

4.0 

6 

AOB 

S207 

99999.9 

0.0 

8.0 

- 

9.0 

4.5 

6 

AOB 

S254 

99999.9 

0.0 

2.0 

- 

4.0 

4.0 

6 

AOB 

S282 

99999.9 

0.0 

4.0 

- 

8.0 

4.0 

6 

AOB 

S310 

99999.9 

0.0 

1.0 

- 

2.0 

2.5 

6 

AOB 

S312 

99999.9 

0.0 

2.0 

- 

4.0 

2.4 

6 

AOB 

SI01 

99999.9 

0.0 

0.5 

- 

1.5 

3.2 

6 

AOB 

SI02 

99999.9 

0.0 

1.0 

- 

2.0 

3.2 

6 

AOB 

SI03 

99999.9 

0.0 

1.5 

- 

2.5 

3.0 

6 

AOB 

SI04 

99999.9 

0.0 

2.0 

- 

4.0 

2.6 

6 

AOB 

SI05 

99999.9 

0.0 

3.5 

- 

5.5 

2.6 

6 

AOB 

SI06 

99999.9 

0.0 

5.0 

- 

10.0 

2.8 

6 

AOB 
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Fig-  73.2.  Configuration  of  the  new  Spitsbergen  small-aperture  array. 
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SPI  Single  channel  corrected  spectra. 
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Fig.  733a.  Noise  spectra  corrected  for  system  response  for  the  Spitsbergen  array  for  6 
vertical  channels.  The  spectra  are  based  on  one  minute  of  data  at  00.00  hours  GMT 
on  day  345, 1992.  The  power  density  is  in  nm2/Hz. 

ARC  Single  channel  corrected  spectra. 
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Fig.  733b.  Same  as  Fig.  7.5.1a,  but  for  ARCESS  data  taken  at  00.00  hours  GMT  on  day 


345, 1992. 


126 


NOMAR  3ci  Rip.  2-92/93 


May  1993 


SPI  hourly  average  of  1  single  channel  spectra. 
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Fig.  7.5.4a.  Corrected  noise  spectra  for  the  Spitsbergen  array  for  41  one-minute  intervals 
taken  daily  at  00.00  GMT  from  10  December  1992  - 19  January  1993.  Each  spec¬ 
trum  represents  data  from  the  B2  SPZ  seismometer. 

SPI  hourly  average  of  1  single  channel  spectra. 

17  somite  fiery  hour  Sundoy  23  Kay  -  Hondo)  07  Jun« 
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Fig.  7.5.4b.  Same  as  Fig.  7.5.4a  but  for  the  time  period  23  May  -  7  June  1993. 
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SPI  hourly  noise  suppression  spectra. 


1  somples  every  hour.  Thursday  10  December 


Frequency  Hz 


Fig.  7.5.5.  Spitsbergen  array  noise  suppression  by  beamforming  for  the  initial  geometry. 
The  plot  shows  an  average  of  24  curves.  To  produce  each  of  these  curves,  an  infinite- 
velocity  beam  is  formed  and  the  spectrum  for  this  beam  is  divided  by  the  average  of 
the  single  sensor  spectra.  The  24  curves  result  from  one  minute  of  data  taken  hourly 
between  00.00  and  23.00  hours  GMT  on  day  345, 1992.  The  horizontal  line  at  -8  dB 
represents  Jn  suppression  for  6  sensors. 
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a)  P  phase 


1993-198:07.36.34.000  +  3.00  SEC 

WIDE-BAND  SLOWNESS  ANALYSIS 
CONTOURS  IN  DB  BELOW  MAXIMUM  PEAK 


NQRSAR20/07/93  1  1  :41 :4Q 


Fig.  7.5.7.  Broadband  F-K  analysis  results  for  the  P  and  S  phases  of  the  event  shown  in 
Fig.  7.5.6.  The  figure  shows  a)  the  P  phase,  and  b)  the  S  phase.  (Page  1  of  2). 
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7.6  An  evaluation  of  the  performance  of  the  Intelligent  Monitoring 
,  System 

Introduction 

The  Intelligent  Monitoring  System  (IMS)  is  a  computer  hardware  and  software  system 
i  which  analyzes  data  from  a  network  of  seismic  stations  (arrays  and  three-component  sta¬ 

tions)  to  automatically  detect  and  locate  seismic  events.  The  events  automatically  declared 
by  IMS  are  reviewed  by  an  analyst  and  modified  as  appropriate . 

The  version  in  current  operation  at  NORSAR  is  described  in  Bache  et  al  (1993).  This  ver- 
t  sion  of  IMS  provides  for  joint  processing  of  data  from  the  NORESS  and  ARCESS  arrays 

in  Norway,  the  FTNESA  array  in  Finland,  the  GERESS  array  in  Germany,  and  the  new 
(fall  of  1992)  array  ne-  Apatity  on  the  Kola  peninsula  of  Russia  (see  Fig.  7.6.1).  In  addi¬ 
tion,  data  from  the  i.ray  installed  at  Spitsbergen  in  November  1992  (Mykkeltveit  et  al, 
1992)  are  made  available  to  IMS  for  interact  e  analysis;  i.e.,  detections  from  this  array 
are  not  yet  automatically  associated  by  IMS.  Based  on  these  data  sources,  IMS  automati¬ 
cally  produces  a  bulletin  of  local,  regional  and  teleseismic  events.  The  key  module  used 
by  IMS  in  doing  this  is  ESAL  (Expert  System  for  Association  and  Location).  ESAL  uses 
the  detection  information  resulting  from  the  signal  processing  (the  SigPro  module,  con¬ 
taining  automatic  signal  detection  and  characterization)  performed  separately  for  each 
array,  and  forms  and  locates  events  using  artificial  intelligence  technologies,  notably  in  the 
*  form  of  rule-based  reasoning. 

One  of  the  basic  objectives  of  IMS  is  to  provide  automatic  event  definitions  of  a  quality 
that  significantly  reduces  the  burden  on  the  seismic  analyst.  Operational  experience  with 
IMS  is  thus  taken  into  account  to  produce  enhancements  to  the  system  that  result  in  per- 
<  formance  improvements.  NORSAR  has  tried  to  assist  in  this  process  by  undertaking  eval¬ 

uations  of  IMS  performance  following  die  release  of  new  versions  by  SAIC,  the  system 
developer.  The  results  of  such  studies  are  discussed  with  SAIC,  and  taken  into  account  in 
the  next  set  of  enhancements  to  the  IMS. 

t  The  subject  of  this  study  is  the  outcome  of  another  such  evaluation  of  the  IMS.  The  perfor¬ 

mance  of  the  system  was  checked  carefully  during  a  recent  one-week  test  period,  and  the 
focus  was  on  the  performance  of  ESAL.  The  performance  statistics  presented  in  the  fol¬ 
lowing  thus  basically  reflect  the  quality  of  the  automatic  results  produced  by  ESAL,  and 
not  the  quality  of  the  final  event  definitions  after  analyst  intervention.  In  addition  to  this, 

(  we  also  point  to  some  potential  for  improvements  in  SigPro  that  would  enhance  the  over¬ 

all  performance  of  IMS. 

Evlauation  procedure 

l  The  evaluation  was  conducted  by  NORSAR  analyst  and  seismologist  staffs,  who  analyzed 

carefully  and  thoroughly  complete  IMS  data  for  the  7-day  period  26  April  -  2  May  1993. 
The  task  was  threefold,  namely,  (i)  to  characterize  the  automatic  IMS  event  definitions,  (ii) 
to  modify  these  as  deemed  appropriate,  and  (iii)  to  check  for  missed  events,  i.e.,  events 
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that  were  not  defined  by  IMS,  but  where  detections  were  available  from  SigPro  that  should 
allow  ESAL  to  associate  phases  and  form  events.  Task  (ii)  is  important  in  providing  infor¬ 
mation  on  the  performance  of  various  IMS  algorithms.  Examples  here  are  statistics  on 
retiming  of  phase  onsets  and  renaming  of  regional  phases.  Task  (iii)  above  was  performed 
through  various  means:  The  bulletin  of  the  GBF  (Generalized  Beamforming)  processing, 
which  uses  the  detection  information  produced  by  SigPro  to  locate  regional  events,  was 
checked  to  see  whether  it  contained  events  not  declared  by  IMS.  Likewise,  the  results 
from  the  single-array  “ep”-processing  performed  at  NORSAR  based  on  the  SigPro  output 
was  checked.  In  addition,  completely  independent  bulletins  like  the  NORSAR  array  bulle¬ 
tin  and  the  PDE  bulletin  of  the  USGS  were  checked,  basically  to  identify  teleseismic 
events  that  IMS  might  have  failed  to  catch,  even  if  appropriate  detections  were  available 
from  SigPro. 

IMS  automatically  produces  “event  plots”  like  the  one  shown  in  Fig.  7.6.2.  These  plots 
give  a  good  starting  point  for  the  review  process,  as  the  inferences  made  by  ESAL  in 
forming  and  locating  the  events  can  be  judged  from  these  plots.  The  actual  assessment  of 
all  events  reviewed  for  this  one-week  evaluation  period  was  made  on  the  basis  of  using  the 
interactive  analysis  tools  offered  by  the  Analyst  Review  Station. 

Each  event  automatically  declared  by  IMS  was  assigned  to  one  of  the  folllowing  five  cate¬ 
gories: 

Acceptable  events:  In  addition  to  events  accepted  without  any  modifications,  this  cate¬ 
gory  includes  events  where  the  anlayst  made  relatively  minor  changes.  The  character  of 
these  changes  was  such  that  the  original  IMS  location  was  not  strongly  affected,  i.e.,  the 
resulting  change  in  location  should  be  less  than  50-100  km  for  regional  events.  For  tele- 
seismic  events  the  requirements  on  “acceptable”  events  were  more  relaxed.  The  modifica¬ 
tions  in  this  category  usually  amounted  to  retiming  of  phases,  association  and 
disassociation  of  phases,  as  well  as  renaming  of  phases,  especially  for  single-array  events 
(e.g.,  from  Sn  to  Lg  or  vice  versa,  from  Pg  to  Px,  from  Rg  to  Lg,  etc.). 

Seriously  mislocated  events:  This  category  includes  events  that  are  real  events,  but 
where  the  event  location  is  too  far  away  from  the  true  location  to  qualify  for  the  “Accept¬ 
able  events”  category.  For  these  cases,  either  the  phase  assignments  were  wrong,  or  phases 
not  belonging  to  the  same  event  were  associated. 

False  events:  This  category  consists  of  those  events  declared  by  the  IMS  that  the  analyst 
rejected,  believing  they  were  not  real  events. 

Inconclusive  events:  For  this  category,  the  analyst  was  not  able  to  reach  a  definite  conclu¬ 
sion  whether  the  events  were  real  or  not. 

Missed  events:  This  category  includes  events  that  were  not  declared  by  IMS,  but  where 
detections  were  available  from  SigPro,  that  should  permit  ESAL  to  associate  phases  and 
form  events.  In  addition  to  such  events,  the  analyst  occasionally  produces  events  by  man¬ 
ually  adding  signal  arrivals  that  were  not  detected  by  SigPro.  An  example  here  would  be  a 
regional  event  for  which  SigPro  has  detected  a  P  phase  but  where  the  S  phase,  though  vis- 
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ible  to  the  analyst  during  interactive  waveform  analysis,  has  gone  undetected.  Such  events 
are  not  counted  as  missed  in  the  statistics  presented  in  this  report. 

All  events  in  these  five  categories  are  divided  into  regional  and  teleseismic  events,  based 
on  their  epicentral  distance  from  the  network: 

Regional  events:  Events  where  the  closest  station  is  within  20  degrees  of  the  epicenter. 

Deleseismic  events:  Events  where  the  distance  to  the  closest  detecting  station  exceeds  20 
degrees. 

Results  and  discussion 

Table  7.6.1  summarizes  the  results  of  the  characterization  of  the  events  automatically 
declared  by  IMS  during  the  one-week  test  period,  and  in  addition  includes  statistics  on 
events  missed  by  IMS.  As  seen  from  this  table,  nearly  80%  of  the  events  declared  are  con¬ 
sidered  to  be  acceptable,  whereas  the  majority  of  the  remaining  ones  are  considered  to  be 
false. 

All  events  in  the  category  “acceptable”  were  considered  by  the  analyst  to  be  real  seismic 
events  with  automatic  location  estimates  that  were  either  not  changed  during  the  subse¬ 
quent  interactive  analysis,  or  that  were  slightly  modified  during  this  process  through  minor 
changes  as  described  in  the  previous  paragraph.  Renaming  of  regional  phases  is  one  exam¬ 
ple  of  such  modifications,  and  Table  7.6.2  offers  statistics  in  this  regard.  Only  phases  for 
which  the  original  arrival  time  was  changed  by  less  than  2.0  seconds  are  included  in  this 
table.  This  is  done  to  exclude  cases  where  a  change  in  phase  name  was  accompanied  by  a 
substantial  change  in  the  phase  arrival  time.  It  is  noteworthy  that  only  for  one  single  case 
was  a  phase  changed  from  P-  to  S-type  by  the  analyst,  and  there  were  no  cases  where  an  S- 
type  phase  was  changed  to  P.  The  extent  of  renaming  can  thus  be  characterized  as  being 
relatively  modest,  meaning  that  the  automatic  phase  identification  in  ESAL  now  works 
very  well. 

As  seen  in  Table  7.6.1,  altogether  105  events  automatically  declared  by  IMS  were  rejected 
by  the  analyst  as  false,  based  on  various  kinds  of  evidence.  We  have  taken  a  closer  look  at 
these  events  to  see  if  there  is  parametric  or  other  kind  of  information  available  that  might 
permit  ESAL  to  automatically  reject  these  events.  The  following  observations  are  made: 

•  For  42  (all  regional)  of  these  105  events,  there  is,  in  the  judgement  of  the  analyst, 
parametric  information  available  that  might  be  used  by  ESAL  to  preclude  the  forma¬ 
tion  of  an  event.  Examples  of  such  parametric  evidence  are  high  or  low  Pn  veloci¬ 
ties,  and  frequency  content  of  Pn  and  Lg  phases  well  outside  the  expected  range  for 
these  phases.  The  formation  of  these  events  might  be  precluded  through  the  addition 
of  new  rules  in  ESAL  in  the  form  of  consistency  checks  on  the  parameters  pertain¬ 
ing  to  phases  used  in  forming  these  events. 

•  15  (also  all  regional)  of  the  false  events  were  rejected  by  the  analyst  because  they 
were  located  close  to  an  array  in  the  network  which  showed  no  sign  of  any  signal 
arrivals  from  this  event.  (It  was  checked  that  the  closest  array  was  operating  nor- 
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mally  in  these  cases.)  This  is  an  expample  where  contextual  information  could  be 
used  to  automatically  reject  an  event  hypothesis,  and  an  appropriate  consistency 
check  in  ESAL  might  rectify  the  situation.  For  several  of  the  events  in  this  category, 
also  parametric  information  indicated  that  the  events  were  false. 

•  Altogether  18  events  (both  regional  and  teleseismic  in  this  category)  were  judged  by 
the  analyst  to  be  so-called  “split”  events,  i.e.,  phase  arrivals  (often  coda  detections) 
belonging  to  a  real  event  were  used  to  define  an  additional,  false  event.  Sometimes 
the  arrival  azimuth  estimates  for  the  phases  used  to  define  the  split  event  deviated  by 
20°  or  so  from  those  of  the  phases  used  to  define  the  real  event,  and  it  would  be  dif¬ 
ficult  or  maybe  even  inadvisable  to  preclude  formation  of  the  second  event. 

We  find  that  27  of  the  false  events  originate  from  detections  resulting  from  bad  data 
(spikes  and  gaps  in  the  data).  For  S  of  the  events,  it  is  not  possible  for  the  analyst  to  see 
any  signal  at  all  (not  even  after  beamforming)  in  the  traces,  and  there  may  be  some  mal¬ 
function  of  the  detector  in  these  cases.  Possible  remedies  for  these  32  events  would 
amount  to  changes  in  SigPro,  but  it  should  be  noted  that  considerable  work  has  already 
been  invested  in  preventing  SigPro  from  declaring  detections  when  bad  data  are  recorded. 

There  are  a  few  false  regional  events  at  distances  around  10°  that  fall  into  none  of  the  cat¬ 
egories  dealt  with  above.  For  these  events,  the  only  reason  for  rejecting  them  is  the 
impression  left  after  close  inspection  of  the  waveforms  The  S  wavetrain  is  impulsive  and 
of  very  short  duration,  and  thus  does  not  match  the  expected  shape  (emergent  onset  and  a 
coda  of  some  length).  It  is  of  course  very  difficult  to  reject  such  events  automatically,  at 
least  until  the  AI  technology  is  able  to  match  the  trained  eye  of  a  skillful  analyst. 

As  seen  in  Table  7.6.1,  there  are  15  regional  events  that  were  missed  by  IMS  in  the  sense 
that  IMS  had  available  detection  information  that  should  allow  ESAL  to  associate  phases 
and  form  and  locate  events.  Only  one  of  these  events  had  detections  on  more  than  one 
array,  and  this  is  a  small  (local  magnitude  1.2)  Khibiny  Massif  event  recorded  at  ARCESS 
and  the  Apatity  array.  Three  of  the  missed  one-array  events  were  associated  with  double 
events  (two  mine  blasts  at  the  same  site,  10  seconds  or  so  apart),  where  IMS  only  defined 
one  event.  For  4  of  the  remaining  11  missed  one-array  events,  the  arrival  azimuth  esti¬ 
mates  for  P  and  S  phases  differed  by  more  than  20°  (they  were  in  the  range  22.8-25.9°). 
This  may  be  the  reason  why  ESAL  did  not  form  these  4  events,  but  it  appears  that  the 
remaining  7  events  (with  azimuth  differences  between  P  and  S  phases  of  less  than  16.4°) 
should  all  have  been  formed  by  ESAL. 

Table  7.6.1  shows  that  14  teleseismic  events  were  missed  by  IMS.  These  were  events 
defined  in  the  reference  bulletins,  and  for  which  ESAL  had  detections  available  (from  2  or 
more  arrays)  that  apparently  should  have  been  associated.  It  will  be  necessary  to  have  a 
closer  look  at  all  of  these  to  determine  why  they  were  not  formed  by  ESAL. 

Conclusions  and  recommendations 

Our  main  impression  after  having  carefully  analyzed  one  week  of  data  is  that  the  overall 
performance  of  IMS  is  now  very  satisfactory.  For  example,  the  rules  used  for  automatic 
phase  identification  in  ESAL  appear  to  work  very  well.  It  is  also  observed  that  problems 
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associated  with  earlier  versions  of  IMS  have  now  to  a  large  extent  been  solved.  There  is, 
however,  still  room  for  some  improvement,  as  discussed  above  and  summarized  in  the  fol¬ 
lowing. 

The  number  of  false  events  may  appear  a  bit  high,  and  we  have  suggested  some  possible 
remedies  that  might  help  the  IMS  to  automatically  reject  some  of  these.  Such  changes 
must,  however,  be  tested  very  carefully  to  make  sure  they  do  not  have  unintended  effects, 
like  throwing  out  real  events.  In  fact,  in  order  to  make  sure  that  the  system  does  not  miss 
real  events,  we  will  just  have  to  cope  with  a  certain  number  of  false  events.  This  number 
can  probably  be  reduced  quite  a  bit  from  today’s  about  20%  of  the  total  number  of  events 
automatically  declared  by  IMS,  and  it  must  remain  a  goal  for  IMS  to  minimize  the  burden 
on  the  analyst  by  defining  as  few  false  events  as  possible. 

Events  missed  by  IMS  represent  a  more  serious  problem  than  that  of  the  false  events.  Due 
to  the  large  amounts  of  data  processed  by  the  IMS  and  the  limited  manpower  resources 
available  for  interactive  analysis  of  the  automatic  IMS  results,  it  is  unlikely  that  the  ana¬ 
lyst  review  process  will  pick  up  the  events  missed  by  ESAL.  Keeping  in  mind  the  basic 
purpose  of  IMS,  it  is  therefore  of  utmost  importance  that  ESAL  captures  all  real  events, 
for  which  there  is  a  solid  basis  for  phase  association  and  event  formation.  We  have  seen  in 
the  aforegoing  that  some  of  the  missed  regional  events  were  associated  with  double  min¬ 
ing  blasts.  Although  this  is  a  situation  where  a  trained  analyst  could  pick  up  the  missed 
event  when  inspecting  the  waveform  traces,  it  will  certainly  present  a  challenge  to  capture 
all  such  events  and  at  the  same  time  avoid  formation  of  split  events.  All  except  one  of  the 
other  regional  events  missed  were  very  small  one-array  events.  Still,  it  is  necessary  to 
have  a  close  look  at  ESAL  to  rectify  this  problem,  as  well  as  the  problem  of  the  missed  tel- 
eseismic  events. 

As  we  have  already  seen,  IMS  performance  could  also  be  enhanced  through  certain  modi¬ 
fications  to  SigPro.  We  will  here  touch  upon  another  aspect  where  changes  to  SigPro 
might  be  beneficial:  Figs.  7.6.3  and  7.6.4  provide  statistics  on  retiming  of  the  phases  Pn 
and  Pg,  respectively,  during  the  course  of  analyst  review  of  the  data  for  the  one-week  eval¬ 
uation  period.  The  figures  show  the  differences  between  the  arrival  time  as  automatically 
determined  by  SigPro  and  the  arrival  time  as  determined  by  the  analyst  during  the  review 
process,  plotted  versus  the  SNR  of  the  Pn  or  Pg  phase,  calculated  from  the  detecting  beam. 
Only  phases  that  after  analyst  review  retained  their  original,  automatic  phase  assignment 
(Pn  or  Pg)  and  in  addition  correspond  to  real  seismic  events,  are  included  in  Figs.  7.6.3 
and  7.6.4.  The  figures  show  standard  deviations  of  the  order  of  half  a  second  for  the  differ¬ 
ences  in  the  arrival  time  estimates,  and  there  are  appreciable  differences  even  for  high 
SNR  phases.  This  indicates  that  there  should  be  quite  some  potential  for  improvement  in 
the  automatic  estimation  of  phase  arrival  times.  Section  7.2  of  this  Semiannual  Technical 
Summary  presents  an  approach  that  might  be  implemented  in  SigPro  and  that  holds  con¬ 
siderable  promise  to  improve  the  automatic  onset  times. 
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The  current  version  of  IMS  makes  use  of  some  region-specific  knowledge.  Further 
enhancements  to  IMS  performance  are  likely  to  be  obtained  through  introduction  of  addi¬ 
tional  such  knowledge,  and  section  7.3  of  this  Semiannual  Technical  Summary  demon¬ 
strates  how  event  locations  may  be  significantly  improved  using  region-specific 
knowledge. 


S.  MykkeUveit  U.  Baadshaug  B.Kr.  Hokland 

T.  Kvserna  L.B.  Lough  ran 
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Acceptable  events 

384 

67 

451 

Seriously  mislocated  events 

5 

4 

9 

False  events 

97 

8 

105 

Inconclusive 

0 

4 

4 

Total  number  of  events  declared 

569 

Missed  events 

15 

14 

29 
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Table  7.6.1  Characterization  of  all  events  automatically  declared  by  IMS  during  the 

seven-day  period  26  April  -  2  May  1993  (see  text  for  explanation  of  the  various  cat¬ 
egories).  Tlie  table  also  includes  statistics  on  events  missed  by  the  IMS. 
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Pn 

Pg 

Px 

Sn 

Lg 

Rg 

Sx 

Pn 

213 

19 

2 

Pg 

11 

203 

24 

1 

Px 

2 

2 

60 

Sn 

34 

5 

9 

Lg 

7 

325 

7 

24 

Rg 

6 

30 

7 

Sx 

1 

20 

75 
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Table  7.6.2  The  table  shows  the  ESAL  automatic  phase  assignments  versus  analyst 

assignments  made  during  interactive  analysis,  for  all  regional  phases  that  were  used 
by  IMS  to  define  regional  events  during  the  time  period  26  April  -  2  May  1993. 
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Fig.  7.63.  The  figure  shows  differences  between  arrival  times  as  determined  by  SigPro 
and  arrival  times  determined  by  the  analyst  during  review,  plotted  versus  SNR,  for 
Pn  phases  used  by  IMS  to  form  events  during  the  evaluation  period  26  April  -  2  May 
1993.  The  mean  value  of  the  arrival  time  differences  is  marked  by  a  solid  vertical 
line,  whereas  the  two  strippled  lines  denote  ±  2  standard  deviations. 
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Fig.  7.6.4.  Same  as  Fig.  7.6.3,  but  for  the  Pg  phase. 


